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1.1. Sistemas, elementos e relações 


Sistema é um conceito que povoa densamente o nosso quotidiano e o nosso imaginário. 
Falamos do “sistema” a propósito de praticamente todos os aspetos da nossa vida e da 


natureza. Desculpabilizamo-nos dizendo que “a culpa é do sistema”. 
Mas afinal o que é um sistema ? 


Poderemos tentativamente defini-lo como um conjunto de elementos, ou componentes, 


relacionados e interagindo entre si de modo a desempenharem alguma função. 


Genericamente diz-se que existe uma relação entre os elementos El e E2 (de um sistema) se 
o comportamento de El é influenciado ou controlado pelo de E2 e se o comportamento de E2 


é influenciado ou controlado pelo de El. 


Um atributo de um elemento é uma propriedade característica do elemento (cor, tamanho, 
peso, etc.). Um atributo de uma relação é uma qualidade característica da relação 
(intensidade, velocidade de comunicação, etc.). No estudo de sistemas dinâmicos a evolução 
dos atributos assume um espeto fundamental, dado que é essa evolução expressa por 
mudanças várias que exprime a dinâmica. Dinâmica e mudança são conceitos estreitamente 


ligados. 


Se entendermos por elemento algum objeto ou algum fenómeno da vida natural ou social, ou 
a sua representação, que tem certos atributos que podem variar ao longo do tempo em 
consequência do comportamento do próprio elemento, teremos uma definição suficientemente 
geral para abranger o objetivo do nosso estudo: os sistemas biológicos e fisiológicos. Um 


sistema diz-se dinâmico se os atributos dos seus elementos variam ao longo do tempo. 


Há dois conceitos-chave nesta definição de sistema: relação e interação. Sem eles não existe 


o sistema. 


A Figura 1 ilustra a importância desses conceitos. Na parte a) temos três elementos: um vaso 
com água, uma lamparina e uma caixa de fósforos, lado a lado sobre a mesa. Não há qualquer 
relação entre eles. Não constituem um sistema. Na parte b), acende-se a lamparina com um 
fósforo, coloca-se a lamparina sob o vaso com água (usando um suporte), e agora os 


elementos estão relacionados e interagem: a chama da lamparina transmite energia térmica à 
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água através das paredes do vaso. Temos um sistema com um objetivo: o aquecimento da 


água do vaso. 


Lo = 


Vaso com água Lamparina Fósforos 


a) b) 


Figura 1. a) O vaso, a lamparina e os fósforos não constituem um sistema. b) montando-os de forma 
conveniente, estabelecendo entre eles relações apropriadas, compõem um sistema de aquecimento de 


água. 


Ao estabelecermos relações entre os elementos do sistema, surgiram propriedades novas, que 
não existiam em qualquer elemento isolado, mas que nascem com as inter-relações entre os 


elementos. Chama-se efeito de sinergia: o todo é maior do que a soma das partes. 


Suponhamos agora um computador em bom estado de funcionamento. Todas as suas partes 
interagem entre si. Entre a CPU e o disco duro circula informação. Entre a fonte de 
alimentação e o processador, ou a memória, circula energia. Se desmontarmos as partes e as 
amontoarmos na nossa secretária, quebrando a interação entre os componentes, deixa de haver 
circulação de informação e de energia e deixa por isso de haver sistema. O sistema não é igual 


à soma dos seus componentes. Ele é mais do que isso. 


Um automóvel em funcionamento contém um número muito elevado de componentes de 
diversa natureza: depósito de combustível, motor, pneus, eletrónica e microprocessadores, etc. 
Entre eles circula informação de múltiplos sensores e atuadores para os microprocessadores e 
vice-versa, circula massa (do depósito de combustível para os injetores no motor) e energia 


(do motor para as rodas, por exemplo, energia mecânica). 
Massa, energia e informação são os três elementos circulantes dos sistemas modernos. 


Uma célula animal ou vegetal é uma fábrica sofisticada e complexa. Cada um dos seus 
componentes tem uma missão específica, e a interação entre eles, conjugando as suas 


missões, permite a vida. Também aí fenómenos de trocas de massa (nutrientes, por exemplo), 
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de energia (por exemplo calor), de informação (através dos genes e das proteínas) são a 


condição primitiva para a existência do sistema celular no seu todo. 


ion of an Animai Celi” g 


att 
Membr 
Membrane, 


Vacuole 


| | 
RO 9, 
| 


a) Célula animal b) Célula vegetal 


Figura 2. Uma célula é um sistema já muito elaborado e complexo (de 


(http://www enchantedlcarning.com/subjects/biology/cellsy, 12 julho 2022). 


1.2. Retroação (Feedback) e (auto)regulação 


A Figura 3 representa uma instalação de ar condicionado numa sala com regulação 
computorizada A temperatura da sala é medida por um termómetro (por exemplo um 
termopar) e enviada ao computador. Este compara com a temperatura desejada para a sala 
(chamada temperatura de referência) e conforme o erro observado manda aquecer ou arrefecer 


até que o erro se anule. O dispositivo que dá ou retira calor da sala chama-se atuador. 


Temperatura 
desejada 


Temperatura 
medida 


Realimentação negativa (Feedback) 


Figura 3. Servomecanismo de regulação: ar condicionado automático. 


Trata-se de um sistema de regulação construído pelo homem com um objetivo preciso. 
Tradicionalmente chama-se um servomecanismo (de regulação), na medida em que há um 
mecanismo servo (escravo), neste caso o atuador. O sistema tem alguma capacidade de 
autorregulação, depois de lhe ser indicada a temperatura de referência. O elemento essencial 


que permite tal capacidade é a informação da temperatura da sala, através do caminho de 
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retroação (da saída para a entrada) ou de realimentação, em inglês feedback. Feedback 
permite a autorregulação. Suponhamos que num dia muito frio abrimos uma janela da sala e 
entra uma grande quantidade de ar frio. Esse facto constitui uma perturbação do ambiente do 
sistema. Com realimentação, o sistema é capaz de reagir a essa perturbação de modo a repor 


as condições de equilíbrio. Mas sem a realimentação tal não acontece. 


Na nossa vida quotidiana contactamos com centenas de servomecanismos variados (muitas 
vezes sem nos apercebermos). Pode-se mesmo afirmar que a nossa civilização tecnológica 
assenta em servomecanismos de (auto)regulação. Os sistemas biológicos têm múltiplos e 


intrincados caminhos de feedback que viabilizam a vida. 


Consideremos por exemplo um sistema biológico de duas populações, o par predador-presa, 
sendo a presa herbívoro e o predador carnívoro. Se a vegetação aumenta, a população de 
herbívoros cresce exponencialmente. Por isso a população de carnívoros aumenta também 
exponencialmente, o que provoca a diminuição dos herbívoros. Assim, o aumento inicial dos 


herbívoros retroagiu consigo própria através da sua relação com a população de carnívoros. 


A regulação do açúcar no sangue pode ser representada pela Figura 4 onde se evidenciam 
múltiplas malhas de realimentação que permitem a capacidade de autorregulação de um corpo 


saudável. 
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Figura 4. Sistema com várias malhas de retroação: a regulação do açúcar no sangue (de Bertalanffy, 1969). 
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Nos exemplos anteriores a realimentação é usada com sinal negativo, e por isso diz-se 
realimentação negativa (negative feedback). Mas pode existir com sinal positivo, e nesse caso 


diz-se realimentação positiva. 


A realimentação positiva pode ter efeitos favoráveis, mas também pode ter efeitos 
catastróficos sobre um sistema. Quando uma pessoa corre, necessita de mais oxigénio e por 
isso aumenta o ritmo da respiração (feedback positivo, mais corre mais respira). As malhas de 
realimentação positiva dominam durante esta fase, provocando um aumento da oxigenação do 
sangue, através dos pulmões, que por sua vez permite a corrida. Numa pessoa saudável, os 
limites da capacidade humana são determinados por feedback negativo, e por isso só se pode 
correr até um certo limite, o que nos impede de exaurir. Se as malhas de realimentação 
negativa se quebram, permitindo o domínio indesejável das positivas, como acontece quando 
os atletas tomas certas drogas dopantes, podem-se alcançar desempenhos sobre-humanos, mas 
por vezes com resultados trágicos, quando os processos biológicos dos atletas não conseguem 
aguentar requisitos descontroladas que levam ao colapso e à morte, como se tem já verificado 


em plenas provas. 


Um exemplo muito mais banal de realimentação positiva acontece quando, num palco, o 
microfone do interveniente está defronte de uma coluna de som: o som vai do microfone ao 
amplificador, deste ao altifalante, deste ao microfone, ao amplificador, etc., até que o 
amplificador satura e se ouve um som agudo muito desagradável. Por isso é necessário que o 
microfone esteja sempre fora do campo de abrangência do altifalante, ou seja, atrás do seu 


plano frontal. 


Uma relação entre elementos do sistema materializa-se, com vimos, através de fluxos de 
materiais, de informação ou de energia entre os elementos. Por isso o conceito de 


comunicação pode ser usado equivalentemente ao de relação. 


A definição que adotámos de sistema, é muito geral, permitindo englobar fenômenos da vida 
natural desde a célula até ao próprio universo. Por isso, em cada caso concreto, temos que nos 
contextualizar distinguindo o que é para nós o sistema em estudo do que o rodeia (ou do seu 


ambiente). 
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1.3 Ambiente, sistemas fechados, sistemas abertos e homeostasia 


Um sistema está sempre envolvido pelo seu ambiente e com ele comunica através da sua 


fronteira que o demarca do resto do mundo. 


Não é sistema Não é sistema 


Figura 5 a). Um conjunto de elementos não compõem um sistema. b) algumas relações entre alguns 
elementos não são suficientes para se considerar um sistema. c) se houver uma densidade de relações 
envolvendo a maior parte dos elementos (ainda que um ou outro se possam manter isolados) teremos 


um sistema com entradas, saídas e uma fronteira com o meio ambiente (de Flood e Carson (1993)). 


A intensidade da comunicação entre o sistema e o seu ambiente é variável. Se é fraca, diz-se 
que o sistema é fechado. Uma reação química num reator fechado em que se misturam os 
reagentes, é descrita, do ponto de vista do exterior, por um sistema fechado. E a este tipo de 


sistemas que se aplicam as leis da termodinâmica. 


Mas se a comunicação entre um sistema e o seu ambiente é intensa, estamos em presença de 
um sistema aberto, que troca materiais, informação e energia com o seu ambiente, através de 


uma fronteira. Uma célula viva é um caso típico de sistema aberto. 


Num sistema aberto, dado que as trocas são permanentes e contínuas, nem sempre é possível 


identificar a fronteira com precisão. 


Os seres vivos são sistemas abertos: estão sempre a receber do ambiente e a enviar para o 
ambiente materiais e energia (e também informação), a sua estrutura interna está em 
permanente mudança, e por isso nunca se encontram em estado estacionário (no sentido de 
que todos os atributos de todos os seus elementos têm um valor constante ao longo do tempo), 
mas em equilíbrio químico e termodinâmico (em que os atributos apresentam padrões 
definidos ao longo do tempo, com possíveis variações dentro de certos limites). A este 
equilíbrio chama-se homeostasia ou homeostase .O conceito foi criado pelo fisiologista 


francês Claude Bernard (1813-1878), https://en.wikipedia.org/wiki/Claude Bernard, com a sua 
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famosa frase "La fixité du milieu intérieur est la condition de la vie libre" que é ainda hoje o 
princípio subjacente à homeostasia.- a constância do meio interior é a condição da vida. O 
termo foi depois fixado pelo fisiologista americano Walter Cannon (1871-1945), no seu livro 
The Wisdom of the Body (1932), onde descreve as características gerais da homeostasia 


através de quatro proposições, http://en.wikipedia.org/wiki/Walter Cannon: 


I- A constância num sistema aberto, como o representado pelo nosso corpo, requer um 
mecanismo que atue para manter essa constância (como por exemplo na concentração da 
glucose, na temperatura corporal, no balanço ácido-base, etc.). Constância deve interpretar- 
se aqui como estando dentro de um intervalo de tolerância saudável, e não como um valor 
absolutamente constante ao longo do tempo; por exemplo a glucose estar a 120mg/l ou 


123mg/I não faz diferença do ponto de vista do estado de saúde da pessoa. 


2- As condições de regime estacionário, no sentido de que as variáveis estão dentro do 
intervalo de tolerância, requerem que qualquer tendência para a mudança (saída do 
intervalo de tolerância) se confronte com fatores que resistam à mudança. Por exemplo um 
aumento da concentração de açúcar no sangue produz a sensação de sede — o corpo procura 


diluir a concentração de açúcar no fluido extracelular. 


3- O sistema de regulação que determina o estado homeostático é composto por vários 
mecanismos cooperantes atuando simultânea ou sucessivamente. Por exemplo o açúcar no 
sangue é regulado pela insulina (produzido pelas células beta do pâncreas), glucagon ou 
glicagina (produzida pelas células alfa do pâncreas, funciona como uma hormona oposta à 
insulina), e outras hormonas que controlam a libertação de insulina pelo pâncreas e a sua 
absorção pelos tecidos. Para regulação da temperatura no corpo humano, o arrefecimento 
do sangue (pelo contacto com o ambiente) estimula o cérebro que ativa mecanismos de 
produção de calor que consomem calorias da alimentação, o que faz com que a 


temperatura do corpo recupere e suba. 


4- A homeostasia não acontece por acaso, mas é o resultado do autogoverno organizado. 


KEY « |molecule 
a-o|labeis on molecules R e 
[F Jtime E o do 


iai celi at 1 - (b) cell at t+s 


Figura 6. Ilustração da homeostasia. Os componentes interiores da célula foram todos substituídos entre os instantes 
te t+s, mas a célula manteve-se num estado estacionário porque a sua estrutura interna manteve-se como 


anteriormente, embora com novos elementos (de Flood e Carson (1993)). 
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1.4 Entropia e neguentropia 


A entropia exprime a tendência que as coisas têm para evoluírem para um estado de grande 
desordem e desorganização. Ela carrega consigo a lei irreversível da degradação da energia e 
da matéria, tal como enunciado pelo segundo princípio da termodinâmica: os sistemas tendem 


para a desordem. Essa foi a conclusão de  Clausius | (https://mathshistory.st- 


andrews.ac.uk/Biographies/Clausius/ 12 julho 2022) ao enunciar o seu famoso princípio da máxima 


entropia (“ a entropia do universo tende para um máximo” (cresce sempre)), ou, 


matematicamente teremos (1) em que S representa a entropia e dS a sua variação. 
dS>0 (1) 


Mas os sistemas homeostáticos mantêm a ordem, tendem para a ordem. Não haverá aqui uma 


contradição na ciência ? 


Foi este o problema tratado por Ilya Prigogine, propondo em 1945 o teorema da mínima 
produção de entropia aplicável a estados estacionários de não equilíbrio, dando uma 
explicação da analogia que relaciona a estabilidade de estados de equilíbrio termodinâmico 
(noção clássica) com a estabilidade de sistemas biológicos, tais como o expresso pelo 
conceito de homeostasia. Prigogine foi Prémio Nobel da Química em 1977 


(https:/Awww nobelprize.org/prizes/chemistry/1977/summary/ 12 julho 2022). Este resultado é dado pela 


equação de Prigogine (2) e (3), 


dS=d,S+d,S (2) 


dS <0 (3) 


em que diS representa a entropia produzida internamente pelo processo irreversível de 
degradação (reações químicas, difusão, transporte de calor, etc.) e deS representa a entropia 


importada pelo sistema vivo, aberto, sob a forma de matéria e energia. Pela lei de Clausius, 
diS > 0 (4) 
substituindo (4) em (2) facilmente se conclui (5) 


deS<0 . (5) 
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Devido a (5) chama-se neguentropia, entropia negativa a deS. Os sistemas abertos importam 


neguentropia. 


ash 


o 


a) b) 
Figura 7. a) Sistema fechado (equação de Clausius) e b) sistema aberto (equação de Prigogine) 


Por isso as leis da termodinâmica aplicam-se apenas aos sistemas fechados. A entropia 
aparece assim como uma força em luta contra a homeostasia. Com a homeostasia, os sistemas 
importam materiais, energia e informação do ambiente, e com eles compensam a tendência 
para a desordem, criando e mantendo a ordem, resistindo ao segundo princípio da 


termodinâmica, e podendo até evoluir para estados de maior ordem e maior organização. 


A estrutura de um sistema exprime a forma com os seus elementos se relacionam uns com os 


outros e estabelece a base de suporte aos processos que nele decorrem. 


Se observarmos um sistema em instantes sucessivos, aperceber-nos-emos que ele evolui de 
alguma forma. O comportamento do sistema é essa nossa perceção. O comportamento do 
sistema diz-se orientado por um objetivo se puder ser apreciado e avaliado relativamente a 


uma finalidade particular. 


Nos sistemas abertos um dado estado final de equilíbrio pode ser alcançado a partir de 
diferentes estados iniciais, seguindo diferentes trajetórias. Diz-se, neste sentido, que o estado 


de equilíbrio tem um valor equifinal e que o sistema goza da propriedade de equifinalidade 


(equifinality) 
1.5 Estado e espaço de estados 


Posteriormente, no Capítulo 4, iremos estudar a noção de estado e de espaço de estados e as 
suas propriedades. Podemos, no entanto, desde já ver uma noção inicial de variáveis de 


estado. Alguns atributos elementares de um sistema são muito importantes e de per si 
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permitem-nos ter uma imagem suficientemente precisa do sistema num dado momento, isto é, 
do seu estado. Por exemplo a temperatura de uma estufa, a velocidade de um veículo, a 
população de uma espécie, o nível de um reservatório de água, etc. Esses atributos chamam-se 
por isso variáveis de estado. Para a representação matemática define-se um vetor x, chamado 


vetor de estado, cujos componentes x; são as variáveis de estado. 


fla (6) 


Como veremos, o estado é uma síntese da história do sistema: o estado atual do sistema 
resulta de toda a sua evolução passada. Num sistema dinâmico o estado variará ao longo do 
tempo. Se o representarmos graficamente ele descreve uma trajetória no espaço n- 
dimensional, chamado espaço de estados. Se pudermos prever com rigor os valores futuros do 
estado, conhecido exatamente o seu estado inicial, estamos em presença de um sistema 
determinístico; caso contrário temos um sistema não-determinístico, ou probabilístico, ou 


estocástico. 
1.6 Autopoiese, adaptação, cibernética 


Um sistema que goza da propriedade de autopoiese tem a capacidade de se reproduzir: uma 
célula produz os seus próprios componentes que por sua vez a produzem. Autopoiese 
significa por isso a capacidade de auto-produção, o que acontece com os seres vivos: eles 
são autopoiéticos porque estão organizados para permitirem que os seus próprios processos 


produzam os componentes que são necessários para a continuação desses mesmos processos. 


Um sistema diz-se adaptativo se tiver capacidade de adaptação. A adaptação surge como a 
necessidade de resposta a mudanças no ambiente com vista à sobrevivência e à manutenção 
de um bom desempenho do sistema. Pode-se dizer que a teoria darwiniana da evolução 
constitui uma teoria de adaptação. A capacidade de adaptação é decisiva quando o ambiente 
muda com intensidade e com frequência. A adaptação requer a capacidade de regulação e 
controlo para resposta às mudanças bruscas e de longo prazo (regulação de curto prazo e 


regulação de longo prazo). Qualquer ser vivo requer múltiplos sistemas de regulação e de 
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controlo de curto e médio prazo para sobreviver; a isto se chama por vezes a lei da variedade 


requerida (requisite variety). 


Por vezes na nossa vida olhamos para as coisas apenas parcialmente, vendo apenas algumas 
partes e esquecendo-nos de olharmos para elas globalmente, holisticamente. Por isso o 
resultado das nossas ações é muitas vezes contra as nossas expectativas, por termos 
negligenciado ou esquecido a complexidade dos fenómenos. Este facto traz-nos a necessidade 
do pensamento sistémico, chama a atenção para a necessidade de metodologias e de modelos 
para lidar com a complexidade. Sem este pensamento formal veremos apenas partes, as 


explicações simplistas, as soluções simplistas. 


A cibernética estuda a adaptação, a regulação e o controlo. Podemos definir cibernética como 
a ciência do controlo e da comunicação nos animais e nas máquinas (proposta por Wiener em 
1948 no livro fundador Cybernetics, or Control and Communication in the Animal and 
Machine; descreve as leis naturais que governam a comunicação e o controlo em sistemas 


dinâmicos. 
1.7. Uma hierarquia de sistemas 


Numa hierarquia de regulação e de controlo poderemos definir sistemas em diversos níveis: 
teremos assim uma hierarquia de sistemas. Uma meta sistema, ou hipersistema, é um 
sistema que se posiciona acima de um sistema numa dada hierarquia. Um subsistema 
posiciona-se abaixo. No ser humano o cérebro/consciência é o meta sistema de nível mais 
elevado e o átomo o de nível mais baixo. Teremos a hierarquia: átomo, molécula, célula, 


tecido, órgão, sistema funcional, ..., consciência. 


A organização hierárquica é assim uma representação lógica dos fenômenos como sistemas 
e subsistemas. Descendo a hierarquia, aumentamos o nível de resolução da análise dos 
fenómenos, vendo cada vez mais detalhes. A escolha de um nível apropriado de resolução 
para análise de um fenómeno é uma tarefa decisiva que determina no fundo quais as questões 
com que lidaremos no estudo. O nível de resolução escolhido chama-se o sistema em 


enfoque (system-in-focus). 


Para se ser um bom cientista de sistemas tem que se ser simultaneamente holístico (olhar para 
o sistema como um todo) e reducionista (ser capaz de compreender o sistema de modo mais 


detalhado, em cada uma das suas partes). 
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Quando se sobe numa hierarquia de sistemas descobre-se uma propriedade fundamental de 
todo a teoria dos sistemas: o todo é maior do que a soma das partes. Trata-se da 
propriedade de emergência: há propriedades do hipersistema que emergem apenas (por vezes 


inesperadamente) quando se combinam os subsistemas. 


As células agrupam-se em diversos órgãos (rim, fígado, coração, pulmão, etc.) cada qual com 
uma função específica e com propriedades diferentes das suas células. Os órgãos, em 
conjunto, formam um corpo com propriedades emergentes: eles organizam-se através da 
comunicação e do controlo (cibernética) numa hierarquia de partes corporais que permitem o 
surgimento das capacidades de observação, audição, olfato, locomoção, emoção, isto é, da 
pessoa humana. Um ser humano não é um agregado das partes corporais. Nem uma sociedade 
um agregado de grupos humanos. Os sistemas juntam-se para formar hipersistemas cujas 


propriedades são diferentes das dos sistemas. 


A emergência é assim uma caracterização dos fenómenos que não conseguimos explicar de 


outro modo. 


A emergência de uma sociedade é a sua cultura. 

A sinergia é um termo usado para descrever a emergência de propriedades inesperadas e 
interessantes, por exemplo em teoria das organizações (a sinergia dos grupos de trabalho, por 
exemplo). 

De acordo com Bertalanffy, poderemos definir a hierarquia de sistemas da Tabela 1. 

Ao longo do último século surgiram diversas metáforas para os sistemas, à medida que se 
procurava aplicar a teoria dos sistemas ao funcionamento dos grupos humanos e das 


sociedades. 


Primeiramente a metáfora da máquina, olhando para um sistema fechado com um conjunto 
de objetivos bem precisos e uma estrutura rígida de controlo para os alcançar (os 
servomecanismos tecnológicos são disso um exemplo). Foi esta metáfora que inspirou o 


desenvolvimento das cadeias de montagem na indústria e a sua gestão. 


Depois de Prigogine essa visão foi ultrapassada, passando-se à metáfora orgânica que 
descreve redes complexas de elementos e de relações com uma transformação afetada por 
realimentação. 

Mas a metáfora orgânica não permite descrever as emoções. Surgiu por isso a metáfora do 


cérebro. 
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Nível Descrição e exemplos 


Sistemas simbólicos Linguagem, lógica, matemática, ciência, morais, etc. 


Populações de organismos (incluindo humanos); 
Sistemas socioculturais comunidades humanas determinadas por símbolos 


(culturas). 


Simbolismo; passado e futuro, o “eu” e o mundo; 
O Homem autoconsciência resultando como consequência a 


comunicação pela linguagem. 


Aumento da importância do tráfego da informação 
Animais (evolução dos recetores, sistemas nervosos, etc.); 


aprendizagem; início da consciência. 


Organismos do tipo planta: aumento da diferenciação 
Organismos inferiores nos sistemas (“divisão do trabalho” nos organismos); 


distinção da reprodução das funcionalidades individuais. 


Sistemas abertos Chama, células, tecidos, órgãos e organismos em geral. 


. Termostato, servomecanismos, mecanismos 
Mecanismos de controlo na | Ê 
homeostáticos em organismos vivos. 


Tabela 1. A hierarquia de sistemas segundo Bertalanffy (1969). 


A metáfora da cultura permite incluir uma rede de valores, crenças e normas. 
Finalmente a metáfora política permite redes de interesses interatuantes, interesses que os 


indivíduos procuram alcançar. 


No nosso estudo da dinâmica de processos biológicos e fisiológicos preocupar-nos-emos com 
a obtenção de modelos matemáticos computáveis (isto é, programáveis em computador e úteis 
para a análise, simulação e predição) de sistemas dos dois últimos níveis da hierarquia de 


Bertalanffy (mecanismos de controlo e sistemas abertos). 


Basearemos a nossa análise nas propriedades dos fenómenos físicos, químicos, 
termodinâmicos, etc. A aplicação das respetivas leis leva-nos a um conjunto de equações 


diferenciais que permitem descrever e simular computacionalmente o respetivo 
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comportamento dinâmico. As equações diferenciais servem também para a obtenção de 
formas de representação matemática muito práticas e de grande utilidade: a função de 


transferência (no caso linear) e as equações de estado (casos linear e não linear). 


Algumas propriedades dinâmicas, como a estabilidade e tipo de resposta transitória, são 


facilmente obtidas nessas representações especiais. 


Alguns sistemas não lineares têm uma importância especial no contexto biológico e 
fisiológico: os sistemas caóticos. Com veremos, um sistema é caótico se tiver um 
comportamento que não se possa prever com o conhecimento possível do seu estado atual. 


Veremos alguns exemplos concretos e uma introdução à sua teoria. 
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2.1 Introdução 


Como vimos no Capítulo 1, energia, massa e informação são os três elementos circulantes dos 
sistemas. Há, como sabemos, muitas formas de energia, e, portanto, sistemas de natureza 
muito diversa. Se quisermos definir uma taxonomia de sistemas, que nos apoie no 
estabelecimento de uma teoria unificadora, poderemos usar o critério da energia e classificá- 
los conforme o tipo de energia que neles circula. Teremos assim sistemas elétricos, 


mecânicos, térmicos, fluídicos, químicos, etc. 


Há sistemas com mais do que um tipo de energia: um robô tem energia mecânica 
(movimento) e energia elétrica (motores) sendo por isso um sistema eletromecânico; um 
acumulador doméstico de água quente é um sistema termoelétrico, numa instalação química 


existem energia química (reação), fluídica (circulação de líquidos e gases), térmica. 


Os sistemas biológicos e fisiológicos são ainda mais complexos. Eles contêm subsistemas de 
diversos tipos de energia, mas têm propriedades emergentes em relação a eles - as que criam a 


própria vida. 


Para analisar sistemas parcelares do corpo humano, por exemplo o cardiovascular ou o 
respiratório, mantendo-nos em níveis de complexidade tratáveis, teremos que lhes aplicar 
métodos de análise fundamentados nos princípios gerais dos sistemas do mesmo tipo de 
energia. O sistema cardiovascular é essencialmente um sistema fluídico: uma bomba, um 
conjunto de canais de circulação, uma certa resistência à circulação. O movimento muscular 


pode ser modelizado por um sistema mecânico. 


Veremos depois que é possível definir um conjunto de características básicas dos sistemas que 
atravessam todos os tipos, permitindo estabelecer um conjunto de analogias entre eles. Existe 
uma analogia entre dois sistemas se as suas equações matemáticas têm a mesma forma. Estas 
analogias permitem-nos usar um dado tipo de sistema (elétrico, por exemplo), para estudar 


outros. 
Na literatura biomédica é muito frequente encontrar-se esta abordagem. A figura seguinte, 
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por exemplo, representa a carga do coração humano por um circuito elétrico equivalente (capa 


do livro de Hoppensteadt). 


Circuitos 
elétricos 
análogos 
do sistema 
vascular. 


Figura 2.1.1 A carga do sistema vascular sobre o coração, pode ser representada por uma simples 
resistência elétrica. Se o efeito de bombagem do coração for representado por uma fonte de 
alimentação elétrica, a corrente que circula na resistência representará o fluxo sanguíneo (da capa de 


Hopensteadt). 


Em (Witten (Ed), An electrical equivalent circuit model of glucose-insulin kinetics during 
intravenous glucose tolerance tests in dogs and in man, p. 1188)) encontra-se um circuito 


elétrico que é análogo ao sistema que descreve a cinética do par glucose-insulina. 


Figura 2.1.2. Circuito elétrico não-linear (devido aos díodos) análogo da cinética do par glucose- 


insulina (de Witten,1987) 


Os estudos dos músculos e do movimento humano recorrem intensamente a modelos 


mecânicos, como ilustra a figura seguinte (de McMahon, 1984). 
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Modelo mecânico de um músculo ativo 


Figura 2.2.3. Modelos mecânicos de músculos. Neste capítulo estudaremos o significado dos diversos 


Modelo de uma perna 


componentes mecânicos (de McMahon, 1984) 


Os elementos dos sistemas térmicos, fluídicos e químicos são ferramentas indispensáveis em 
estudos de modelos compartimentais (muito usados em farmacologia e medicina). Por 
exemplo a figura seguinte mostra um sistema fluídico que permite modelizar a ingestão oral 


de um fármaco, a sua absorção pela corrente sanguínea no intestino e a sua excreção no rim 


(de Bruce, pp. 107). 


x(1) a(t) 


y(t) Z(t) 


+ Gu! |-— Blood |— Kicney |— 


O 


a(t) | E : | 
| | vit) 
—— — >» 


———>» 


Z(t) 


Figura 2.2.4. Sistema de três compartimentos modelizado por um sistema fluídico (Bruce, 2001). 
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Qual a vantagem das analogias intersistemas ? Normalmente um Engenheiro Biomédico tem 
facilidade em lidar com circuitos elétricos porque os estudou em Física, e em refletir sobre 
eles, mais do que com sistemas de outro tipo. Reduzindo um sistema qualquer a um circuito 
análogo elétrico, é possível estudar o sistema original como se de um circuito elétrico se 


tratasse. 


Nos parágrafos seguintes estudaremos as noções básicas de sistemas elétricos, mecânicos (de 
translação e de rotação), térmicos, fluídicos, e químicos. Aplicando os princípios básicos de 
cada tipo, obteremos as respetivas equações diferenciais que descrevem as relações entre as 
suas variáveis. Comparando as equações diferenciais de dois tipos de sistemas deduziremos as 


respetivas analogias. 
2.2. Sistemas elétricos 


Um sistema elétrico representa-se normalmente por um conjunto de componentes associados 
às suas propriedades físicas: resistência, capacitância, indutância, fonte de tensão, fonte de 
corrente. Ligam-se esses componentes de acordo com as relações concretas no sistema, e 


obtém-se um circuito elétrico. 
Cada tipo de componente tem a sua missão específica. 


a) Resistência R 


SARA Figura 2.2.1 Resistência elétrica 


Para a resistência R, 

v(t) = Ri(t) (2.2.1) 
A energia dissipada é dada pelo efeito de Joule: 

We =Ri(t” (2.2.2) 


A resistência não armazena energia, dissipa-a apenas. Se se tirar do circuito nada permanece 


nela. E por isso um elemento dissipador de energia. 
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b) Capacidade (condensador) C 


c it 


+ E Para o condensador de capacidade C, 
V 
1 t 
Figura 2.2.3. vt) = = ficodt 
Capacidade elétrica C Cs 


A energia armazenada num condensador é dada por (2.2.4) 


(2.2.3) 


(2.2.4) 


Esta energia não depende da corrente, e permanece mesmo quando esta se extingue, mesmo 


quando se retira o condensador do circuito e se leva para casa. 


A corrente no condensador é dada pela variação temporal da sua carga, ou seja, (2.2.5) 


(09-42 
i(t) = dE 


(2.2.5) 


Mas como Q=Cv, substituindo, aplicando a regra da derivada do produto e atendendo a que Q 


é constante, temos (2.2.6) 


dv 
(O)=C— 
(t) E" 
Ou, resolvendo para a tensão, virá (2.2.7) 
dv 1. 
—=+i(0 . 
dt C O 


Integrando ambos os lados da equação, dá (2.2.8), 
vt) = a Í i(t)dt = 4 fode +Vv 
Ce. Cs É 


sendo O o Instante inicial e vo a tensão inicial. 
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A tensão v no condensador depende da corrente nele acumulada. Se se retirar o condensador 
do circuito, depois de o carregar, a sua energia permanece nele. Por isso se pode chamar 
energia potencial, que se poderá usar colocando o condensador num circuito (tal como a 


energia hídrica contida numa barragem). Uma pilha vulgar é um condensador. 


c) Indutância (bobina) L 


mo 

RE oo oa 

e E Para a bobina L, temos a equação (2.2.9) 

v 
Figura 2.2.4. 
Indutância elétrica 
(t 
v(t) = OO A (2.2.9) 


Integrando agora ambos os lados da equação, dará (2.2.10), 


(t) = - Í v(t)dt = fuga +iy (2.2.10) 


sendo O o instante inicial e io a corrente inicial. 


A energia armazenada numa bobina de indutância L é dada por (2.2.11), 


ea 
W, =—Li (2.2.11) 
2 
e existe, portanto, apenas quando há circulação de corrente. Ela é armazenada no campo 
magnético criado pela corrente quando circula na bobina. Retirando a bobina do circuito, nada 
nela permanece. Por isso a sua energia pode-se chamar de cinética, por depender do 


movimento (da corrente de eletrões). 


Muitos sistemas contêm os três tipos de elementos do circuito RLC: dissipador de energia 
(R), armazenador de energia potencial (C), e armazenador de energia cinética (L). Estes 
elementos recebem energia das fontes (de energia) do sistema. Em muitos casos poderemos 
reduzir os fenómenos no sistema a estas trocas (e dissipações) de energia. Um sistema requer 


normalmente, para ser estabilizável, um dissipador de energia. 
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O e 


Figura 2.2.5. Um circuito RLC. 


Para se obter um modelo matemático do circuito, aplicam-se as leis dos circuitos elétricos: de 


Ohm, de Kirchoff, do condensador e da bobina. 


Lei de Kirchoff das tensões: A soma das tensões é nula em qualquer malha fechada de 


um circuito elétrico. 
Exemplo 2.2.1. 


Considere-se o ponto negro no circuito da Fig. 2.2.6 e uma malha fechada que dele parte e a 


ele regressa. Aplicando a lei de Kirchoff das tensões : 


Vr+tv,+v.—u=0 (lei de Kirchoff) (2.2.12) 


u(b= , PRE NR - calculando agora cada um dos v. 


di = 


RIOEL qe licoar- u()=0 (2.2.13) 


b) 
Figura 2.2.6. Lei de Kirchoff das tensões. a) circuito; b) lei de Kirchoff e sua aplicação. 


Se o nosso interesse for a tensão no condensador C, então essa será a saída y do nosso 
sistema. Aplicando a cada elemento do circuito o que vimos anteriormente obtemos a parte b) 


da fig. 2.2.6. 


dy(t) nº 
(D=C— E: =Cy (2.2.14) 


Derivando ambos os lados da equação (2.2.14), 


: 2 .. 
dO YO cy 
dt dt 


(2.2.15) 


Substituindo na equação da malha e simplificando a escrita obtemos (2.2.16) 
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RCy+LCy+y-u=0. (2.2.16) 


Ou dividindo por RC e simplificando, resulta em (2.2.17) 


ce . 1 1 
+ = u. ZodNT 
PR Ra ( ) 


Obtemos assim uma equação diferencial de 2º ordem que relaciona a entrada u do sistema (o 


circuito) com a sua saída y . Poderemos representar o circuito pelo diagrama de blocos 


Sistema 
(circuito elétrico) 


Figura 2.2.7. Sistema, entrada, saída. 


seguinte 


Entrada 


Os elementos armazenadores de energia deste sistema são o condensador e a bobina. A 


resistência não armazena, dissipa. 


Se quiséssemos guardar informação sobre o passado do circuito (a corrente que nele circulou) 


quais os elementos do circuito que o permitem fazer? 


A carga atual do condensador é o efeito acumulado de todas as correntes que o atravessaram. 
Por isso se pode dizer que o condensador exprime, através da sua carga, o que é o mesmo que 
dizer através da sua tensão, (Q=CV) toda a história passada do ponto de vista da corrente 
(veja-se 2.2.8). O mesmo se pode dizer da corrente na bobina: ela é o efeito acumulado de 
todas as tensões ao longo do tempo (de 2.2.10). Quer dizer que esses dois componentes 
armazenadores de energia, o condensador e a bobina, contêm a história e, portanto, a memória 


do sistema. 


Exemplo 2.2.2. 


Figura 2.2.8. Exemplo 2.2.2 
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Seja o circuito da Fig. 2.2.8. Para o analisar aplique-se a técnica dos nós. 
1) Seleciona-se um nó de referência (em geral um nó de massa). 
11) Escolhem-se como incógnitas as tensões nos outros nós. 
Hi) Escreve-se a equação das correntes de Kirchoff em cada um dos nós. 


Iv) Exprimem-se as correntes de 1ii) em função das incógnitas e dos parâmetros do 


circuito. 
No caso da figura 2.2.8, 
1) O nó da massa é o 4. 
11) As incógnitas serão as tensões v2 no nó 2 e vs no nó 3. 
Hi) Equação no nó 2 : ii = iz+is. 
Equação no nó 3: ixHis=0. 
(Note-se que o número de equações de nó deve ser igual ao número de incógnitas). 


iv) Calculando agora no circuito os valores das correntes em função das tensões nos nós, 


i+i=i o is = = se = (2.2.18) 
1 2 

bti,=062>3 C, no “) q (2.2.19) 
2 


obtêm-se duas equações diferenciais de primeira ordem lineares de coeficientes constantes. 


Escuutudo- ar em fimcão de UEM ede Vi SW-1 
(agora W= Vi-Vê,) 


4 
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Vi — Ve — Ve -(W Ve) => 6 dv, 
R, e de (2.2.20) 

-(VW— 

Ve, ( - E) C2 d (Voe Nes ) =0 

Ra Ad 

. (2.2.21) 
ER À Nas ] 
C, = €, (ata) c, CANA ea +) Ve 
ç 4 
e, 3 -— VW, - DES RR DO (2.2.22) 


2.3. Sistemas mecânicos de translação 


Os sistemas mecânicos estão relacionados com o movimento de corpos. No caso do 
movimento de translação (retilíneo), temos um sistema mecânico de translação. Se o 
movimento é curvilíneo, rotativo, teremos um sistema mecânico de rotação. O corpo humano 
usa os dois tipos nas atividades quotidianas. O projeto de próteses médicas requer 


conhecimentos de sistemas mecânicos e por isso veremos algumas noções elementares. 


Quando um corpo com uma certa massa se desloca sobre uma superfície plana, por exemplo, 
opõe-se-lhe uma força de atrito entre ele e a superfície. Esse atrito absorve parte da força (da 


energia) aplicada ao corpo. Tem por isso uma função de dissipação de energia. 


Certos corpos têm a capacidade de se deformarem quando se lhe aplica uma força: por 
exemplo uma mola helicoidal comprime-se ou expande-se. Muitos corpos são deformáveis, 
alguns facilmente outros muito dificilmente. Para que essa deformação se verifique é 
necessário consumir energia. Mas se libertarmos o corpo ele retoma a forma inicial e devolve 
a energia de deformação. Na adolescência atirámos papelinhos aos nossos colegas de escola 
usando uma fisga de borracha: fornecíamos energia para deformar a borracha, libertávamo-la 
e ela entregava a energia ao papel imprimindo-lhe um movimento retilíneo de grande 


velocidade inicial. A propriedade que determina a deformação é a elasticidade dos corpos. 
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Aplicando forças a corpos teremos que considerar os seguintes elementos essenciais de um 


sistema mecânico: 
- as forças aplicadas, 
- a massa dos corpos, 
- os atritos entre os corpos e as superfícies de deslocação, 
- os coeficientes de elasticidade dos corpos. 


Da aplicação de uma força resulta um movimento, definido em relação a um ponto de fixo de 


referência e de posição nula, representado pelas linhas oblíquas e vertical na Fig. 2.3.1. 


- e 9t4) 
E lt) ; Ee 2lt- 
5 K 
att ou friccalB mols ou 
efestuudeda ke 
(ovjmpad Enfim) 
Aplica -se Ama fria the), produz -"e ANA dadeeamudo gt). 
(em repouts, autis de aplicar LH), 0 dedicando "pulo, o 
fosco 1" a de aajuimua v8o) 


Figura 2.3.1. Os quatro elementos de um sistema mecânico de translação 


Imaginando o movimento de um corpo sobre uma superfície, existe atrito em todos os pontos 
de contacto. Para simplificar a análise, representa-se o atrito em todos esses pontos por uma 
força concentrada num ponto aplicada ao corpo opondo-se ao seu movimento. Representa-se 
por um êmbolo (Figura 2.3.1a)), que funciona de modo análogo a uma seringa: existe uma 
fricção entre as paredes exterior e interior. Essa fricção representa o atrito. Uma extremidade 
do êmbolo não se desloca e por isso liga-se à referência; a outra extremidade liga-se à força 


que lhe é aplicada (normalmente por um corpo em deslocação). 
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O êmbolo, devido à força f(t), sofre um movimento de translação de posição y(t). A 


velocidade desta translação é proporcional à força. 


k Rei te) 
E my PP” e 
8 3 ad” 

atrito cu rca B mol ou 


estuudeda ke 
a 


Aplica-se uma flar the), juodag- arma daditamudo g(t). 
(em pousa, autis de aplicar fr), O dedicando 2º vulto, e 


posicao 1" a de aejuimua vêo) 


- Lua (3 
(2.3.2) 
att B mede K mess M (2.3.3) 


flo = 640 pO=KIO  qt)= dat 


Se analisarmos os sistemas mecânicos de translação do ponto de vista da energia, 
concluiremos que há um dissipador de energia - o atrito B - e dois armazenadores de energia — 


a mola K e a massa M. 


Para comprimir ou expandir a mola é necessária energia. Enquanto está deformada, a mola 
contém uma certa quantidade de energia, que devolve, sob a forma de movimento, quando é 


libertada. Portanto a mola, analogamente ao condensador, tem energia potencial. 


A massa em movimento tem energia cinética. As quantidades de energia de cada um dos 


elementos são indicadas pelas fórmulas seguintes (2.3.4), (2.3.5), (2.3.6). 
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posnidas: Bk,M 
ar maznadars du 4 e 
pala : NE arena) 


(2.3.4) 
mass : La WMM p= 4% ng (emilio) 

alisa padar du enarpãs see su (23.5) 
atrito : wW= B q = B(y Jú (Mam/n ou Was) 

Cprilêwwa) efw-dk (2.3.6) 


O princípio fundamental de análise dos sistemas mecânicos é enunciado pele Lei de Newton: 
a soma das forças aplicadas a um corpo de massa M - a força resultante - produz neste um 


movimento de aceleração. A força resultante é igual ao produto da massa pela aceleração. 


lu A Mais tim (- da mavimade) 
F-=M.a Ba 


Exemplo 2.3.1 Deslocamento de um veículo. 


Considere-se um automóvel de massa M, apoiado num plano com atrito, elasticidade nula (o 
que é de facto praticamente verdade). Aplicando-lhe uma força u ele desloca-se no plano 


segundo a Lei de Newton. 


Qual às fréso aplicado, o vuaner 7 
eu, «pisa do ectercor 
o — afud opsê-sa às ymavi- 
prsiaho, su SAR, aplica MAY 
forca da Brhal fone oo és de U 
A 8y (2.3.8) 


Pode-se representar pela figura seguinte e escrever a lei de Newton (2.3.9). 
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veiha AL 
E. do da Mewlan 


u-By =Hy 


(2.3.9) 


Exemplo 2.3.2 Deslocamento de uma massa com atrito e elasticidade. 


Quous as foras ojlicadm a nº 
o uu, Aplicada de gelirim 


- ad, daRy dC 
as e muda sob 
"diva Es = ky (2.3.1) 


A elasticidade é representada pela mola K ligada à referência. 


Y r—— 


ua dn da Newton (2.3.12) 


U- Isy- 89 = ny 


K sy 


B 
My+By+ky = 44 


Exemplo 2.3.3 Suspensão de um automóvel 


Quando a roda de um automóvel cai num buraco é sujeita a um movimento vertical devido à 
força da gravidade. Esquecendo, para análise, o movimento horizontal, o automóvel pode ser 
representado por uma massa M com um atrito e uma mola a ligá-la à referência (o solo, que 
não se desloca). A mola e o atrito são característicos da suspensão. A arte de construção da 
suspensão está na determinação desses valores. Uma suspensão sem atrito oscilaria para 
sempre, uma suspensão sem mola seria muito “dura” pois não oscilaria. Uma suspensão 


confortável deve ser oscilatória, mas não muito. 


Calculando a equação diferencial do carro-suspensão, 
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ME Mass do axtomeval 
ut fera valia) aplicada as 
vultnt. ( pry acea ol psmiclada) 
M apro ofuntia am burRCa 


My +By+ky = 


as [74 . 4 al MM 
CS BI RI A (2.3.13) 


Obtém-se uma equação diferencial de 2º ordem linear de coeficientes constantes. 
Exemplo 2.3.4 Duas massas interligadas por mola e atrito. 


5 «4 


tu 


4 
Sn elrck a 


Que falas aduama sra as 
eos E dé osflrtada qdênna 
- A fnta da mda Sujaik 


£& yrove 


Aqui é um pouco mais complicado, dado que as massas estão ligadas entre si por atrito e 
elasticidade. Tem que se considerar uma massa de cada vez, e considerar todas as forças que 
lhe são aplicadas. Os deslocamentos das massas são x e y , diferentes entre si. O que tem que 


se considerar para calcular a força do atrito e a força da elasticidade é a diferença (x-y). 


Depois aplica-se a lei de Newton a cada uma delas, como se faz de seguida. 
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E isto Pur 8 c-9) 


ag posa de adrike, toa a ma, 


H, x+8 (x-y d+ K(x-9)=44 


do 


Sue Anta, atuam Sobra H,? 
cado fis, exlrriora 
“E - a fais de moh 


eme! 
bo ota Mevitma e Na | pega Pau K(x-y) 
asi Kbe-y)+8 bi-y)= h29 


(2.3.17) 


a fim da sa E 
a Er ; (65) os opa E fº o (2.3.18) 

2320) | %+ à 9h (x-9) e 

3) [1 a (y ja bn)=o 


Obtêm-se duas equações diferenciais de 2º ordem de coeficientes constantes. Note-se que as 


equações estão acopladas, isto é, x entra na de y”” e y entra na de x”. 


Veremos que uma equação diferencial de segunda ordem pode-se reduzir a duas equações 
diferencias de primeira ordem. E em geral uma equação diferencial de ordem n pode-se 


reduzir a n equações diferenciais de 1º ordem, o que é útil como veremos no Cap. 4. 
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Exemplo 2.3.5 


kzY2 


Be Ya 


“Yya 


Keldo- fe) 6, 8 (9-2) G2-h,9 (2.3.22) 
ki(g Je) tê (3%) 
4 (2.3.23) 
ua, add E 


doe de nwitm a HM, 


ut) mg = My 96) + ByL GE) a 60) + 4 [an/6)-9o(6)] (2324) 


dee do Nintw a Ho 

Nag + 8 (ido) (gd) Boa de = Ma, 03 
(od fa cimo smart Yam dog) de) O 

E dtlici) E (ste) = eg (2329) 


de mes (9 =ya) + B> Ya — ki(9-%a) +k29 = 9 (2,827) 
M, Ma Hz Ma 


Sistema de duas equações diferenciais de 2º ordem, lineares, de coeficientes constantes, 


redutível a um sistema de quatro equações diferenciais de 1º ordem. 
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Comparação com os sistemas elétricos 


Se observarmos com atenção as equações dos elementos dos sistemas elétricos e dos sistemas 
mecânicos de translação verificaremos que têm uma forma semelhante. As variáveis que 


ocupam posições análogas chamam-se variáveis análogas. 
Y : 
tg te 2 
k 
f=8(ffo) f = ln-fe) add | (2.328) 
2 velbudade v qe * acolaraca a; 


aay 


f = B(w-w) f= K(G-I) £=H+ (2.3.29) 


= K (AG -4R) 


er 
lo 
> 
s 
Do 
“Ce 


Vi - 
f= E = ln) Vil = | E .bf fengg= nt 
Ze e (1-0) (2.3.30) 
a) b) c) 


Figura 2.3.10. Analogias entre os elementos dos sistemas elétricos e dos sistemas mecânicos. A 


corrente elétrica está representada por f, para evidenciar as analogias das equações. A velocidade e a 


tensão estão ambas representadas por v. 


alsttrios 


Tabela 2.3.1 Analogias elétrico-mecânico. 
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2.4. Sistemas mecânicos de rotação 


Se o movimento de um corpo é rotativo, ele será descrito por um sistema mecânico de 
rotação. Teremos aqui binário em vez de força, momento de inércia em ver de massa, atrito ou 


fricção rotativa. Os símbolos usados são os da Figura 2.4.1. 
No êmbolo do atrito o pistão interior gira, e entre as duas superfícies existe uma fricção B. 


A elasticidade está relacionada com a torção (e não com a expansão ou compressão). O 


momento de inércia J está relacionado com a inércia que o corpo opõe ao binário. 


- componerdis 


friceas totehve B 
CNemfrasta) 


Aplica - at aum brnavio T(t) » prada -ar same rtaca (desk - 


Comet onpalea) Bem a hos 
rehusuço. Eis cd ig 


Figura 2.4.1. Elementos dos sistemas mecânicos de rotação 


- Luus (2.4.1) (2.4.2) 
friceê [> Corças x aneras J 
Ttt)= & Eai  Te)=ko(t) Th)=T dot 


(2.4.3) 


= B wt(t) = Jx(t) 


Do ponto de vista da energia, teremos uma situação semelhante à dos sistemas mecânicos de 
translação. O atrito é o elemento dissipador, a massa e a mola armazenam energia cinética e 


potencial, respectivamente. 
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- parênito. 8 k,7 
- armezenadars ata ig sa 
Porção : =Whke trend 
nur: QL = Jw2 zw T(0)* (cxinitia), ao 
au duss pad de CAMASAR (Tolo qua Nm ) 4.5 
poco nodticma - AoAçÃ 8 w'- (6) (Mira) 
2=fw.dk (2.4.6) 


(2.4.4) 


O princípio fundamental da análise é a Lei de Newton (para o movimento rotativo): o binário 


resultante é igual ao produto do momento de inércia pela aceleração angular. 
1] 
(o de Muvtim 
= : (2.4.7) 
2 Sd 


Exemplo 2.4.1. Disco suspenso por um eixo fixo na outra extremidade. 


Um disco é suspenso através de um eixo com uma certa torção. Ao disco é aplicado um 


binário, rodando o disco encostado a uma superfície com um certo atrito. 


Vw K2 copo du torção de sixo 
e, 83 Jia , evta o dia 247) 
du La Sputico sxtercor 
TE binawo exteria aflind, 
TS moment: de cnerul em 
relacs Os Xe. 


T 


a: 


Ste) 


Figura 2.4.2. Disco suspenso. 


Quais são os binários plicados ao disco ? 
- binário exterior T(t) 


- binário resultante da torção do eixo 
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Trorção=K 0 (t) (2.4.8) 
- binário devido ao atrito B, 


T. =BO(t) = Bw(t) (2.4.9) 


der da Nilo aplicada as disto 
Tags Ls = 6 eo) 
T-Kro-89-760 ->» To+R6+k9 =T 


ser 2.4.11 
- (2.4.11) 


E em (2.4.11) obtemos uma equação diferencial de 2º ordem, linear invariante. 


Exemplo 2:4.2. Motor com carga acoplada. 


Tm É bnauio do meta 
Tem E momenk da cnaus 


de vmotor + 2x 


ar E Bm Carpe À 
5 7%, By, 8 Cut: da fuicem 
Figura 2.4.3. Motor e carga. AmÁdunS le mé 
O motw ablia Tm ao uxo. Estr k £ Cotf. du torças 
opox-s-lhe com o búnairio da Corças ds MXs 
mars o brnalrco da afrcda. A 2 moment dá snareur 
| da sorga 
O utam afdnica & Ca O seu Om 2 desk comete 


benario da Portão. E sk yu produ ana pulam do mo 
a mtacio da torga(s o tspunia & 2 desle. an. da caga 
da torço KK foste ado, nós ara fnsalnl 
tromanaidar quanto brinanio a! Crpr) 
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a % 


boni d E 
pe dl Figura 2.4.4. Os diversos binários na montagem 
om  F=kB-0,) 


de da Mit oplrada ar ma 
Tm- BK (6-8) = Im Sm (2.411) 
br de Nuntm aflrcada a Cora 
K (0m=6, ) =J 8, (2.4.12) 


“A B . x Tr 
Om El m + z( m L) Tn 


O, — = (6m-8) = O 


(2.413) 


Em (2.4.13) temos duas equações lineares de 2º ordem, de coeficientes constantes, redutíveis 


a 1 equação diferencial de 4º ordem ou a 4 de 1º ordem. 
Exemplo 2.4.3: Orientação de um satélite 


As antenas de um satélite de telecomunicações têm que estar bem orientadas em relação à 
estação terrestre. Se por qualquer razão o satélite gira indevidamente em torno do seu eixo, 


perde-se o contacto. 


Por isso os satélites estão equipados com mecanismos que lhes permitem (sob controlo 
terrestre) orientarem-se segundo as três coordenadas espaciais. Em relação a um eixo 


perpendicular à página, passa-se como ilustrado na Fig. 2.4.5. 


Os propulsores aplicam um binário de rotação, em ambos os sentidos, conforme são 


disparados. 


Admite-se que não há atrito e que o corpo do satélite é rígido (não se deforma), tendo o 


satélite um momento de inércia J. 
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Assim a Lei de Newton aplicada ao satélite é (2.4.14). 


Figura 2.4.5. Orientação de um satélite. 


=iro (2.4.14) 
2.5 Sistemas térmicos 


Nos sistemas térmicos a energia existe sob a forma de calor. Armazena-se em alguns 
elementos, transmite-se entre os elementos. Para os analisar usam-se os princípios de 
conservação da energia: a energia não desaparece, é armazenada no sistema ou transmitida 


entre este e o ambiente. 


Exemplo 2.5.1. Cilindro doméstico 


Num cilindro doméstico, funcionando a electricidade, aquece-se e armazena-se água. Na 
figura 2.5.1 representam-se as transmissões de energia térmica (calor) entre as diversas partes 


do sistema e o ambiente. 
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Figura 2.5.1. O cilindro doméstico e seus 
mecanismos de acumulação e transmissão de 
calor 


O 9, E qudidade du calor frracido pala renina dednia 
O WE quantidade de calor absnida pula pus no exbindor 
(D 208 Quantidade de palor qusar de «lindas (va ajua ) 

O qr8 quantidade ou cabo apa iudir pa cilinda (va ajual 
(E 4% E qualidade de pala au fra (pule pair) pau fra. 


perdendo-se no ambiente 


Aplicando o 
PRINCÍPIO DA CONSERVACAS DA ENTRA 


a e + 
Go +4%: KWt Po Ga | 
color qua sular belas ac sido 
Dé ea 7 


(2.5.1) 


Substituindo agora as expressões de cada uma das parcelas, à esquerda e à direita desta 
equação (2.5.1) obtém-se a equação diferencial da variação da temperatura da água do 


cilindro. 
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a à 
colluds da diverso 24 


-4: frneudo pata resitimio alailania 
em cada a'vatamia, (2.5.2) 
dp = RE 2 uu 


5 %  absnide jrela e el prs no. 
durania «ama apos ou CUnpo E) (2.5.3) 
e M.ec.aT 
ia Ena 


nau VomiaLço ola 
apra espirita decpa (up. de altus 


a o” «Nec AT (54 
7, E Stilo DE = >» %, - de 
dim 6t-»o 


e e TT (2.5.5) 
-% : AGA COMA A apus do Autundro (2.5.6) 


=V.e.T ad milinde o Lomsd V 
to ion dieta de ediado) 


=: pur da -me fuas o onlridi aliases ola farsa 
R$ resina firma oa panda EK, (2.5.7) 


T-Ta fa 
Ja = “e. 
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(2.5.1) 
subida nde na sd bolanto de trerpa 
usVetT: = Me. = Tev eT ad ai: (2.5.8) 
Pol Ur A aa MT e 
T + a M ) = [atá ” M hd Me a (2.5.9) 


Figura 2.5.2. O termômetro de 


mercúrio 
E md 
tmp: Tm 


dp à Uimfunatuna Ta 


Camata calo da erjana Co Vrantiuaa 


guslidada da pala dsmimlia dra nsatudao amas 
cada swAtamÃT : 


q() = Te dic nt e) “JE Ta 


(2.5.10) 


E do termómetro de 


resistência térmica R 
iuusmenk als fados do pustunÃo 
Tm(o) E tim pustua tuial cla merukrio 
Tm (tl — Tm(o 


E suuumenk As Enputuas 
C Ecspanwdedo irma de mucinco (Estal) 


C-Tmlt) £ Qualidade dh delm confia na mmairio 
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de falo mucius dd O ok É 
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(2.5.10) 
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balaio de Lu pa 
ate). at = a[cTmltl] = C.4mbt) CsD 
o STE Na 


: Apa Tdads de 
Atraves de fpamdo cal CrÃE oo ne 
NAMDALIAA ATO 
E Tm o Sm » RC TytIm=Ta 
e ot 4€ -»o 
(2.5.12) 
Tm + e Im = a Ta (2.5.13) 


2 = E (-Te) ind (2.5.14) 


Figura 2.5.4. Componentes análogos eléctricos e térmicos 
e ” > 
Termito det Elietrio 


puiinu. Cluwa R restivor eudna 
capamdoa Cha C elo copiada etica 


fim pinotiaas , potential 
fluie da coln 1? Corruda 


Tabela 2.5.1. Analogias elétrico-térmico 
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Dois sistemas dizem-se análogos se, sendo de natureza energética diferente, têm a mesma 
estrutura (os seus componentes têm funções análogas e sinais análogos) e as suas equações 


diferenciais têm a mesma forma. 


Dois sinais são análogos se, embora de natureza diferente (corrente, calor, ...) têm a mesma 
evolução temporal resultante da resolução das respectivas equações diferenciais: se as 
graficarmos obteremos as curvas com a mesma forma, eventualmente noutra escala. As 


unidades das grandezas são as próprias de cada sistema. 
Dois componentes são análogos se exercem sobre os respectivos sinais o mesmo efeito. 
2.6 Sistemas fluídicos 


Nos sistemas fluídicos circulam fluidos (líquidos e gases). Iremos ter em consideração 
sobretudo o caso dos líquidos. O princípio fundamental de análise é o princípio da 
conservação de massa: líquido não desaparece, ou é armazenado no sistema ou trocado entre 


este e o ambiente. 


Exemplo 2.6.1 Dois tanques comunicantes através de uma válvula. 


fa H,$2,F3 $ como, 
mê/a 


q,= | GEE 06)  h,h 8 mvus da liquido 
k, nf, 


Durante um intervalo de tempo 4t, em cada vaso entrou e saiu uma certa quantidade de 
líquido. No vaso 1 entra gie sai q1; se o primeiro é maior do que o segundo, o nível aumenta. 


No vaso 2 entra q1 e sai q2; se q2 é maior do que gro nível diminui. 


Os canais de comunicação tendem a criar dificuldades à passagem do líquido; por isso têm 
resistência fluídica R1 e R2. O funcionamento de uma válvula de abertura variável baseia-se 


nesse facto: diminuindo a abertura da passagem, aumenta a resistência e diminui o caudal. 
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Para uma certa abertura constante, o caudal depende da diferença das pressões entre os dois 
lados do canal de comunicação. E as pressões são proporcionais ao nível h, se a secção do 


vaso for constante. 


Por isso os níveis variaram 4h; e 4h> e os volumes variam A; 4h: e AzAh> 


far o veso L: A Ok, =q.0t-9 46 (2.6.3) 


Quan da — AO dá 
volume ps se 


for o Voo 2: Az NA = 9,9€ —9,4t (2.6.4) 


Dividindo as equações por 4 e aplicando o limite quando 4t>0, 


A, —. = 9-2, him À, A - 9-9; (2.6.5) 


bt-so ot 
(2.6.6) 


Mas atendendo agora à definição de derivada, 


Vai flt)= Get) - Golt) (2.6.7) 
Rz À, e) = 9ult) = 220) (2.6.8) 


obtém-se um sistema de duas equações diferenciais lineares de 1º ordem de coeficientes 


constantes. 


Estas duas equações podem-se reduzir a uma equação diferencial de 2º ordem. De facto, 
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q, + E / + - Re (2.6.9) 


Derivando agora (note-se que Ri e R> são constantes), obtém-se sucessivamente 


R, q, = h, - h, Ra E < h, (2.6.10) 
k, =R, 4, + 22 3, (2.6.11) 


af Rq + Aa 3, º mt (2.6.12) 


(2.6.13) 


AA 2, Rz 4, + (AR, tha: +A,8) 4, + = (2.6.14) 


até que (2.6.14) é a equação diferencial que relaciona o caudal de saída q2 como caudal de 


entrada qi 


Reparando na expressão 


h-h (2.6.15) 


vemos uma analogia com a lei de Ohm, se a tensão for análoga ao nível, a corrente análoga ao 


caudal, a resistência elétrica análoga à resistência fluídica. 


Num sistema fluídico líquido os elementos armazenadores de energia são aqueles em que o 


líquido se pode acumular, como numa barragem. Um vaso com um certo nível h tem energia 
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potencial: se abrirmos uma saída no fundo sai um jacto de água capaz de produzir trabalho (é 


esse o princípio das barragens). Portanto o nível h exprime energia potencial. 


2.7. Analogias entre os diversos tipos de sistemas 


Entre os diversos tipos de sistemas estudados até aqui existe uma relação de analogia, no 
sentido de que se pode representar, por exemplo, um sistema mecânico por um sistema 
eléctrico cujos componentes são numericamente iguais (nas unidades respectivas) e cujas 


variáveis seguem trajectórias iguais na forma. A tabela seguinte sintetiza essas relações. 


As variáveis circulantes têm valores instantâneos medidos num local. As variáveis 
diferenciais medem-se tomando em consideração dois locais, isto é, requerem um ponto de 


referência. Elas entram ao quadrado na determinação da energia existente no sistema. 


Tabela de analogias 


DISSLPIdO re 


VER. 


Tabela 2.7.1. Analogias entre os cinco tipos de sistemas estudados até aqui. 
O conceito de energia unifica a análise de sistemas. 


Veremos posteriormente que estas analogias são muito úteis, por exemplo nos sistemas 


fisiológicos, facilitando a sua análise. 
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2.8 Sistemas químicos 


A modelização matemática dos processos químicos envolve os princípios básicos de 


termodinâmica, cinética, fenómenos de transporte, etc. 


Um sistema químico recebe materiais e energia do exterior, fornece materiais e energia ao 


exterior e dentro dele acontecem reações químicas. 


Para escrevermos as suas equações aplica-se o princípio geral de conservação de massa e de 


energia (nada se perde, tudo se transforma). 


Imaginemos um sistema químico a funcionar. Durante um certo intervalo de tempo At entram 
produtos, saem produtos, geram-se produtos dentro do sistema, consomem-se produtos dentro 
do sistema, acumulam-se produtos dentro do sistema. Admitindo que At é pequeno, poderemos 
escrever para um componente S do sistema o princípio de conservação, dado pela equação de 


balanço seguinte (2.8.1) 


Durante o intervalo de tempo At, 


dentro do sistema fora do sistema consumida no sistema 


acumulação de S | | fluxo de S para fluxo de S para quantidade de S E quantidade de S A 
dentro do sistema | 


gerada no sistema 


(2.8.1) 
Se dividirmos ambos os lados da equação (2.8.1) por At, obtém-se (2.8.2). 


(2.8.2) 


acumulação de S 
dentro do sistema | | fluxo de S para fluxo de S para a quantidade de S quantidade de S 
“| dentro do sistema | | fora do sistema consumida no sistema 


At gerada no sistema 


Aplicando o limite quando At 50, vem (2.8.3) por definição de derivada. 
(2.8.3) 


acumulação de S 
dentro do sistema | | E de S para | ae de S aa A Ro de S | E de S | 


dt dentro do sistema | | fora do sistema | | gerada no sistema | | consumida no sistema 


S pode ser a massa total, a massa de um componente individual, a energia total. 
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Assim, para um sistema químico genérico, poderemos escrever 


Balanço mássico total: 


Sm SD 0FE-D pr, (2.8.4) 


ientradas j:saídas 


Balanço mássico para o componente A: 


d(n,) e: d(C,V) = x Cali E 5 Ca F, +rV (2.8.5) 


1 
dt dt iientradas j:saídas 


Balanço energético total: 


dE Eno > oEh-D pFh+QtW (28.6) 


dt ientradas j:saídas 
em que 


p & densidade do material do sistema 

p; É densidade do material do fluxo de entrada i 

P; É dendidade domaterial do fluxo de saída j 

V à volume total do sistema 

F, É caudal volumétrico da entrada i 

E: & caudal volumétrico da saída j 

n, É número de moles do componente A no sistema 

c, é concentração molar (moles/volume) de A no sistema 


c, & concentração molar de A na entrada i 


A; 

c, & concentração molar de A na saída j 

r É taxa de reacção por unidade de volume para o componente A no sistema 

h, É entalpia esecífica do material na entrada i 

h, à entalpia específica do material na saída j 

U,K,P & energia interna(U), cinética(K) e potencial (P) 

Q £ quantidade de energia trocada entre o sistema e o seu ambiente por unidade de tempo 
q g p p 


W, & trabalho (shaft work) trocado entre o sistema e o seu ambiente por unidade de tempo 


Aplicando a um sistema qualquer obtêm-se equações diferenciais de 1º ordem, que constituem 


o seu modelo matemático. 
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Exemplo 2.8.1. Aquecimento de um tanque bem agitado (stirred-tank) (de Stephanopoulus) 


Seja um tanque que está a ser alimentado por um caudal de líquido Fi (m?/min) à temperatura 
Ti (ºC). No tanque é aquecida por uma serpentina a vapor onde entra um caudal de vapor de Fy 
kg/min. Do tanque sai um caudal F de líquido à temperatura T. O tanque está bem agitado e por 
isso a temperatura da saída é iguala temperatura dentro do tanque. O vapor de aquecimento, ao 


perder energia na serpentina condensa, saindo o condensado para o exterior. 


Figura 2.8.1 Um tanque bem agitado (de Stephanopoulus) 


Aplicando os princípios genéricos indicados acima 
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Massa total no tanque Mt: 
Mt=pV=pAh 


sendo p a densidade do líquido, A a secção do tanque, h o nível do líquido. 


Energia total do líquido no tanque Et: 
Et=U+K+P 


Não há variação da energia potencial nem da cinética porque o tanque não se move. 


Por isso a variação da energia será apenas a variação da energia interna U. 


Para líquidos, sendo H a sua entalpia, 


(2.8.7) 


(2.8.8) 


U H 
du dH (2.8.9) 
dt dt 
H=pVe(T-T.) (2.8.10) 
sendo 
E: à capacidade térmica do líquido no tanque 
Ter à temparatura de referência à qual se assume que a entalpia do líquido é nula. 
Aplique-se agora o princípio da conservação da massa e da energia. 
Conservação da massa: 
aa da | 
a total massa total | | massa total | 
o Ages e ORA ali Da (2.8.11) 
dt entrada saída 
ou seja 
d(poAh) PAdh dh 
E >0F-oF > >—=p0F-oF 5DA--=E-F > 
dE PE-P ER PE-P AR 
dh FE F (2.8.12) 
dt A A 
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Conservação da energia 


ã E de 


energia total | * | total de energia total de energia ha energia fornecida (2.8.13) 
dt entrada saída pelo vapor di 
ou seja 
d[pAhe, (T=T.)] 
= PR (TT) =PFe(T To) +Q 
PAC d[h(T =T,)] 
> > dt ' = pFc,(T-Ley)-pFc(T-T)+Q 
Ad[h(T —T 
= SSD e pl La) F(T Tipos 
dt pc, 
= AdlhT] pr rT+ 2, assumindo T, =0 (2.8.14) 
dt pc, 
Ad[hT]  anYT AT dh  anIT (E = E) 
dt dt dt dt 
= Ah T(E-P)=FT, pr 
dt PC, 
ss amil ra T+ LS 
dt pc, 
RELA E F(T —-T)4 o, (2.8.15) 
dt Ah Ah pc, 
Juntando as duas equações 
fita no (2.8.16) 
dt A A 
Ss EM -T)+ Too (2.8.17) 
dt Ah Ah pc 


obtém-se um sistema de duas equações diferenciais de 1º ordem não lineares acopladas. 


Note-se que este sistema químico é a conjugação de um fluídico com um térmico. 
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2.9. Sistemas biológicos 


Os sistemas biológicos descrevem o funcionamento dos seres vivos. Vejamos o exemplo da 


modelização da evolução da população de espécies. 
2.9.1. Crescimento de uma população de uma espécie biológica 


Suponhamos uma espécie biológica cujos indivíduos nascem e morrem na mesma estação 


(isto é, vivem um número inteiro de anos). 


Inicialmente existem No indivíduos. Ao fim de um ano, N; indivíduos, etc., ao fim de k anos 


Nk indivíduos. 
Em princípio poderemos admitir que 
Ni= ANo (2.9.1) 


em que A é uma constante positiva que depende das condições ambientais (existência de 
alimentos, água, clima, etc.). Se A> 1, a população aumenta; se A<1, diminui até à extinção. Se 
A>1 é constante ao longo dos anos, a população crescerá ilimitadamente, levando à explosão 


populacional Malthusiana. 


Ora se a população crescer sem limite, por um lado não haverá alimentos suficientes, e por 
outro lado os predadores terão mais facilidade em caçar indivíduos da espécie em análise, e por 
isso o crescimento da população será limitado. Este efeito de limitação pode ser introduzido na 
equação através de um termo subtrativo que seja insignificante para pequenos valores de N 


(pequenas populações) e muito influente para grandes valores de N, por exemplo 
Ni= ANo — BNo? (2.9.2) 
Se B<<A, o termo subtrativo só será relevante para grandes valores de N. 
Nos anos seguintes teremos 
No= ANi — BN? (2.9.3) 
N3= AN> — BN? (2.9.4) 


Haverá uma altura m em que a população não poderá crescer mais, isto é, em que atinge o seu 


valor máximo, Nmax. No instante seguinte 
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Nm+1= ANmax — BNma? (2.9.5) 
para que Nm+1 seja positivo, 


AN BN. 20>N, (A-BN, )>0=>A-BN,(>0 


max max 


A 
SA BN => Naa É = (2.9.6) 


Encontrámos assim o limite superior da população, A/B. 


Exprimindo a população como uma fração da população máxima, por xk, 


N 
Xe N... (2.9.7) 
teremos, da equação anterior, substituindo, 
Nua - AN, —- BN; > Kat amas = AXN ax —-BxN a > Xys1 = Ax, BN a 
A 
> Xe AX, “Br 2 Mei = Ax(1-x,) (2.9.8) 


Temos assim encontrado o modelo de crescimento populacional, na forma de uma equação de 


diferenças não linear, 
Xu = Axl x)= f(x) xe[0,1] (2.9.9) 
em que xk é a população no ano k, como uma fração da população máxima. 


Haverá estados de equilíbrio, em que a população se mantém constante ao longo do tempo? 


Nestas circunstâncias teriamos 


Xys1 = 
ou seja, 
Xks1 =Ax(l-x,)=x, = Ax, — Ax, ni Ad o) aii (2.9.10) 
pe en Ko - 
> — Rates s 
x,=0 
X, =— Xk =0 


Se A>1, temos assim dois estados de equilíbrio. Se A<1, só a origem pode ser estado de 


equilíbrio, dado que xk tem que ser maior do que zero. Note-se que A>0. 
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2.9.2. Modelo Lotka-Volterra de predador-presa 


Este modelo descreve a interação biológica entre duas espécies, na relação particular em que 


uma (https://services.math.duke.edu/education/ccp/materials/engin/predprey/pred2 html, 22/07/2022 ) 


, à predadora, come a outra, a presa. Há muitos exemplos na natureza: leão-gazela, pássaro- 


inseto, pandas-eucaliptos, etc.. Para não complicar o modelo façam-se as seguintes hipóteses: 
- a espécie predadora é totalmente dependente da presa (é o seu único alimento), 

- a presa tem uma quantidade ilimitada de alimentos, 

- a presa só tem um predador, o que está a ser considerado e não tem outras limitações. 


Na natureza as coisas são mais complicadas, mas assim poderemos desenvolver um modelo 


aceitavelmente complexo. 


Se existirem zero predadores, a presa não tem qualquer limite ao seu crescimento e por isso 
cresce exponencialmente. Sendo x(t) a população atual no instante t, a sua taxa de crescimento 


será 


ato 


É ax(t) (2.9.11) 


Mas havendo predadores, o crescimento será menor: os predadores introduzem um fator 


negativo nesta equação. 


Se y(t) for a população de predadores no instante t, cada um deles tem certa probabilidade de 
encontrar uma presa. Essa probabilidade é tanto maior quanto mais predadores existirem e tanto 


maior quanto mais presas existirem. Por isso se pressupõe que 


- a taxa de encontros entre um predador e uma presa é conjuntamente proporcional ao tamanho 


das duas populações. 
- dos encontros realizados, uma parte fixa resulta na morte da presa. 


Assim teremos a equação (2.9.12) para as presas 
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a = ax(b)- bx(Dy(b) . (2.9.12) 


Quanto à população de predadores, y, se não tiver alimentos morre a uma velocidade 


proporcional ao seu tamanho, ou seja, neste caso, 


DO 
O (2.9.13) 


por não haver energia para novos nascimentos. 


Ora felizmente para os predadores existem presas, e por isso eles vão aí buscar energia para se 
reproduzirem. Por cada encontro predador - presa há uma parte que resulta em alimentos para 


o predador e por isso a sua população variará de acordo com (2.9.14) 


eo = =cy(D + pX(O (O (2.9.14) 


sendo p uma constante. 


Juntando as equações obtidas obtém-se o modelo de Lotka-Volterra 


o = ax(t) — bx(t) y(t) (2.9.15) 
DO ey + pro (2.9.16) 


composto por duas equações diferenciais de 1º ordem acopladas, que não podem ser resolvidas 
separadamente e para as quais não é possível encontrar uma solução analítica (mas apenas 


computacional). 


Pode-se mostrar que existe um par (xs, Ys) para o qual dx/dt=dy/dt=0, isto é, as populações estão 


em equilíbrio e não variam ao longo do tempo. Para isso basta fazer 
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0 = ax(t) — bx(t) y(t) 
O = —cy(t) + px(t) y(t) 


Y(t) = 


” [a =by()()=0nx(050 I [a -by()]=0 (2.9.17) 


[-e+ pro) =0n 40040" |t-c+ prol=0 "1 o- 


lo cla 


um tal estado chama-se regime permanente ou estado estacionário. O seu valor depende das 


constantes a ,b,c ,e p, os parâmetros do modelo. 


Voltaremos a este modelo mais tarde (trajetórias de fase, estados de equilíbrio, ...). 


Existem modelos mais completos de duas populações. Por exemplo em 


https://www.montana.edu/screel/teaching/bioe-370/documents/L-V competition populus.pdf 


encontra-se um modelo que inclui competição (luta pelos mesmos recursos) entre espécies, 
interação agressiva (combatem-se ativamente, embora não se comam), cooperação ou simbiose 
(ajudam-se mutuamente), relação predador-presa, interação forte-fraco (uma espécie está mais 


capacitada para sobreviver). 


2.10 Sistemas fisiológicos 


Os sistemas fisiológicos incluem os organismos dos animais em geral e do homem em 
particular. São dos sistemas mais complexos que a natureza criou. Eles compõôem-se em grande 
medida de sistemas dos tipos que vimos (mecânicos, fluídicos, químicos, etc.) que ao incluírem 


órgãos vivos adquirem, por sinergia, novas propriedades. 


Exemplo 2.10.1. Um vaso sanguíneo 


Considere-se um vaso sanguíneo ilustrado na figura. 
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V ê volume do vaso 

Q, 2 caudal de entrada 

PB 2 Pressão à entrada do vaso 
Q, é caudal de saída 


P, É Pressão à saída 


P, (Pressão exterior) & O 


Figura 2.10.1 . Representação de um vaso sanguíneo (de Hoppensteadt , p. 9) 


Para encontrarmos uma relação entre Q, P1, P> e V consideram-se duas propriedades do vaso: a 
sua resistência ao fluxo e a sua complacência (compliance) para com a pressão (que tende a 


aumentar a sua secção). A fim de facilitar a análise idealizem-se duas situações. 
(1) Um vaso sem complacência, rígido, de volume constante e bem conhecido: 


Admita-se que o vaso está a funcionar em regime permanente em que todas aquelas quantidades 
são constantes. Nessas condições, necessariamente Q:=Q», porque não há acumulação de 
sangue no vaso por ser constante o seu volume. Seja Q= Q1=Q» esse valor do fluxo. Para ele 


teremos, como nos sistemas fluídicos, o caudal Q dado na Fig. 2.10.2. 


Sp q=4E 


Ressdimna R 


Figura 2.10.2. Um vaso rígido, puramente resistivo 


A R chama-se resistência do vaso e teremos aqui um vaso de resistência. 


(11) Um vaso elástico sem qualquer resistência fluídica: 


Nestas condições a pressão na entrada é igual à pressão na saída, para qualquer valor de Q. A 


relação entre a pressão P e o volume V do vaso pode ser aproximada pela relação da Fig 2.10.3, 
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/ ” e V 

/ ! V=CP 5 C=— 

aopot dya P 
Combla unua Cc Q, * Q, 


Figura 2.10.3. Um vaso puramente elástico, sem resistência 


em que a constante C se chama a complacência do vaso. Estamos a admitir que o vaso tem um 
volume Va nulo quando a pressão é nula (volume “morto”, dead). Se assim não for teremos 
V=V,+CP. 

Estas situações 1) e 11) idealizadas exprimem a realidade com alguma aproximação. Um vaso 
real tem resistência e complacência ao mesmo tempo. A relação entre as variáveis não é 
exatamente linear, mas a linearidade é uma hipótese simplificadora. 

Consideremos agora o vaso simultaneamente com resistência e complacência, num regime não 
estacionário em que o caudal de entrada Q1 é diferente do caudal de saída Q2, como na Fig. 
2.10.1. Nestas condições ou há enchimento (o vaso aumenta de volume, se Q: >Q2) ou 
esvaziamento (o vaso diminui de volume, Q1<Q2). Durante um pequeno intervalo de tempo, At, 


a variação de volume será (2.10.1), pela lei de conservação da massa, 


AV=QiAt — QrAt (2.10.1) 
ou seja (2.10.2) 

AV 

rd, a 2 (2.10.2) 


Ou ainda, aplicando limites, (2.10.3) 


lim S=T=V0=40-Q(. (2.10. 3) 


At>0 At 


Pela definição de complacência C, como relação entre o volume e a pressão, tem-se (2.10.4). 
V(t) = CP(t) (2.10.4) 


A pressão de referência Po é a do exterior do vaso, e P(t), a pressão a que a complacência está 
sujeita, é medida em relação a Po . Mas levanta-se agora a questão: em que ponto devemos 
considerar a pressão P(t) ? Podemos considerá-la à entrada, como Pi(t), ou à saída, como Pa(t) 


(note-se que estas duas pressões são diferentes), ou num ponto intermédio qualquer, por 
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exemplo a média das duas ? Esta questão será discutida mais à frente. Por agora aceitemos a 
escolha de Pa(t). 
Assim teremos (2.10.5) 

VC) =CP(t) . (2.10.5) 


E agora derivando ambos os lados da igualdade , 


VC) = CÊ(t) (2.10.6) 
ou seja (2.10.7) 
Bp='e.. (2.10.7) 


Substituindo agora (2.10.3) em (2.10.7) obtém-se (2.10.8) 


. Q Q 
Pt) = Rae — a (2.10.8) 


O vaso tem uma resistência fluídica R. A relação entre a pressão e o caudal é, à entrada 


(2.10.9) 


Pi(t)-Po 


GD)-— (2.10.9) 
e à saída (2.10.10) 
Qu(t) = 2d (2.10.10) 
Substituindo (2.10.9) e (2.10.10) em (2.10.7) obtém-se finalmente (2.10.11) 
[A (a pe cido PER RS to (2.10.11) 


RC RC RC RC 


Vejamos agora a analogia entre o enchimento do vaso e a carga de um condensador. A 
complacência controla o aumento do volume do vaso e a capacidade controla o aumento da 
carga elétrica do condensador. Conforme a complacência C, o vaso alarga-se até chegar a um 
volume fixo, V, que depende da complacência, sendo V=CP, em que P é a diferença de pressão 
entre os dois lados da parede do vaso. Conforme a capacidade C, o condensador carrega-se até 


chegar a uma carga fixa, Q, que depende da capacidade C. Aí Q=CV. 
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Por isso um vaso sem resistência e com complacência pode ser analogamente representado por 


um condensar, como na figura 


V É 


Compladnwa € 
Ve CP 


ú C. 


Q =CV 
(coma) 


Figura 2.10.4. Analogia entre a complacência e a capacidade 


A resistência do vaso à circulação do sangue é a resistência fluídica que vimos nos sistemas 
fluídicos. Se o vaso for rígido (Fig. 2.10.5) não tem complacência e o caudal de entrada é igual 
ao caudal de saída. Há uma diferença de pressão entre a entrada e a saída e o caudal depende 


dessa diferença. Como sabemos a resistência fluídica é análoga da resistência elétrica, tal como 
na Fig. 2.10.5. 


' E Lo 
Sp 8 vo 
Ressdimna R 
pf o aÇh 
> ———— L>——+ 
8 R R 


Figura 2.10.5. Analogia entre a resistência fluídica do vaso rígido e a resistência elétrica. 


Juntando as duas situações obteremos a seguinte tabela de analogias, que se extrai da 


comparação das equações dos dois sistemas. 


Vaso Circuito elétrico 
Pressão, P V, Tensão 
Volume, V Q, Carga, 

Caudal, Q L, Corrente, 
Resistência, R R, Resistência 
Complacência, C C, Condensador C 


Tabela 2.10.1. Analogias fisiológico-elétrico. 
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Vejamos agora como poderemos obter um circuito elétrico análogo a um vaso que tenha 
simultaneamente complacência e resistência. Em primeiro lugar deveremos constatar que estas 
propriedades se distribuem uniformemente ao longo do vaso, desde a entrada até à saída. 
Poderemos dividir o vaso em secções pequenas, tal como na Fig.2.10.6, e para cada uma delas 
considerar a sua complacência e a sua resistência. Teremos para cada uma delas uma resistência 


e uma capacidade como análogos elétricos. 


Figura 2.10.6. Analogia entre um vaso e um circuito elétrico. Cada secção é análoga a um 


circuito R-C. 


Trata-se de um sistema de parâmetros distribuídos, e não de parâmetros concentrados. De facto, 
na natureza, as propriedades dos elementos distribuem-se normalmente pela extensão dos 
elementos. Em engenharia, por necessidade prática, simplifica-se a análise para a tornar tratável 
e supõe-se que elas estão concentradas num ponto. Só assim poderemos desenhar, por exemplo, 
um circuito com R-L-C no caso dos elétricos. Trabalhar com parâmetros distribuídos torna-se 
praticamente impossível. Fazendo essa aproximação ao circuito aos sistemas da Figura 2.10.6, 


obteremos finalmente a Fig. 2.10.7 para um circuito elétrico R-C análogo a um vaso sanguíneo. 


P. =0 Ro pan cd v2 


Figura 2.10.7. Analogia entre um vaso e um circuito elétrico no contexto de parâmetros 


concentrados. 


Esta analogia torna-se evidente se escrevermos as equações diferenciais para a variação da 


pressão P> na saída do vaso e da tensão vz aos terminais do condensador. 


Pivariável: vi variável : 
dp PB n E (2:10:11) dv, MY dv hM 
dt RC RC dt R dt RC RC (2.10.12) 
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Para o vaso temos a equação (2.10.11), que obtivemos anteriormente, e para o análogo elétrico 
a (2.10.12). Fica assim justificada a escolha da pressão à saída do vaso no desenvolvimento de 


(2.10.11). 


No corpo humano a circulação produz uma divisão entre os vasos de resistência e os vasos de 
complacência. As grandes artérias e veias são sobretudo de complacência, porque bastam 
pequenas diferenças de pressão para produzir o fluxo de sangue, e o seu volume sofre variações 
muito significativas. As pequenas artérias e veias nos tecidos (arteríolas e vénulas) são 
sobretudo vasos de resistência: o seu volume varia pouco e verificam-se aí grandes diferenças 


de pressão. 


Exemplo 2.10.2. Músculo suspenso por uma extremidade 


Imagine-se um músculo de rã fixo por uma extremidade e no qual se suspende uma corpo de 
massa M (de Bruce). O músculo, devido à sua elasticidade, alonga-se. Pode-se desenhar um 


equivalente mecânico de translação com uma mola K (elasticidade), atrito viscoso R, massa M. 
A força aplicada à massa é a da gravidade e igual ao seu peso: F=Mg. 


O atrito e a mola reagem com as forças (b) 


; dy - 
atrito: F=B--=B 2.10.13 K 
B dt y ( ) R 
mola: F.=Ky (2.10.14) 
weight 
(mass = M) 
A força resultante é a soma das três e produz (Fig. 2.10 Bruce) 
Mg 
piissincamCa tao Figura 2.10.8 (de 
. o. o. . Bruce) 
Mg-Ky-By=My > My+By+Ky=Mg (2.10.15) 
e E vei (2.10.16) 


M íM 
Obtendo-se assim uma equação diferencial de segunda ordem linear 
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Exemplo 2.10.3 Modelo biomecânico da inspiração/expiração pulmonar 


O ato de inspirar e expirar pode ser descrito por um sistema biomecânico ilustrado pela figura 


seguinte. 


Figura 2.10.9. A biomecânica da inspiração (Bruce, p. 31). 


No estado de relaxamento, sem qualquer fluxo de ar, a pressão dentro dos pulmões é igual à 
pressão atmosférica, PB. Durante a inspiração os músculos respiratórios — isto é, o diafragma e 
os músculos da caixa torácica contraem-se e baixam a pressão (Ppi(t), pressão pleural) em 


torno dos pulmões dentro do peito. 


Devido à sua elasticidade (descrita pela complacência C, o inverso da elasticidade), os 
pulmões expandem-se devido ao abaixamento da pressão transmural dentro deles e por isso a 
pressão interior dos pulmões (representada por Pr(t)) decresce momentaneamente. Gera-se 
assim uma queda de pressão na resistência fluídica Raw das vias respiratórias, donde resulta 
um caudal de ar para dentro dos pulmões. Enquanto os músculos se contraírem, com força 
crescente, este processo repete-se e a inspiração continua. No fim da inspiração os músculos 
expiratórios do abdómen e da caixa torácica. contraem-se, forçando Ppi(t) a ultrapassar Pp 
(pressão atmosférica) e expulsando o ar dos pulmões. Na respiração de descanso os músculos 
de expiração não se contraem, basta a relaxação dos músculos de inspiração para que Ppi(t) 


ultrapasse Ps. 


Este processo biomecânico pode ser descrito matematicamente a partir das propriedades 


físicas das estruturas envolvidas. Trata-se de um sistema misto fluídico-mecânico. 


O fluxo de ar depende da diferença de pressão através da resistência fluídica Raw e da própria 


resistência. 


(DADC/FCTUC/2022 Licenciatura em Engenharia Biomédica/MCPF 66 


Dinâmica dos Sistemas Biológicos e Fisiológicos Cap. 2 Modelização por equações diferenciais 


Se I(t) representa o caudal de ar de inspiração/expiração, teremos, 


P,- Pt) 


HO = E 


ouseja P.-P(D=I(D)R,y (2.10.17) 


Sendo V(t) o volume dos pulmões, durante o intervalo de tempo At a variação do volume será, 


admitindo I(t) constante durante At, 


Vt+ o -VO=I(LA (o ADO =I(t) 
= I(0=V.(0) (2.10.18) 
e por isso, substituindo em (2.10.17), 
EPI VIORG o (2.10.19) 


A complacência dos pulmões Cr é definida pela relação entre o volume dos pulmões e a 


pressão transmural 


Ra (OE (2.10.20) 
P — Pa(t) 
donde 
p-B()-4O > pop (+40 (2.10.21) 
C, C, 


Substituindo agora Pr na expressão acima, 


P, SPL) =V(O Ray > P, por vio Rs, 
E (2.10.22) 


> B-P(D=V(ORy + = 


L 
Como a pressão atmosférica é constante, define-se a pressão diferencial 


P()=P,-P,(t) (2.10.23) 
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e 
PO)=VORy +40 > 
C, 
e 1 1 
Vi (D+ RG, O “Ra, P(t) (2.10.24) 


ou 


V(O+ av, (t) =bP(t) 


obtendo-se uma equação diferencial de 1º ordem, linear, que relaciona o volume dos pulmões 
VL com a pressão relativa intramural P(t). Do ponto de vista sistémico pode representar-se 


pelo diagrama de blocos 


P(y) Vitt) 


Sistema de respiração 


Figura 2.10.10 
Note-se que 


VO=IO) S V(O=V(O) + [rode 
o (2.10.25) 


em que Vi(0) é o volume em repouso. Como é constante (característica do indivíduo), pode 
considerar-se nula e VL(t) representa por isso a variação do volume dos pulmões em relação a 


Vi(0). 

2.11 Outros tipos de sistemas 

Podem-se desenvolver modelos para sistemas financeiros, de gestão, sociais, etc. 
Exemplo 2.10.1. Empréstimo bancário 


Considerando x 
= valor do empréstimo total 


o) 


A ag A 
U, = prestação entregue no mês k 


k 


r É taxa de juro, 100r%, mensal 


A a A 
Ve = montante em dívida no mês k 
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o que se fica a dever para o mês seguinte é igual ao que se deve no mês corrente mais os juros 
sobre os que se deve no mês corrente, menos a prestação uk que se faz no mês corrente. Ou 


seja, é uma equação de balanço financeiro dada por (2.11.1). 


Vu AO Upa UA (2.11.1) 


2.12. Simulação das equações diferenciais e de diferenças. 


As equações diferenciais podem-se simular em circuitos eletrónicos com componentes 


apropriados. 


No caso de equações lineares de coeficientes constantes bastam blocos que executem as 


operações de adição, multiplicação por uma constante e integração. 


Três bloco dutiimi botes : 


Y, V; Ve 
o E 
é LS 


1; Vo = AY; vo = % digo, 
Semader Am bl: hador Integrados Ee, 


Figura 2.12.2 Os três elementos de simulação de sistemas dinâmicos lineares invariantes. 


Com eles poderemos simular qualquer equação diferencial linear invariante. 


Vejamos alguns exemplos. 
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Exemplo 1. Uma equação diferencial de 1º ordem 


Fig.2.12.3 


eg - dt. da (“oveana — tw neussaArs sum ariana do 


Estes modelos de simulação são implementáveis fisicamente , com o chamado computador 


analógico . 


Permitem analisar o comportamento de sistemas físicos através da observação das variáveis 


do circuito analógico análogo às do sistema físico. 
Permitem ainda encontrar outra representação do sistema (equações de estado a ver no Cp. 4). 


Exemplo 2. Simulação do posicionamento de um satélite 


Em 2.4.14 vimos o modelo de um satélite como um sistema de rotação, que aqui se repete. 


à ” Lg 
Jo=T s 9=5T a O A 


8 +0.0 + 0.9 = zT 


Fig.2.12.4 
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Exemplo 3. Sistema de duas equações diferenciais de 1º ordem 


4+Skh=u G=-Sfa+4 
Y-3y =o Ja = 39, 


Fig2.12.5 


Exemplo 4. Equação diferencial de 2º ordem 


J+93sy-29=4 STS 


Y=-95y+2y+4 y=/ydt= [q det 


Fig.2.12.6 
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Exemplo 5. Equação diferencial de 3º ordem 


Trageby+ eg 4 


id / / Fig.2.12.7 


Exemplo 6. Sistema de duas equações diferenciais de 2º ordem: 


], E? 3h + NH =M; 


e. . 4» 
Ya +91%-4 -Sy +iSh =24u, 


f=-28 429 +M, 
No =—919,+9, +54, - 4S À +? Uz 


podendo-se desenhar o diagrama de blocos: 
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Fig.2.12.8 


2.13 Sistemas contínuos com atraso puro ou atraso de transporte 


Em muitas situações transporta-se energia ou massa entre dois pontos sem qualger 


transformação de permeio 


| - 
: “CG Ywtt-U) 


Tá: temp. raassero psre o transhovk dada a vluula 
al a! calda de Can ça. 
Fig.2.13.1. Atraso puro ou atraso de transporte 


Uma tremonha de cereais apresenta um atraso de transporte entre a entrada e a saída dependente 


do comprimento e da velocidade do tapete rolante. 


Se aquecermos uma barra metálica numa das extremidades, o calor demora algum tempo a 


chegar à extremidade oposta. Admitindo que a barra está termicamente isolada (não se perde 
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energia térmica), a temperatura em T> seguirá a mesma evolução da de Ti, mas com um atraso 


Mt) | Ee 
6 € 
Ta (tb) 
Fig 2.132. E açãe 
Td — 


Se y=, Ê Telbla Ti(t-ta ) 
U(8) 42, Gt) rasylt)= Cute) 


puro, 
cselanha labmere 


T, Bee metilu. | E 


Se V=le, | 
dr Gltlea,s te)ra gt) = bw (E-R) 


ue À 
Sferamuto 17 Ai feremiat com elis bar. 


O atraso puro (diz-se assim porque o sistema só aplica atraso à entrada, abstendo-se de qualquer 
outra ação) existe em muitos processos e normalmente é um problema. Aparece geralmente 


combinado com outros efeitos. 
2.14 Sistemas não lineares 


Para sistemas não lineares, em que por exemplo as variáveis se multiplicam umas pelas outras 
(como no modelo de Lotka-Volterra), são necessários blocos adicionais que executem essas 


operações. 


O SIMULINK (Matlab) simula blocos eletrónicos (pode-se dizer que implementa blocos 
eletrónicos por programação). Contém um grande conjunto de blocos que permite implementar 
facilmente qualquer operação matemática e portanto simular qualquer sistema (linear ou não 
linear). Por exemplo para o modelo de Lotka-Volterra tem-se o seguinte diagrama de blocos no 
SIMULINK, em que as constantes a, b, c, p, são implementadas por amplificadores (quadrados 


na figura). No capítulo 3 veremos porque 1/s é um integrador. 
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Diagrama no SIMULINK 


Modelo de Lotka-Volterra 


Osciloscópio 


a = ax(t) — bx(t) y(t) 
DO ey()+ pro Fig 2.141 


A simulação de sistemas discretos, definidos por equações de diferenças, também de pode 


fazer no SIMULINK. Agora o bloco integrador é substituído pelo bloco atraso (unit delay). 


De facto, para a equação de diferenças por exemplo do empréstimo bancário, linear 


Ye = Ye try — Um = U+r)y,—u, 


Poderemos definir o operador atraso z!, e o de avanço z tal que 


2 Deal=) Ay )= Yi z[z [yall= 2y]= Ye 


Poderemos representar aquela equação do empréstimo pelo diagrama de blocos seguinte 


= + 
Prestação od 
Visualização 


Figura 2.14.2 Diagrama SIMULINK do empréstimo bancário 


Para o modelo de crescimento populacional, não linear, teremos que usar o bloco de produto. 
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Modelo do crescimento de uma Diagrama no SIMULINK 


população isolada 


Visualizaçãe 


Xu = Ax(I-x,) = Ax, — Axi 


Fig.2.14.3 


Qualquer equação diferencial se pode reduzir (aproximar) a uma equação de diferenças, 
aproximando a derivada (pelas diferenças para a frente, por exemplo) ou o integral (pela regra 


trapezoidal, por exemplo). 
2.15. Conclusão 


Na generalidade dos sistemas podem-se aplicar os princípios básicos (do respetivo domínio) 
para se obter um conjunto de equações diferenciais (nos sistemas contínuos) ou de diferenças 
(nos sistemas discretos), que são uma representação matemática do sistema e constituem por 


isso um modelo matemático. 


Esse modelo servirá, como veremos, para o estudo do comportamento do sistema (a sua 
evolução temporal, por exemplo) e contém características matemáticas que exprimem 


propriedades internas essenciais do sistema. 


Obtido o modelo, importa saber como se vai usar, e é esse o tema dos próximos capítulos. 
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3.1. Introdução 


A obtenção das equações diferenciais estudadas nos capítulos anteriores é uma etapa 
importante no desenvolvimento de modelos computacionais. A partir delas desenvolvem-se 
modelos matemáticos em dois espaços conceptuais: (1) função de transferência e (11) espaço de 


estados. 


A função de transferência representa o sistema em análise de uma forma compacta e intuitiva. 
Há parâmetros dessa representação que estão estreitamente relacionados, de modo conhecido, 
com propriedades estruturais do sistema, e por isso conhecida a função de transferência, é 
possível deduzir facilmente essas propriedades. O seu desenvolvimento decorre no domínio 
complexo, pois trata-se de uma função de variável complexa. Representa facilmente o sistema 
no domínio da frequência (substituindo s por jw). Tem, no entanto, duas fortes limitações: não 
é prática para a programação computacional e aplica-se apenas a equações diferenciais 


lineares e de coeficientes constantes. 


TRANSfONNADA DE CAPLACE & 


Figura 4.1.1 Ilustração da utilidade da transformada de Laplace 


No espaço de estados definem-se as equações de estado, que são nem mais nem menos do que 
um arranjo em forma padronizada das equações diferenciais. Estudá-lo-emos no Capítulo 4. E 


de aplicação mais geral, para o caso linear e não linear, variante e invariante. 
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3.2. Aplicação da transformada de Laplace às equações diferenciais 


A transformada de Laplace é uma boa ferramenta para a resolução de equações diferenciais. 
Por exemplo no Capítulo 2 vimos o sistema mecânico da suspensão de um automóvel, para o 


qual obtivemos a equação diferencial de 2º ordem 


Imagine-se um músculo de rã fixo por uma extremidade e no qual se suspende uma corpo de 
massa M (de Bruce), como no Cap. 2, pág. 65. O músculo, devido à sua elasticidade, alonga- 
se. Pode-se desenhar um equivalente mecânico de translação com uma mola K (elasticidade), 


atrito viscoso R, massa M. A força aplicada é u=Mg 


M y+By+Ky=u 


Suponhamos, para fins ilustrativos, que M=2, B=3, K=1. Substituindo obtém-se 


2y+3y+ly=u 


Para um dado u conhecido, por exemplo um degrau unitário, qual será a evolução da posição 


muscular y(t) ? A resposta obtém-se resolvendo a equação diferencial. 


| ut? 


Figura 3.2.1. Para uma entrada em degrau com será a saída ? 


2y+3y+ly=u 


Ma ut 


Aplicando o a propriedade da derivação no domínio temporal (ver Apêndice), obtém-se 


sucessivamente 
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a(nYio)-n dle) (o) +3 La Yto)- y66)] + Yin) = Vlo) 
(20"+30 +) Yo) = 290) + 29/0) 4390) + Ulo) 
vlo) = (22427900) tez gt) a 4 Ul») 


ZA" +30 +1 204"+30+414 
ni my 
dada: Yol=-4 Fa fas 
G(6)=L 


à 2 
x aii as 4 Mo 
2024304 (2026304) 2 


Verifica-se que a saída Y(s) tem duas parcelas: 
1) a devida às condições iniciais, Y7i, e entrada nula (zi de zero input) 
11) a devida é entrada, Y7s, e condições iniciais nulas (zs de zero (initial) state) 


Considerando as duas em simultâneo, depois de reduzir ao mesmo denominador, 


Yo) = ZAP — Atas 
Alenteso) A(a+i)(aros)da 


Procurando as raízes do denominador e decompondo-o em fracções parciais, obtém se 
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25"430+4 = (4 qin+t01)2 = (244) (o +07) 2 
«x e em pm - À 
As a=trtViria «br Lo 


2 
Yo) = o E cia A 


A (041) (0495) 


Calculando agora os resíduos Al, A2 e A3 (pela técnica do Apêndice), 


À, = Yo) = caoin+9 - 4 
p5o (0H)(0+01) 
As =(2+] (O) - a(nP- ostras — =0 
pe (=) (-1+0,17) p 
dc qi « 2snartat MO BE od 
= 


-0j (S5N-9 81) QU 


obtém-se 


Calculando agora a transformada inversa obtém-se y(t) 
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Verifica-se que y(0)=-1, como era de esperar. 
3.3. Polinómio característico, modos e estabilidade 


Se calcularmos a resposta da suspensão do exemplo anterior, quando apenas as condições 


iniciais são não nulas (não aplicamos qualquer entrada), obtemos, 


Yty) = (2443) fls) +2 yo) 29/0) 4 + [3 Yo) + 2966), 


e 24,84 4d 212405) (044) 


e invertendo, decompondo previamente em fracções parciais, 


Ay Aa 


et 
A+0,5 A+ 


-05€ ' 
Ja: (0) = ME” + ha a 


Yeda 


Calculando Al e A2, 


A, = UAI ae] III ado 


2 (2+1) A <-0,8 4 
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A = 216)4 + Days) + 29(0)1 Lt) to) +29l0) Oy 
TT  ——— neem TE 


2(ntes) AZ- = 


podemos constatar os seguintes factos: 


o a” 
Pat) se Soma da a. de 4 frisa 
for A, £ A? 


“ A, L Às duprsa dan dos Cd. patas 


E por isso podemos afirmar que as duas exponenciais 
e 


são características intrínsecas do sistema, independentes de toda a realidade externa e das 
condições iniciais. Chamam-se por isso, aos seus expoentes -0,5 e -1, modos do sistema e 


são como que o seu DNA. 


Analisando mais em detalhe, obtemos 
-8,5t - 
oct) = 290 +969] 0º" 4 [-gto)-a go] 


Para certas combinações particulares de condições iniciais o sistema exibe um comportamento 


estranho. Vejamos. 


Se 
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() =-Y(s) 
J61= -9(0), 


fact) = Eq) cet) 5 x gt) EE 


a resposta só depende de um dos modos; o outro modo não aparece, não é excitado. 


fe) = -2 lo), 
-g t cú 
Tlél= 2 plo ag) e = gh) 8 


Peça 
a 


a resposta também só depende de um dos modos, não sendo o outro excitado. 
Podemos assim constatar que, 

e Os modos são excitados pelas condições iniciais. 

e Existem condições iniciais especiais que excitam apenas alguns dos modos. 


e Para um sistema de ordem 2 existem 2 modos. Para um sistema de ordem n existem n 


modos. 


Qual a origem dos modos? 
Se no exemplo da suspensão automóvel do Cap. 2 (2.13.3) considerarmos diferentes valores 
dos seus parâmetros, a que corresponde uma constituição física diferente, por exemplo, M=2, 


B=4, K=20, u=0, t>0, vem 
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Hy-Birkg=u t=p fetdroga4 =0 


A*Ya)- A Ylo) = Yo) + 20T/9) - 2Ylo)+toTh) =» 
Uº tas ero) ro) = AY+ Jlo)+29 (o) 


Vo) = Le) +9e) ee9/o) 


2 
A +22+40 


Calculemos a resposta temporal, invertendo este Y(s): 


Yo). LS go) s2 900) (a 9ho) fo) +(%o 49, ) 


far)2+ q (a +1)* + 3º 
Desenvolvendo, 


— Jola+1) p Do s%o) 
(+14 )" +32 (A 41 + 3º 


A) cd “+ A s+Yfo) is 
(a +41) +32 54 poa) + 32 


Ou seja, consultando a tabela da transformada de Laplace, 


HE) = Yho) . Em 20 + L (Goto) E dub 


e graficando, 
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los | "q t 
oo Lgt)= Cestos comat 


1 envio 0 


a 


auspiraol Aommrara ad osciltas | 


obtemos o movimento oscilatório da suspensão. 


Se aumentarmos B de 4 para 14 (maior atrito, maior amortecimento) a suspensão deixa de 
oscilar (verifique, refazendo os cálculos, ou use o Simulink, com condições iniciais nos 


integradores). 
Portanto os modos derivam das propriedades construtivas do sistema. 


No caso geral temos, sendo p o operador de derivação, 


w) + j : j 
q (t) + END Pe e, y(el Eta LEI = Bout) 
ca 2 


(m mm = 1) é 
= Emultis ba (He. + 5, ut) boult) 


Aplicando a transformada de Laplace (a ambos os lados da equação) 
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Fa 40) - TA) Ata) - To 


ta, Ao) . Pulo) aa fe) 


+ab ol y6)] + Gollhn) = 


m nao! “2 pm =! 
= bm A Ulo) - bm Ato) 0 Cla) AT) 


ici 2) E la 
+5,49 Uta) * bm A” E a uilo)- bo d h) 


+ bolvlo) - db, Ulo) + bo Uls) 


ou seja, rearranjando, 


(4 


n n-1 
(» + And + +44 +40) Yla )+ [srito) 236) 


- | m mo ml 
> ( bmA + bm? 4.4 A 4bo) Ulo) + [-by, 4 ulo) .. =] 
e se a entrada u(t) for nula para todo o t 


Va) = ATO). = aço) 


zL n nA 
a ... 
ção E 


verifica-se que Yzi(s), e portanto yzi(t), depende das condições iniciais 
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3.4. Estabilidade em relação às condições iniciais 


Se um modo do sistema tem parte real positiva, a sua exponencial tende para infinito. Dada 
uma condição inicial que excite esse modo, a saída do sistema tende por isso para infinito, 
ainda que a condição inicial seja muito pequena. De facto, a soma de uma exponencial 


crescente com exponenciais decrescentes tende para infinito. 


Se um modo do sistema é nulo, ele contribui para a soma com a exponencial de zero, ou seja, 


com um termo constante. 


Se existirem dois modos nulos, poderíamos pensar que cada um contribuiria com um termo 
constante e portanto a soma das suas contribuições seria um termo constante. Mas tal não é 


verdade. 


Temos aqui uma demonstração do efeito de sinergia: o todo é maior do que a soma das partes. 
A contribuição conjugada de dois modos nulos, ou seja, de um modo nulo duplo, é um termo 
que tende para infinito com o tempo. Veja-se no Apêndice a transformada inversa de um 


modo nulo 1/s?. 


No caso de modos puramente imaginários (com parte real nula) temos que considerar 
sinusoidais de amplitude constante, quando se trata apenas de um, ou de sinusoidais de 


amplitude crescente até ao infinito se se trata de um modo imaginário duplo. 


Sempre que uma dada condição inicial excita um modo duplo com parte real nula, a saída do 
sistema tenderá para infinito, mesmo que a entrada seja nula. O sistema diz-se então instável. 
No caso contrário, em que todos os modos têm parte real negativa o sistema é estável. Um 
modo de parte real nulo mas de multiplicidade 1 mantém a estabilidade do sistema (em 
relação às condições iniciais): a saída não tenderá assimptoticamente para zero mas manter- 


se-á dentro de um certo intervalo finito. 


Assim podemos enunciar uma condição necessária de estabilidade em relação às condições 


Iniciais: 


e E condição necessária de estabilidade que todos os modos se situem no semi-plano 
complexo esquerdo, ou seja, tenham parte real negativa, ou, no limite, se situem no 


eixo imaginário com multiplicidade 1. 
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Note-se que os modos imaginários aparecem aos pares conjugados: se existe por exemplo 2i, 
existirá também necessariamente -2i. Mas o modo é ainda de multiplicidade 1. Só no caso em 


que haja dois modos em 2i e outros dois em -2i é que o modo é duplo. 


A resposta às condições iniciais pode representar-se por 


em que 
Aln) = p"+ E +a,4+4 
T(n) = Ba.,9 + MR di + Bro+Po 
Po : dlejumes da, (onda pmias, Yo) 9h): 
Chama-se a 


A (pn) = poliweimio conaderadis de sesama 


dutumine q peoproda Levar, 
eu vas focada qu uatiral. de Silinma 


raio da fla) 2 mods Às ssiima 


A resposta livre, ou não forçada, é a que se obtém com entrada nula (e condições iniciais não 


nulas). 


3.5 Função de transferência 


Seja o exemplo 
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as N 
2906) +30) +94) = ult) 
Calculando a transformada de Laplace, com todas as condições iniciais nulas, 


(202430 +)Yhl=vlp) 


Tlol= 1 Ulo) 
23 A+! 


Ou seja 


Yo) = G (0) Uta) 
y 
FUNICAS DE TRANFERÊNCIA 


da sanida e de minado dra 
Cd: tie Aviz vuls. 


G(s) “transfere” a entrada para a saída, ou de outro modo, a entrada transfere-se para a saída 


através ou pela acção de G(s) e por isso chama-se a G(s) função de transferência. 


A função de transferência obtém-se a partir da equação diferencial de ordem n que descreve o 


sistema, aplicando-lhe a transformada de Laplace e resolvendo para Y/U. 
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Se tivermos a possibilidade de aplicar um impulso de Dirac ao sistema, a sua resposta será 
Y(s)=G(s)U(s)=G(s).1=G(s) 
e portanto 
HO = L'[Y(s)]= L'[G(s)] 


ou seja, a função de transferência também se pode definir como a transformada de Laplace da 
resposta do sistema a um impulso de Dirac. Por isso se conhecermos a função de transferência 


e quisermos calcular a resposta impulsional h(t), basta fazermos 
At)= [et 20660) ]=L"[660).1]=L Ven] 


Por exemplo 


| JE a 
Gl). 2044 414,42 sMteta2e 


ACaên st A+4 As2 


Note-se que a saída tem duas exponenciais, cada uma delas produzida por uma das fracções 
parciais. E é sempre assim: a resposta impulsional é a soma das contribuições de cada fracção 
parcial da função de transferência. Se forem fracções simples, resulta numa soma de 


exponenciais. 
Exemplos de cálculo da função de transferência: 


e sistema de dois tanques (ingestão de um fármaco no Cap. 2), com Aj=Az=1) 
.. “ 
A ha fito Do + (ABr Mola Ao) É, 49 270 
| A,hr (4 (47 + (A, R, +Az Lr+A, R2) +1)0/4)= Qulo) 


ln) Fi Qeln) e” 4 
Qulp) AAA R, A +A (Ay +hz Ra +A, 2) 
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0 Sucjumado andomsvel Mq+By+ky=u 


(nem +nB + k)Flod= Vln) 


= MD). Ao 
Vl») 420408 +K 


No caso geral, sendo p o operador de derivação, 


0 Caso qual A(p)ylt) = Blb) Ult) 
G (A ) = Bla) 
Ala) 
Polos e zeros da função de transferência 
Os polos de G(s) são as raízes do seu denominador 


Os zeros são as raízes do seu numerador. 


G(p)= Bl) 0(n)>0o —> 2ern 


Ala) fla)>o —> poho 
Se Gl) 2” probris( au Az an 8) Borin 
Polo. Um uimito Pit À», Atad ua Copias 


eum fuel da Gl) se Gfh)=00 


Bero. Um vmar fin natal 84 
Conflito e! dm 31mm do Gn) se A (a)=0 
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A função de transferência pode escrever-se com o numerador e o denominador factorizados, 


dizendo-se que está na forma polo-zero. 


£,,2u,--, Em » ero 


Pr, Pe, e, pn g pelo E ia 


cla) « +(0-8) o-êe) (0-2) 
(A-br) (o — foz) -- (0-f) 


A função de transferência é assim completamente caracterizada por 
- os polos 
- OS zeros 


- o ganho K, no numerador, uma constante característica do sistema. 
N(n] =o => (A-% Há Be) ce (o-tm)=o 
EN (p) >o x (0-8) (0-22)--: (2 -&m)= ” A reje 


o 2200 fm sm no Lofirim Lero plonra Io o yuamasasola, 
Se por exemplo tivermos uma função de transferência G(s) com 
e Zeros: -2, -3 
e polos: -1,-0.5, 4 
e G(0)-5 


Dos polos e zeros poderemos escrever 
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G(o) = Me (nst)lo+3) K? 
(a + (a 49,7) (ou) " di 


Para calcularmo o valor de K usa-se o terceiro dado, obtendo-se 


Go) = 5 G/o)=S = ia > 3K 
Lo s.y 
=) K = (2 


% (n42)(043) 
[044 (0495) (044) 


elol= 


Qual é o significado físico de um zero de uma função de transferência? 


Se aplicarmos ao sistema uma entrada cuja transformada de Laplace U(s) tem um polo no 
denominador igual a um zero da função de transferência G(s), quando multiplicamos 
G(s)U(s), esse polo vai cortar com o zero de G(s) e o seu efeito na saída fica assim anulado. 
Se por acaso a entrada tem só esse polo, ela não influenciará a saída. 
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3.6. Cálculo da saída, conhecida a função de transferência 


Conhecida a função de transferência de um sistema, definida uma entrada para ele (no 


domínio complexo), a saída no domínio temporal calcula-se simplesmente por 


Rd EI pa 


Ytn)= €la).U(a) 
y68) <L [tt] euade ca cond 


(UA ls Vas 


Vejamos um exemplo, de um sistema com uma entrada em degrau unitário. 


G (n)= 


(240, 5) (041) 
ulo= Llgin]= Lo daapom umltaizo 


Tola & (o). Vin) = pi 
ETR VE 


“ A3 
P+49j 94 él PAGE A+ Pe) 


Os resíduos foram calculados por 
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Aq = ato] > z-2 


2 =-04 2(0+1) 


» 0,5" 
AZ feto Yo | E 4/2 - 1 
di (p+497) 3 Io <= 
hs - A Vo | - “fa = 4 
p= (2 +05) (241) 


e finalmente, consultando a tabela de transformadas, 


no -0ST 
Jia = TO a O ca 


Vejamos as diversas componentes de y(t) 


iba q sao do 
aa É dind om der fo da salame 
wubade S (9) 


Poderemos extrair a conclusão de que a forma da resposta com condições iniciais nulas e 


entrada U(s) não nula depende de: 
e polos de G(s) 


e polos de U(s) 


(DADC/FCTUC/2022 Licenciatura em Engenharia Biomédica/MCPF 


97 


Dinâmica dos Sistemas Biológicos e Fisiológicos Cap. 3 Função de transferência e dinâmica de sistemas 


Exem ls (Chen, 184) 
Glp) = (ato) (ama) (4 -1)* 


A [a s2) (4 +3) (n+42 -[2) (» +2 t+/2) 


pas ayalpu edrsde, Ulo), 
bend À . k É 
ha, (2+4]'2) 


Ylt)= Ay + trt, Ay E a A+ 


— a = 
A Sn at th, Etna 


+ tamo caviar U (4) 


— 4, ) Uta) mos tim plo fra a gerw de Gl») 
Vulrma pelo /€3 trás Cama GA) 


Uls) = (a +3) 
(141) (042,5) 


ma durdo ao fude ela Vl») 
As 2 + Ay Pc 
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— 44) Uln) mos tim polo Aprois & zero do (ln) 


Uln)= fp=3) 


A (042) 
(A410) (nm2) fá AJ (5-3) 
E (ne) lo 43) [24 E-8) (+72) 


Yhl= Gn) U(o) = 
AUAMMAdaÃ BA 4, 
ira E Hide 
qo 
nt 


q 
-20 - 
se ita 2 ap 22 


YM= 4 +hzt esto A 
+ kg 2 snat 4 hq E do 


a 
lima doando aU(o) : me, he te 


NE, Un) Ama fudoo a prai « Zero de 6 ln) 
2 nas formais Cm 6(n) 


dia) + Sea 


(A +10) nmez) 


2-3 


(A +40) (042) (2-1) - 
pH n+2)(043) ([n42-12) (n 4242) fon) 


Ylo)= 65). Uln)= 


trmo devida x Uln): AS 
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iv) U(n) tim gera apurar a Judo de 
G(n) emos tm pio comun cj 60) 


Ulaja 243 
(+40) (044) 


tnl= Elo) Ulo) = fa 410) puz ) (0 É “2 
nave) (248) [042 -J2 Ma segja) (241) (044) 
Us) Jaru E te 
ama Êo pol -3 nos eueilado 


Chama-se cancelamento polo-zero ao que acontece em 1ii) e iv). O cancelamento elimina da 
resposta (do sistema) características intrínsecas do sistema, e por isso é um acontecimento que 


não deve acontecer. 
Através deste exemplo poderemos generalizar e concluir que: 


e Quando U(s) tem um polo igual a um outro da função de transferência, resulta em Y(s) 
um polo duplo. Este polo duplo introduz um efeito próprio na resposta temporal 


devido à sua duplicidade. 


e Quando U(s) tem um zero igual a um polo de G(s), então em Y(s) dá-se um 


cancelamento de um polo por um zero que elimina o efeito desse polo. 


e De modo semelhante quando G(s) tem um zero igual a um polo de U(s) dá-se também 


um cancelamento que elimina o efeito desse polo na resposta temporal. 


Como vimos anteriormente, a resposta temporal com entrada nula e condições iniciais não 


nulas é determinada pelos modos do polinómio característico. 


Considere-se o exemplo 
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Ya) = (Coto tag), 4 ul») 


LAº+ 30 +1 205º+ 30 +1 
(ço 
fer Tas 


protinanis ponadaridis : 2543404 
fracas du transfira (Yo) =9(0) =6) 
Clo)j= A 


20243044 


dunmeinada Dio): 24º43041 
e rpral om pdivivas paadicdir dep 
o peto pos A prsis dm ido, 
Vejamos outro exemplo 
24 + 3y do = uL+u 
(2.42 +34+1)Y(0)- A2y(o) ap (o) = 39h) =(A +13) Ulo ) - 4/0) 
vt). 2H) +2Goezylo) alo) ar qu) 


Z2NA2+304+1 2042 0+) 
E 
RA 1 
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- 22 Glo)+29le)+39h)-Wlo) 4 AÁ Ulo) 


Z(o+)(2+0,5) 2(24)(2-+95 


4 Me Vl» 
AGE 


Pine pardindo: 2024304) 
pure : =4,-95 


dameninasa da Gl) : Atos 
polo =0,00 
Assim: 
e todos os polos são modos, mas nem sempre todos os modos são polos 


e se há um cancelamento polo-zero no cálculo da função de transferência, há um 
modo que é eliminado, e a função de transferência deixa de ser uma 


representação completa do sistema. 


3.7 Resposta temporal dos polos e modos em função da sua localização no plano 


complexo. 


O plano complexo pode-se representar pela figura seguinte. 
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Ro 


peer] À À jocmeme == 
0 — = 
> U)-— A pescar 


pmapindria 


O SPE é o Semi-Plano Esquerdo e o SPD o Semi-Plano Direito. Qualquer polo e qualquer 
zero se situam neste plano. Como sabemos os polos podem ter natureza diversa. Eles podem 


ser 


Recs Tma ponairivo” Com plitio” 
Semplo, HErmadols, Sverjrha, 
Hu bio, hu lkALo múltho 
"am fare conjugado 
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Vejamos o que acontece nessas diversas situações. 


D) Polo ral Simas 


Uls clos cms qul= LE 


Tie do —s qtg = 128 


Fido ão ca Ud 


A =1 t.0€ 


— 


Ya) = à —» qt) al Fa 


J“ 


- ht A Z T 
RR 
lim yé)a O 


lim 4 (tl) =00 
€ -s00 


t-s00 
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1) - Polo red mulkplo 


Yta) = 4 é Rs (W= to 
(n+1) 


Ta) = z a EE (e) = E 
+ 


Yo)jz=+A Ss ye jo 


(2-1)? 


Yoda Lo sys tott , 
(p=)? 


Conclusão: 


e os polos reais que se situam no SPE contribuem com termos decrescentes 


e os polos reais que se situam no SPD contribuem com termos crescentes 
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na Pole na or tem 
>implu 
Tips = e” ylt)=t À 
múltoplos 
Tale Alis ut) = É EM 
A já 
no a 42 / 
Th) = 75 cs YOE) ..€ . 
Conclusdy 


Polo na orcgtm 
simpl, - contrcbu com um termo condante 
multthlo - contribw' tom um tirmo cresanie 
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iv) Polos imaginário confugado 
Siva pita 


Va) a í ——» y(£) = sen. É 
A" +41 


(= Lo + yit)=bmat 
A + 4 


conclusás 
o polos úmashamio sino, cortrebuam cum 
sms das duran 


Mútt 
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v) Pole complxm dompugadr» (simples) 
- multiplicam o alerts das reais falo efeds 


Amaanaarn 
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Sintetizando todos os casos que vimos, obtemos a figura seguinte. 


RESUMO Dk LOCALUIZACÃO Do! POLOS E Dos! moDos Am 
' 
Plano S E suas Cow TRI BuI Cos TEmPormis 


3.8. Estabilidade BIBO- Bounded Input Bounded Output 


Um sistema diz-se BIBO- estável se a uma entrada u(t) limitada (BI-Bounded Input) 


corresponde uma saída limitada (BO-Bounded Output) 


No geral vimos que 
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Ya lt) = ( Aavida aafdo 4 fo) 
+ ( Avda om fdo d vo)! 


Se os polos de U(s) têm parte real negativa, ou, sendo imaginários puros têm multiplicidade 


um, a entrada u(t) é limitada e a sua contribuição para a saída é também limitada. 


Se os polos de G(s) têm parte real negativa (com qualquer multiplicidade), a sua contribuição 


para a saída extingue-se com o tempo. 


Assim 


É condiços necssarnz e guficial elo 
esta bobilada BiBo qu nm pato A Elo) 
AR palina no da 


Considerando a relação entre os modos e os polos, poderemos dizer que 


e é condição necessária e suficiente de estabilidade BIBO que os modos se 


situem no SPE. 


Quando há um cancelamento, isso não quer dizer que o modo do sistema cancelado 
desaparece. Acontece simplesmente que ele deixa de ser visível na função de transferência, 


mas o seu efeito permanece. 
3.9 Influência dos zeros na resposta temporal 


Sejam o casos seguintes e as respetivas respostas temporais 
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A 4 -0" —bê 
Mo pra e q (EA 
Ulaja MS aj [ecoa mat] 
(Ata )(n+h) dé 
O zero altuws 0x Letfiwewiis 
da, Le ponenudi, 


Se por exemplo a=1, b=2, sem o zero teremos 


UM) = Fal ge 


Com o zero será 


| Em 
X=0,8 = JL) = raid +48 L 


d=3 > yt)= gg" qr 


ytt) é uma combinaca bina ole eeporenudis 
eos prolos dofinem & ma do eeboreuuiai 
(acgas Jatina ) 


o 68 Cert defivum 05 coefriuenko dessa 
Comp ts 
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3.10. Resposta temporal a um degrau dos sistemas lineares e invariantes 


) 
U(a) tta 


com função de transferência 


Seja o sistema 


G(A) = Bl) — bm ne aa es +biA+bo 
Alo) Me a, AP A 4a A+ 


2 (0-2) (A-20) -(A-2m) 
(A -pola-pe).- (o-py) 


man 


Consideremos o caso 


e men, Gn) etita mea bropri 


o todo À pelo pe miuam we SPE 


que é suficientemente geral para abranger todos os casos práticos. 


A ordem do sistema é n, e considere-se uma entrada em degrau unitário. 
ordem de sidima = n, Gram dr P(A). 


Uls) = 1: mu6 attl pars conadira 
A 


— ssima,. 
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1) Adema da 1º ordem 


G (5) = aaa e K Ep 


(aso) 
donde, 
Tlol= ss MK. aa Az 
nata A A P+A 

Calculando, 

AM = Ao) “a = 

ATO PPA o 
K 
ho = (a teto) s Fr =. É 
AZ-A a AZ-a & 

obtém-se 


e graficando observam-se as características gráficas (ganho, constante de tempo) 
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valor fival (em Taco) 


Plana s 


e Ganhs + 
e constant de Tempo te 


1)  sutimas da 2º ordem 


G (n) = - 


Ara s+ão 


polm da G(n) em furto de 2,246: 
À & + VaZ-4ao 


RÉ OT 

Neste caso temos várias possibilidades para o tipo de polos: 

» VGA —2 poli reas simbo 

; 2" =-44, ada À pot past duplo 

à Mtlhas— 2 polo complio Cojupdo 


E frequente, por ser mais conveniente, usar a seguinte notação para sistemas de 2º ordem: 
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G(n) = k 
pr+ 23 VA + 2 


Asa = —S Um El UV $24 
e Z>1 > 2 pos rearr simbla, 
, 4=4 51 pol ral dujlo 


coc4 <1 > 2 hdo complio tmpudo 
, Aya == *) tomVt- 3 


O tipo de raízes depende do parâmetro & 
4) 431, 2 polo reais ditas 


At — — É, W, - 4 wu (42-14 
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Exemplo 


G(n] = k = E wn=V2 
AÉ+3A+2 P+D(0+2) Go A 
PA 
LUlo)= = 
Vinte tis q a E PE da, 
(o4)lo42) 0 00 At Apa 
a E rag 
a 9 A+] A+2 
e ate 
q(t= o us K LA & a yll=o 
2 k 
geo) = > 
Graficando, 
Plawo » > 
Jo 
pé pomk ou impluxas 
+ 
gare dana, 
e Gar K/2 
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VA) Gal, uma pio poal duplo 
ExewatÃo 


Gn) = K = PM tô, al 
pÉ+204  (ax4 E É L 


Por exemplo, 


Yo) = E MR é A 
p+1)> A A (0+1) A+1 
ax K k 


Jt= K- et oo nã” a k [4 SE] 


ylo)= [1] 
Y (00) = k 


cujo aspecto gráfico é o seguinte. 


y(t) 


Plana 


F 


es retÃipaÃoan 
e Goro e, 
4 froma ola cursos  dipumots de tu, 
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tusj0c 4 <1t, 2 polo complixo conjugado 


Ed 
Mr 2fi, A+ tum 


Gl») = 


A saída, para uma entrada em degrau unitário virá, 


K A 
t . A 
A +23 ba +44, 


Yo) = 


-fut 
ES L e = 


JM) = Ko. 
wu,2 tu Vi-sE 


= cas” 'q 


Graficamente teremos 


Piaom t Ef a Pre DE 


|fe | = 


qto) = o 


de! y00) = a 
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Juntando agora os três casos numa só figura, 


2 >2.2> 
G(n) A aa ca ] 
428 UA 42 ia ldado 


$<1 


A resposta y(t) classifica-se como 


oc Gar: esposo sub amerteuda 
Ê =1 rechrdo Lom «metêcmek cribiur 


421 É rh sabre -am tectos 


o use 


44 cep dh ama - 


E ha 
Sta = — Sum | tw V1-42 4 fuga la oscila 
= — 2] Ea” d 
cia af adç 
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O coeficiente de amortecimento e a frequência natural são propriedades intrínsecas do 


sistema. 


Representando no plano complexo, won é o módulo das raízes complexas e & o cosseno do 


ângulo com o eixo real negativo (ver figura seguinte). 


Ju 
4= Va. $=0 » Md = in 


ug dewmiya 


D) 
WA = um VI-$2 


A frequência natural amortecida depende desses dois parâmetros. Como se vê na figura, se & 


aumenta, os polos aproximam-se do eixo real e oa diminui. 
A frequência natural também se chama frequência de ressonância do sistema. 


Há várias características da resposta temporal de sistemas de 2º ordem que têm interesse 


prático, e são ilustradas na figura seguinte. 
Por exemplo o tempo de subida exprime a rapidez do sistema a responder a uma entrada. 


O tempo de estabelecimento mede o tempo que o sistema leva para alcançar um regime final 


(constante) quando a entrada é um degrau. 


A sobre-elevação mede uma característica frequente de sistemas dinâmicos: a resposta 


ultrapassa o valor final antes de se estabelecer. 
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é 
sá ty 
Gl) =: - " , 


e 


Caradenttio, princpars 4º, wm 

T= gun) constant du temp 

T= É0g. tempo de pubida (do o aqi) 
e: Rs timpo da furos 

PO = 100 wp(çEs )sobreceliuacod em Yo (fear ovushagt) 


To = 3T (só) fâmpo au estrblecimenk até” Yo evro 
yT (2h) 


Efeito da introdução de um zero: 


Se introduzirmos um zero no sistema de 2º ordem anterior, teremos 


A >o(em) 
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O zero está expresso como um múltiplo de &wn, a parte real dos polos, por conveniência de 
análise. Fazendo simulações para diversos valores de a, obtêm-se os gráficos da figura 


seguinte, para o caso de €=0,5. 


[v) 2 4 é g lo tt 


O zero altera fundamentalmente a sobre-elevação e o tempo de subida. Quanto mais influente 


é, menor é o tempo de subida, mais rápido é o sistema) e maior a sobre-elevação. 
O zero tem assim um efeito de aceleração. 


A localização do zero no plano complexo determina a intensidade desse efeito. 


Plays De) 


A Us, 


diminui a sua influência 


Quanto mais longe dos polos estiver o zero (em direção a -c0), menor é a sua influência. 


E se se introduzir um polo adicional, passando-se a um sistema de 3º ordem ? 
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3.11. Sistemas de 3º ordem 


Um sistema de 3º ordem tem três polos e pode obter-se de um de 2º com um polo adicional. 
Para manter G(0)=1, a fim de simplificar a graficação com a saída a tender para 1, como nos 


casos anteriores, seja G(s) 


G(n) ” x 4 tu) 
(A+ ag ums (+ 248,0 +02) 


Para uma entrada em degrau, e simulando para vários valores de a, obtêm-se as curvas da 


figura seguinte. 


[2 
((B 
, Pa 
As 
04 X=0,1 
LIXA 
[RA 


o 2 4 E &º tow 


O terceiro polo aumenta o tempo de subida, a resposta fica mais lenta. Quanto mais próximo 


estiver da origem, mais influente é. Em -c0 não tem influência. 


A posição do terceiro polo depende do valor de o. 
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Quanto mais à esquerda estiver o polo, menos influente se torna. 


3.12. Sistemas de 4º ordem e superior 
Do que vimos anteriormente podem-se enunciar os seguintes princípios gerais: 


e um zero adicional acelera a resposta e aumenta a sobre-elevação (torna o 


sistema “mais leve”, “mais ágil”) 


, 


e um polo adicional trava a resposta e diminui a sobre-elevação (torna o sistema 


22 66 


“mais pesado”, “menos ágil”). 


Graficamente, teremos os seguintes efeitos sobre a resposta temporal, para uma entrada em 


degrau: 


, t 
a et PA Gls) k E (o+3:) 


“(esp (246) ovbo [espe oah) 
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3.13. Sistemas de fase não-mínima (zeros no SPD) 


Para que um sistema seja estável, não pode ter polos no Semi-Plano Direito, e por isso se diz 


que o SPD é a região de instabilidade. 


Os zeros não influenciam a estabilidade e podem situar-se, desse ponto de vista, em qualquer 
ponto do plano complexo. No entanto influenciam o regime transitório da resposta. Podemos 


dizer que os zeros não influenciam o longo prazo, mas apenas o curto-prazo. 
Vejamos um exemplo de uma função de transferência com um zero no SPD. Sejam 


À o s-2 " E —s—2) ad 
Ve ane O a PS 


Note-se que as funções de transferência têm ganhos de sinal contrário. 


Para a primeira teremos 


is ÃO... mus sd sã 
pls 41) (42) dá PA DM A 04 ne 


gltla ma 4302 0 


Para a segunda-feira, e dado o sinal contrário de G(s), teremos uma saída simétrica. 


Glnl= gs ge (x==1) 


(pAn)(lo+2) 
La) = = lo=t) 
p (n+1) (n41) 


yíl= 24-37" 0 


Graficando, obtém-se a figura seguinte, que ilustra o comportamento típico dos sistemas 


chamados de fase não-mínima ou de resposta inversa (porque a resposta inicia-se em 


direção contrária à final). 
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3.14. Sistemas com atraso puro 


A equação diferencial 
qllragit)zbu(t-t) 


diz-nos que a entrada do instante t-ti influencia a saída no instante t. A entrada necessita de 
um intervalo de tempo ti para ter efeito sobre a saída. Diz-se que ti é um atraso puro (pure 


time delay). 


O diagrama de simulação desta equação diferencial é 


Aplicando a transformada de Laplace à equação diferencial, tendo em atenção a propriedade 


de deslocamento temporal, e invertendo, obtém se 
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(442 )Ylo) = 4(0) + ba Up) 
plel= LU [SOIL Lap 


Ou seja, a resposta às condições iniciais não é afectada pelo atraso (como era de prever, 
atendendo ao diagrama de simulação, dado que elas entram directamente no integrador). Mas 


a resposta à entrada é afectada pelo atraso. Considerando-a separadamente, 


G(»)- Yto) = b At 
U(p) E Ay Pã tinal 


obtém-se uma função não-racional, e, portanto, para ela não é possível calcular a 


transformada inversa de Laplace. 
Com calcular y(t) nestas circunstâncias adversas ? 


Tomemos o seguinte dois sistemas em paralelo, com a mesma entrada, 


Na parte de baixo, não havendo atraso puro, não há qualquer problema em calcular y(t). Mas 


na parte de cima aparece o termo irracional do atraso puro. 
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Mas haverá alguma relação entre a saída de cima e a de baixo? Se houver, e se a 
conhecermos, poderemos obter a resposta de cima desde que tenhamos obtido a resposta de 


baixo. 


Ora atendendo à natureza do atraso puro, um efeito que apenas retarda as coisas, poderemos 


concluir que a relação entre as duas entradas é ilustrada pela figura seguinte. 


4 


afreio 


puro 


Em geral, para a mesma entrada U(s) 


Exemplo de um sistema com atraso puro (Murray, Mathematical Biology, 42). 


Há uma teoria segundo a qual a concentração de dióxido de carbono no sangue controla os 
níveis de respiração através de uma malha de feedback com atraso puro. Pode-se adoptar um 


modelo simples, adimensional, na seguinte forma 


e OE v(g) = 
dt l+x” 


Sendo a e m constantes positivas e Ta o tempo de atraso puro. 
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Em Murray encontram-se modelos populacionais com atraso puro. 
3.15. Resposta em regime final 


A saída em regime final, ou o valor final da saída, é definido por, 


€ seo 


No caso em que u(t) é um degrau unitário, 


Ylol= 66). Ulo)= Glo). L 


Aplicando o teorema do valor final (note-se que tal só é válido se G(s) for estável), 


Teo = hum yt)= Lim A Yla) 


tao AO 
=ltm AL. 6(n) = Um é (4) 
A Do dá li 
= 6 (n)| -— 6/0) 
A zo 


Loo = 6 (0) 


Tendo G(s), para calcular a resposta em regime final para uma entrada em degrau unitário, 
basta substituir s por zero em G(s) e fazer as contas. O valor constante obtido é o ganho do 


sistema, isto é, o valor final da saída quando a entrada é um degrau unitário. 
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3.16. Síntese 
s e (n-)) ” (m> (m=1) 
q + Qn=1] eee + ay+ão = bm U + bt ++ bolo 


cha 


(9 [condiçoes Lnauors = Dia É E A (on «ns 


W ne! w ne! 
A + 04.44 + ft As A+LA, A esses +&o 


ac Yes 
funcço ou trauckrimuz 


DRE na jfro 
Gh)= bm pe ««tbo sé bm (o-20) 00-29). (om) 


Praia (BIO) (abr) 
ia A 


Polinimio faradiifve 
Ne A Re +AA+RO (de Vac) 


Modn : rouzs el pin dana dani 


em relacos a Px 
Etabilideds : esta e Cdr o mods estar no sPE 
BIBO cstovel se tel o polo ez: wa SPE 
Flo 2 made, atais au Completo Crgupados ,vrimpls 
Qa ruas Hi lho. 
fovabrinate, duiuar da, Lreha netas dra 


Pesjurda sas ALfira Fral fome usa um depor ; Elo) 
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RESPOSTA A DEGRAU 
(analhm gua lilalane) 


LUNÇHO DE TRANSFERÊNCR pesfosta (t) 


EA 
* A+a 
4 


bolo rara Eunpho, am 


pol red dubl 
pola comfl. Conpugaa mm 
furl ma SPD 


Gln)= Gn). 6º 


ATERRO Pulo | FASE NãO MINI EFEITO Dos -Z ERAS 
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3.17. Vantagens e limites da função de transferência. 


A função de transferência é um conceito e uma ferramenta muito útil em estudos de dinâmica 
de sistemas. E fácil com ela calcular a resposta do sistema a qualquer entrada, incluindo a 


entrada impulsional (impulso de Dirac, cuja Transformada de Laplace é simplesmnetel. 
Com ela pode-se descrever matematicamente a composição de subsistemas. Por exemplo: 


1) sistemas em série: produto das funções de transferência 


vb sra Eb 


G(s)= 6/( ade G, k ) 


11) sistemas em paralelo: soma das funções de transferência 


Gln)= 6,0) + Go. (0) 


Hi) sistemas com retroacção negativa 
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iv) sistemas com retroacção positiva 


Esta forma de representação tem também sérias limitações. 


A mais grave resulta do facto de só se poder aplicar a sistemas lineares. 


Quando existem cancelamentos no seu cálculo, como vimos anteriormente, ela não é capaz de 


representar completamente todas as características relevantes do sistema; diz-se então que é 


uma representação incompleta. 
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Capítulo 3 Apêndice 
Transformada de Laplace 


3A.1. Introdução 


No capítulo 2 obtivemos equações diferenciais de diversas ordens, lineares ou não lineares, 


como modelos matemáticos de sistemas. 


) (m) 


o (rm=!) q 
o ssnih ,9+40]= bm 4 + bug +. tb +b,4U 


ra gr 


Ali cando puma tulicda wu qualpun, punto Yobbide! 
Lo resolucas de aqua dfrenoal 


Se as quisermos usar para prever o comportamento do sistema em determinadas situações, 


n) (ns 
+ & 1 
n-, 


teremos que encontrar a solução dessas equações. No caso linear, para ordens baixas, 
poderemos encontrar uma solução analítica (pela técnica da solução geral mais uma solução 
particular, por exemplo); mas se a ordem for elevada esta técnica pode tornar-se de difícil 
aplicação. Também para sistemas de equações de 1º ordem nem sempre é fácil aplicar esta 


técnica. Importa por isso uma técnica alternativa mais simples. 


No caso linear a Transformada de Laplace cumpre com êxito esse objectivo. Tem além disso 
a vantagem adicional de, aplicada ao sistema, permitir exprimir facilmente importantes 
propriedades do sistema (dinâmicas, estabilidade, etc.) que são postas em evidência pelos 


parâmetros da própria transformada de Laplace. 


TRANSFORMADA DE CAPLACE 4d 


Utilidade da Transformada de Laplace 
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No caso não-linear temos que recorrer a métodos numéricos para transformar as equações 
diferenciais em equações de diferenças e depois usar o computador para iterar. Trata-se dos 
métodos de integração numérica (rectangular, trapezoidal, Runge-Kutta, etc.) que estão 
implementados por muitos programas científicos, incluindo no MatLab/Simulink, como 


sabemos. 


Por isso neste Apêncice iremos reveros elementos da Transformada de Laplace úteis neste 


contexto e a sua aplicação à resolução de equações diferenciais. 


Consideraremos que todos os sinais temporais existem no intervalo [0, +00 ], sendo nulos à 


esquerda de zero. 


3A.2. Definição de transformada de Laplace (unilateral) 


A tunmtormada de Lopes de uma firmam Pt) 
dutwoda tm [º, +00[ 2 dada pr 


cus Lltle]e [Pte)-ett de 
A=t+aw E operado du Laplace 


Eu Sei fl) = EE Gude Flo) 


eo CO duo at] 
AX [ dé artes | 4! nc O tia 
É ó A-2 a 
j io -(n-2)t | 
= — L - 2 
A-2 | =00 t=o 
O a 
=4, ta 


-(n-27t (junto -(2)€ que 
4 =L L a -—- 


- Er (cos tot-j sen wê) 
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= (7-2) t 
L 


(cn E =p sem nm 


linulads Mas odindla 
tomar ma boda E1,1] 


) Com t t-»oo 
1) se Gl2' 
+x.00 


E = (x>0) 
at) Se [=2 


“Ea 


tua) 52 


Considerando as duas partes em conjunto, conclui-se que 


-(A-2)t to , x [<2 
=/ indupmds, Se MC=2 
O, se (>2 


fntamk a Lam po mada au Lapla ou dis 
funda posa Pes Clulands-a 


Fla)= =t [0-1] 1º pow 52 


A-2 
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34.3. Definição de transformada inversa de Laplace 


Dada uma função F(s),a transformada de Laplace inversa é 


E) Co 
os Lentes de antipal. definem PT SS MA 
O uxo Aral Ava q=c, red ad yu 


Em gpa star denho as du comun têmua 
e astim mas fr, Lda! puedo svlsrredos. 


Geralmeda nox Le AR à defina forr calulos 


tt) a A de Fls), cons vemmo, for mo 
a repos ou Comverpivia fruta Gama fa ola 
su qsnirsee «nã tc tim pamsnenta mm 
cia. 


C +00 + 
tt) E Leto) & de | Fl) e da fam €zo 
2 


3A.4. Transformada de Laplace de funções elementares 


Impulso de Dirac 


doe Í s(Ddt=1 
= e = 
0,t=0 a 
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AL 2 a 


Elo) = LEO] = [ 175 Rg a 


[4] 
(todo o SU) si! dudro dn lirahs ol imfipsem) 


t=o 


i) Degrau unitário 1 
GE Er=b 0 
TOC lo,e<o 


Fln) = 2 [aeil = [é 2 o. - [ Ea 


0 


-A00) 2º*| E E -1) E sa 


ii) Exponencial 


Exponencial > 


f(n) = E] =[ Dede [38% Eu 


Õ 
qt, Gta en [= 
DR TAa - 


o 


/ 


Ata 


n 


Se a<0, a exponencial é crescente, tendendo para infinito. 
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iii) 


Rampa unitária 


Rania mito tH)= E 


Ftp) / CÊ / Paes 


(o) 


1) 


SAD 


22 
giga oo o o 
4 SA E o au | 
4 A. ô õ A a 
= "Az 


Para funções mais complicadas ver a Tabela de Transformada de Laplace. 
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3A.4. Propriedades da transformada de Laplace 


i) Linearidade 


Flo)= Eq [t(8) É lp)= 2 (fot) 
Plegtore bm] = «hl)+ Bolo) , up das 
x. Í [onwt) = 


we vw) 
wo wyt = sé + e* 


AS CE) LA [Ta (e 


, j 


-! qdo ll qu DES 
' ) e Aê 2 
Au calor” God tw 
| 
= — 
a2p ut 
de mama fama, 
2 
2 [om vt] 
neo 


ii) 


derivação no domínio temporal 


Z[SHe]=2E0]=4Lge]-to) = 40) ho) 


“E gm) = 8Fla) 0" lo) nf) Ea) 
A 
dininan darwin Cponat pulhiblia 
ai de ist À ni pr Lia Ea 
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iii) integração no domínio temporal 


Lytirpacos we denmino (ndinal — destas pr 4 


UV tear] = 466) 


j do do 
0 Pá 
iv) deslocamento no domínio complexo 


Psloca med vor danaimao Crraflisto a ralada — 
cmo po à Eyr demminmo Mnpnal 


f [2 *ete)) = F(n+a) 


“> Z [e nen vt] = * Lonas] , A EIR 


AxAta 
ly 


+ 


. (n+a)Prco 


[2 cut). E [cos wt | 


A-A 


(A «aj, Lya 


AdDA-A 


HH 
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iv) Teorema do valor final 


Este teorema relaciona no comportamento de f(t) no infinito com o comportamento de 


F(s) na origem 


AA PÉ) = AA A Éla) 


t De N-Do (1x Fla) mes Em pola me SPD) 
v) Teorema do valor inicial 
Relaciona o comportamento de f(t) na origem com o comportamento de F(s) no infinito. 
hom Lt) = tem AFA) 
t —pO AA TPOO . 
(melo) 4” alontomda fofa? 
vii) Convolução 


Para calcular a convolução de duas funções, num intervalo de tempo [0,t], inverte-se no tempo 
uma delas, multiplicam-se uma pela outra e integra-se o produto, obtendo-se o integral de 
convolução. Num sistema linear, se calcularmos a convolução da resposta impulsional h(t) 


com uma entrada qualquer u(t), obtém-se a saída y(t) correspondente à entrada u(t), isto é, 


JC) f “hlt=) uto) dt 


O teorema da convolução da transformada de Laplace diz-nos que a transformada de Laplace 
da saída é o poduto da transformada de Laplace da resposta impulsional pela transformada de 


Laplace da entrada, se o sistema estiver relachado (com condições iniciais nulas), 


Tl) = Ha) lo) pt! 


u(t) EE y() 
fume dr Tamofentmia OS leds 
rantfent a tina fas a sala ). 
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Veremos que H(s) se chama a função de transferência do sistema. 


viii) deslocamento no domínio temporal 


Desesmeak wo daminio denfunal — muliflcneg 
Ama dmnio Copia 


LÍftt-a) = 2" F(4) 


3A.5 Inversão da transformada de Laplace: expansão em fracções parciais. 


Uma transformada de Laplace é em geral uma função racional, dada por 


Fl). Mio) | frmiço pacinal 
D(») Ma): prt nec dn A Ceefricintts Novo 


Dlp): peinims den a + ” 


Diz-se que ela é: 


Pe pr ap ap 

sa Ãa meÃT nar : vlo) « D/a 
E) bepréb prva, : dg Ee qem a 
Han) atdaRlpra, : Gram Nv(a) > gra D/6) 


se do do “ PR RA oe Gods Ela) *” — Gram d Dig), 
pi term ho manada 2 
fo dino nm 


Se F(s) é própria, pode decompor-se numa soma de fracções simples, cujas transformadas são 


fáceis de calcular (consultando uma tabela, por exemplo). 
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Flp)= Ao +A, flo) + A, Ela eee ha lo) 


A, É oonkamtã Pl, mula se Flo) e sAintamedo próbaãs 
nas vuda pr Flo) fr bprópria.ho=f(0): 


flo), ml) untado, da Lodela tugas Inamatr - 
mada não BO, PM), es fa lE. 


As, hr, Am pão constant pras. 


Agora 


pue)= MS) + Mill) Me falta + mbbt), too 
Pode constatar-se que 


£l») fuer então — Foo) = o ou umband 
Flo) Antomuda frsgpnam, ano pia =0 

En) begun F(00) £o0, cla 
£ (n) am própria —— +flo)- too 


Definição de pólos e zeros fe F(s) 


D (») 


N(s) e D(s) são polinómios co-primos (não têm raízes comuns) com coeficientes reais. 


[hetol. AO A epuaeso Dla)=0 
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Er * vous da Lpuato N/n) =» 


Como sabemos, uma equação de coeficientes reais pode ter raízes 
e reais simples 
e reais múltiplas 
e complexas aos pares conjugados (conjugadas) 

e por isso os pólos e zeros podem ter qualquer dessas formas. 


Quando F(s) é própria, sendo N(s) e D(s) co-primos, o 


ou Jordem] ou Flo): numero ds bolos 
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Exemplo: 


3 
lp) = O antena Á ; biprob res 
214 P/+ IE nt420p +8 


) fadirzecp Aus meras 
den = [4,%, 18, to, 8]; de roots (den) Matlab 
+Paja, 2,2,2) 


Ate PO HIS A 4 20 A 8 = (a e(p42) 


| ABRA rea du rel Dlidada 3 
Elo) = - got en3-2a tn B-24 pe res simple, 


aq) (n42? 
CL) dacomposcas sm fracçods Ra 
 Fln)= Als Ay Do sho e DE PL 


Gi = 


tdo sdão linmo pivtim ma Elela 


14) ninversas 
fM)= Ao ECA, 4 O a DO > + Agta pe: 


coma alta. Ao, As, a, As, Bag 2x Cala da ratiBlnis 


O cálculo dos resíduos A; depende do tipo de polo. 


Ao é fácil de calcular, bastando para isso substituir s por co e levantar a indeterminação. 
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A, =F(00)=2 


A1 corresponde a um polo simples (-1), e calcula-se por 


1,3 
— À = aero] = 2(c1) + 0) maloa), 3 


per (1+2)2 


O processo de cálculo facilita-se se fizermos a redução de F(s) 


443 2 
F(n) = En A é o 2A+ro (A tTAeIito +200+8 
20% -Jhp? Bop" Voa é E) 
-43p? 36h". 424 —16 


3,2 | 
Flo): 2 q AAT36N-M2A-46 | . 2 Filo) 
nº +20 18p2 424048 


E agora trabalharmos sobre Fi(s), decompondo-a em fracções simples 


Flnls 13/3602 420 +16 R Eco A As Pia 
(241) (042) At (+42) (042) (042 


na = 
ad 


No caso de polos múltiplos, têm que se introduzir termos com denominadores com expoentes 


desde 1 até à multiplicidade da raiz (neste caso 3) 


Agora, para cada um dos polos simples, faz-se 


Av 
Asp Aga(o4p)F(n) 
A=-py 
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No caso de um polo múltiplo, é um pouco mais complicado. 


bol, mula. multilcAsd n, 
da, Aron 


A+p (osp)* (a tp)” 
dz (A4b) PRO 


Apa E s [esto] 


Do] 1 fe pPFto] mi 


13134360" puza +16 


- L+2)3 
fc O a ps! ( 
AS . o 
24 +1E (021436 (-2) +uz(-2)+18. ,26 
A+ . 
Au q apto t e 
pai pa-2 
d. 1iaó stato 1] 
As É Ea 
á dn A+) pa 
2 +16) 
A 22 1442 (249) -(13p/+360 4424 
34) + 


= -26 
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dê (ot eoiree ud 


e portanto 


«Ut 2 < 
-á “e - ut “o 
s4)= 264)+32 —162 126 EL 


Caso de raízes complexas no denominador de F(s) 


2(pl= ta) E folina nato da tom crftuela Mas 
Do) pénis de A um itficinis, ataio 


Como o polinómio denominador tem sempre coeficientes reais, as suas raízes, se foram 


complexas, têm que aparecer aos pares conjugados, isto é, 


ii f = Tag My Aona 
. E 
então pº -(G- p! ty tam We 
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Exemplo (Chen, 166) 


fl). Do = cá 


(a40)(n2,2. 045) (244) (pe 1-2) )(p 40 +27) 


pot : 1505 


Flo] - A, pode + fa 
p4t (pat-2p) (n41+427) 


pra deimendlarar- Le pas A =4) 


cha dewmpoicgo mas + cmuenienii porcpu obra 
a tCabalha. ima ptrduro (Completos 


di L Edo +42 - 2Ajt041 Adao 


a ArZA +! tl =(4+1)S +22 


tm ; ado qu 
di a ay 


A-p)(o-b*) =(140 060 +a+yb) 


2 ,2 
= 4º4 2044 +aº +b = (144) +6 
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Fl») = a io iii = 
(040841) 42º] 


No Cola 
O pen it ds [8 
(s40)+ ui 


“at 
À cf vê e—s 
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A, Az A+ 43 

Ed 

A+1 (Aa) 42 
et 


É umveniuda cur go 


Luvvo Um dumemnada 
yradaçtia. hevham 


ALA 
RES Sá tha 
(ata +u? pa 
Fla) =A, pg LD pe ds as 
2,,* .,.2 
A+d (n+44)2+2 (n+)+2 


77 


Falta agora calcular Ai, Az, As: 


+ « 
+ Az 2 o2t + Age snat 


- rsllunde de dy (nrmidua de uma puto cast drepho,) 


Ag = (A+) fp) — 


A=-4 (244) 2º fran)” 
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- vealludk da ha 2 Rs 


Fln) = A EAR = Aa (n+1) + 283 
(241) [a 49 2427] p +43 (ato r4 Defesa 


: fine CU) + hz (249) (041) + 2A3 (049) 
(2+1) [++ 4] 


=» (244 445) (Ih) + Az (0242444) +2A3/0+44) =A 
(é + Ant + (-L +2A2 +283)A + (Era +283)= A 


Igualando os coeficientes de iguais potências de s de ambos os lados da equação: 


o hz = Ya 


N 


LA 
q the 
Lethc+2As = As = 1h 


E rhotihs = ÕÔ 


Invertendo, 


fit) > E pel Ea dt +) Pai 


obtém-se f(t). 


34.6. Aplicação da transformada de Laplace à resolução de equações diferenciais 


Seja por exemplo a equação diferencial de 2*ordem 
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aq +bq+Cy= u 


v Ê [ay+byecy] = 2[u] 
1) pra Lentes ds É 
a £T5]+b El] +e Elg)= É [4] 
cui) jul direna ds decimarso, 
ETg)= AP Liga yo) - lo) = 0º Yo) agto)-jfo) 


LT = 2 tfy]-gt) =4Yo)-gto) 
L[y)= Yo) 
Substituindo, 
sta Yo) - Aa Yo) - ato) + ab lo) - byho) eo) = (o) 
(a 4º + bp +c)Yln) = (an +)ylo)+ ago) + Ulo) 
Ya) = (0 n+b)96) +ayto) (Uta) 


antebara 2/24 bar 
la ——————s DE A 
sonda ouvida às Condia parda eunda 
Aura , Var à imlanda, Yes 


Veremos o significado das duas parcelas de Y(s). 
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Tabela das Propriedades da Transformada de Laplace 


Property Time function Laplace transform 
Linearity ou fi(t) + oo fo(t) oyFi(s) + ooF o(s) 
s-shifting etF(t) Fí(s+a) 
Time shifting f4-Tag(t-T), forTz0 esTF(s) 
=d 

s-differentiation (1) ae (s) 
Time differentiation Es f(t) sF(s)-— f (0-) 

dt 


d? 
Ex) fc) 
Time integration 


Time scaling 


Convolution 


Final value 


Initial value 


(DADC/FCTUC/2022 


t 
Jo fioar 
flat), fora >0 


y h(t- Tu (vdr 
f(eo) = tim sF(s) 


f(0+4) = lim sF'(s) 


SF(s)-sf(0)- f(0) 


Fís) 
s 
Ares) 
a 


HO(s)U(s) 


(if sF(s) has no open RHP 
poles) 


(1f F(s) is strictly proper) 
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Tabela de Transformadas de Laplace 


f(1),120 F(s) 
ot) l 
| A 
or q(t 
q(t) E 
p o, 
s? 
E ] 
t” (n: positive integer) n: 
gr+ 
e“ (a: real or complex) 1 
s+a 
tera l 
(s+a) 
meat n! 
(s PP à agia 
sin Qot Do 
s2 + 05 
cos Dot PER re 
s2 + 
sin Oo! — SS 
(2 +09? 
2 2 
É Cos Mot Ss — O 
(57 + 05) 
e sin wot RE. + 
(s+a) + Ds 
e“ cos wot — Sta 


(s +a) + q 
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4.1. Espaço de estados, variáveis de estado, equações de estado 


Contrariamente à função de transferência, definida no domínio complexo, a representação no 


espaço de estados desenvolve-se no domínio temporal. 
Exemplo 1 (Bruce): 


Imagine-se um músculo de rã fixo por uma extremidade e no qual se suspende um corpo de 
massa M. O músculo, devido à sua elasticidade, alonga-se. Pode-se desenhar um equivalente 
mecânico de translação com uma mola K (elasticidade), atrito viscoso R, massa M, como já 


fizemos no Cap. 3. 


A força aplicada à massa é a da gravidade e igual ao seu peso, F=Mg. Ou, se se desejar, pode- 


se-lhe aplicar uma força u (de distensão ou de compressão). 


O atrito e a mola reagem com as forças 


atrito: F, -BÉ By 


t muscle 


mola: F.=Ky À B 
A força resultante é a soma das três e produz o deslocamento x il 
inte 
Mg-Ky-By=My > My+By+Ky=Mg a 
Fig. 4.1. Músculo da rã suspenso 
é Ba R por uma extremidade (Bruce). 
ás M” Ea M” =9 


Uma equação diferencial de segunda ordem, linear. 


Pode-se reduzir a duas de 1º ordem, por substituição de variáveis, 


E 


X=y W=y 


Agora A . - 
=P M=Y=X% 


N e e . . B K 1 
MED) st 
B K 1 


X é cl 
M M M 


ou seja x,= 
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Portanto 
X =*% 
E B K 1 
re A x +—u 
M M M 
O) 1 1 


Ou ainda, mais simplesmente, 
x = Ax+ Bu 
Exemplo 2 


No modelo de inspiração-expiração pulmonar, definindo a variável de estado x1 como o valor 
da pressão VL, a entrada u como P(t) e a saída do sistema como o que queremos observar- o 


volume pulmonar, teremos 


Vit) av (t) = bP(t) 


x 
H> 
< 


uê P(t) 


x, =—ax, +bu 


uma equação diferencial de 1º ordem, linear, de coeficientes constante (por isso invariante), 
completa por haver uma entrada externa, u(t). Deve ser escolhida como estado uma variável 
que seja capaz de exprimir a memória do sistema, em geral, como vimos no Capítulo 2, 
associada à quantidade de energia ou de massa armazenada. Neste caso essa variável é o 


volume pulmonar (associado à quantidade de ar). 
Exemplo 3 


No modelo de Lotka-Volterra que também vimos no Capítulo 2, 
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dx(t) = ax(t) — bx(t) y(t) 
= = —cy(t) + px(t) y(t) 


definindo as variáveis de estado x1 e x2 como o número de elementos de cada população 


(acumulação de massa), substituindo, simplificando a notação, 


x, é x (presas) x, 2y (predadores) 
1 
dx,(t E 
eua =ax(D—bx (Dx,(t) x =x, — bx,x, 
dx,(t) . 
o = —ex,(t)+ px (Dx, (t) X =X + PX,X, 


obtêm-se duas equações de estado, invariantes (porque os seus coeficientes são constantes 
com o tempo). Não há entrada externa, neste caso, e por isso as equações são homogéneas. 


Com há produtos de variáveis de estado, as equações são não-lineares. 


Se nos interessar observar a evolução das duas populações, diz-se que elas são as saídas do 


sistema (o que se vê do lado de fora), e assim 


saída 1:y, =X, 


saída 2: y, =X, 


Exemplo 4. Metabolismo da glucose 


O metabolismo da glucose (ou glicose) no sangue pode ser aproximado pelo modelo 


simplificado seguinte (Ackerman et al. (1969) (em Linkens)) 
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d 
ao Tm —m,h+J(t) 
Se=-mh+mg + K(t) 


g: desvio do nível da glucose do seu valor recomendado 
h: desvio do nível da insulina do seu valor recomendado 
J: taxa experimental de infusão de glucose 
K: taxa experimental de infusão de insulina 


m,m,,m,,m,, constantes características de cada indivíduo. 


Podemos fazer o mesmo desenvolvimento para todos os exemplos que vimos. Por isso 
podemos generalizar, dizendo que a estrutura geral de um modelo de um sistema contínuo, 
independentemente da sua natureza, é constituída por um conjunto de equações diferenciais de 
primeira ordem, com m entradas externas. Usando uma notação genérica, vem, para n 


variáveis de estado, 
dx. j 
eo fi CO, ..., x (O, u(l, ..., u (bt, t) E Laval 


em que f; são funções contínuas dos seus argumentos. As condições iniciais, necessárias para 


definir o estado inicial do sistema (a sua memória), serão, 
Xxi(t)=xo i=1,...,n 
Teremos também um conjunto de r equações de saída (r saídas) 
x (O = 9, (lt, ... x,.(0, u(t), ... uti), t) i=1,..sr 


em que as variáveis de saída medidas se exprimem em funções das n variáveis de estado do 
sistema e das m variáveis de entrada. Não há nenhuma relação geral obrigatória entre as 
dimensões dos vetores de estado, de entrada e de saída, mas na maior parte dos casos n >> 


máx (m, r). 


Se f; ou g; dependem explicitamente de t, para algum i, como por exemplo em 


x=-2xe'+3u 
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diz-se um sistema variante. Quando tal não acontece, chama-se sistema invariante. 


As funções fi têm como argumentos apenas xi, quando não existe qualquer entrada externa 


u(t). Chamam-se neste caso sistemas autónomos (dependem apenas de si próprios). 


Nas funções gi é pouco frequente existirem entradas como argumentos: as entradas 
influenciam as saídas através dos estados, e não diretamente. As saídas são por isso em geral 
funções apenas dos estados. As saídas são, no fundo, o que queremos observar, medir, 


calcular, no sistema. 


Em notação matricial poderemos escrever 


X XY: u, Xio 


X Y, X Xno 


(vetor de estado) (vetor de saída) (vetor de entrada) (vetor de estado inicial) 


e os vetores de funções 


If 9, 


fogu)= |. gu) = 


fi 9, |] 


O modelo no espaço de estados é então definido por uma equação de estado (com condições 


iniciais) e uma equação de saída. 


f(,u), x(to) = x0 
9(x,u) 


II 


Rd 


Se as funções fi e gi são lineares e invariantes, teremos 
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OGU) = 4,%X + Bd Fe ao AD do Dia os TS Las 
JAR) = CD ta Pd e Do DE dat 


em que aijj, bij, Cij e dj são coeficientes constantes. 


Ou seja, a equação de estado toma a forma matricial 


x= Ax + Bu x(ty) =X, 
y = Cx + Du 


que se pode representar pelo diagrama de blocos da Fig. 4.2. 


Condições 


iniciais 
I (memória) 


Integrador 


Figura 4.2.Diagrama de blocos da representação no espaço de estados. Normalmente D é nula. 


As matrizes chamam-se 
A=nxn matriz de estado 
B=nx m matriz de controlo, ou de entrada ou de comando. 
C=rxn matriz de saída (ou de observação) 


D=rx m matriz de avanço ("feedforward") 


u = vetor de controlo, ou de entrada, de dimensão m 


y = vetor de saída de dimensão r 
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No caso SISO (Single Iutput, Single Output) B é um vetor coluna, que se pode anotar por b, e 


C é um vetor linha, anotada por ct, D é um escalar d. Assim vem 


x= Ax + bu x(ty) = x, 


y = c'x + du 
mas neste caso SISO também se pode utilizar a notação mais geral (A, B, C, D). 


Se as matrizes A, B, €, ou D contêm elementos que variam com o tempo, diz-se um sistema 
variante (neste caso algum dos coeficientes ajj, bij, Cij ou di; variam com o tempo) e anotam-se 


por A(t), B(t), C(t), e D(t). 


A matriz D exprime a influência instantânea das entradas sobre as saídas, ou seja, aquela 
influência que atua imediatamente e antes da devida à variação dos estados. Por isso se chama 
matriz de avanço. Na generalidade dos sistemas físicos, providos de inércia, esta influência 
instantânea é nula e, portanto, a matriz D é a matriz nula. Lembremo-nos, por exemplo, do 
efeito da posição do acelerador de um automóvel sobre a velocidade por ele atingida - o tempo 
de reação é sempre finito não nulo (a não ser que o automóvel possua um número infinito de 


cavalos...) não existindo por isso ação de avanço. 


Definição 4.1 


Estado. O estado de um sistema é uma estrutura matemática constituída por um conjunto de n 
variáveis x;(t), x>(t), ..., xn(t), chamadas variáveis de estado, em que os valores iniciais x;(t0) 
desse conjunto (i=1, ...,n) e as entradas ujct) (=1,...,m) para t > ty são suficientes para 
descrever univocamente a resposta futura do sistema. Existe um conjunto mínimo das 


variáveis de estado requerido para representar o sistema. 


O instante inicial ty toma-se usualmente como a origem do tempo, ty-0, mas pode ter um 


valor qualquer. 


Definição 4.2 


Vetor de estado (ou simplesmente estado). O conjunto das variáveis de estado x;(t), i= 1,...,n 


representa os elementos ou componentes do vetor de estado x(t) de dimensão n, 
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xi(t) X 
X(t) pa 
também chamado simplesmente estado. 
Quando todas as entradas ur, j=1,...,m, são especificadas para t>ty, o vetor de estado 


resultante descreve univocamente o comportamento do sistema para t > ty, dados os valores 


de x(to). 


Definição 4.3 


Espaço de estados. O espaço de estados é o espaço XN de dimensão n, subespaço de IRN em 
cujos eixos de coordenadas se definem os componentes do vetor de estado. O estado x(t), em 


qualquer instante t, é um ponto de IRN. 


Definição 4. 4 


Trajetória de estado. Quando o estado varia com o tempo, chama-se trajetória de estado ao 
percurso (curva) por ele produzido no espaço de estados. Alguns autores chamam plano de 
fase ao espaço de estados bidimensional e trajetória de fase à trajetória de estado nesse caso 


(mais tarde veremos porquê). 


Definição 4.5 


Equação de estado. A equação de estado de um sistema é um conjunto de n equações 


diferenciais de primeira ordem, em que n é o número de variáveis de estado. 
4.2. Escolha das variáveis de estado 


A primeira etapa na aplicação destas definições a um sistema qualquer consiste na escolha das 
variáveis (do processo) que devem ser variáveis de estado. Deve notar-se que não existe um 
modo único de fazer essa escolha. Na prática existem várias escolhas possíveis, e teoricamente 
uma infinidade. Algumas escolhas têm interesse particular. Quando as variáveis de estado 


representam grandezas fisicamente mensuráveis e observáveis chamam-se variáveis físicas. 


A escolha das variáveis de estado físicas (tal como se fez nos exemplos anteriores) baseia-se 


em geral nos elementos do sistema armazenadores de energia ou de massa. A tabela 4.1, já 
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estudada no Cap. 1, indica alguns elementos armazenadores de energia (mais comuns) 


existentes em sistemas físicos e as correspondentes equações energéticas. 


A variável física da equação de energia de cada elemento armazenador pode ser escolhida 


como variável de estado do sistema. Em alguns sistemas pode ser necessário escolher outras 


variáveis para além das dos elementos armazenadores de energia. 


Os elementos armazenadores de energia e de massa sintetizam toda a história passada do 


processo. De facto, sendo todos os fenómenos físicos transformações de energia e de massa, o 


estado atual da massa e de energia de um dado processo é o resultado de todo o que ele 


passou. Daí a pertinência desta forma de escolha das variáveis de estado. 


Tabela 4.1. Elementos armazenadores de energia 


Sistemas Elementos Energia Variável 
2 
Condensador C “e tensão v 
Elétrico 3 
Indutância L i 
É) velocidade de 
w 
Massa M M 
2 translação w 
. 2 loci taçã 
Mecânico Momento de inércia J J o MELPsAERS duna ação 
w 
2: 
Mola K K a deslocamento x 
Gás 
Compressibilidade de v Pp E 
E ar pressão 
| um fluido V/Ks, KB KB 2 
Fluídico 
constante de Boltzman 
qa ne C si altura (nível) h 
o u v 
Capacidade de fluido C 2 
g 
Térmico Capacidade Térmica C A temperatura O 
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4.3. Resolução da equação de estado linear 


dado 
Xlo)=% , tfade aricar 


ult),tsz0 , enfiada afluns 
pode - se Laluda X(t),a tasnjedais da edad, 


netolvendr à efrncgo du aÃeas ; ve 2, X00) 
2 à soluncça da equaçp de estadas. 


Para a variável de estado x; teremos, 


o ZH z AX) +.+ A, UE) + GU EI + est Plim tt) 
Lau] = É [yu | 


2 av, L[xt))+ ot Ante) + bes dl [uitt)) 
Aplicando o teorema da derivação da Transformada de Laplace 


A Xilal-Xito) = Acy XlA) ++ Ag, Mala) + be, Uno) +: + ba, Uma) 


Fazendo o mesmo para todas as variáveis de estado, e escrevendo em forma matricial, vem, 


para a equação de estado 


— X= Ax+Bu Z [x). A 2[x)+8F[u) 
A Xla)-Xlo) = A X(n)+B Ulo) 


e para a equação de saída, 
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— ye Cx+Du L ly] = cflx]+DÊ[u) 


Yo)= Cxto)+D Ula) 


Aplicando a Transformada de Laplace à equação de estado (propriedade da derivação), 


A Xlo) = eloa AX) + BUG) 
A Xla)- AXin) = Xo + BULA) 


ALXu)-AX) = Ze + BUG) 


[aT—A] X 4) a X + Bula) 


| X (a) = [aa x, +[oI-A] BU(A) 


A oqecniininadcasai reencontra —s 
Xi tfrajedeirinas Xsç: Cumpedtrin ou 


obtém-se a solução no domínio complexo. 


Substituindo X(s) na equação de saída, 


Ylola CXlol+ DU(S) 


Yhl = CLnEA'y + chc-s] Bus) + DU(») 


Tlsla [aI- AJ o +[c (aI=]"'B +D] Uln) 
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reencontramos aqui as duas partes da saída do Cap.3: a devida às condições iniciais e a devida 
à entrada. 


Resposta com condições iniciais nulas: 


Yu)= [cCoT-a!'B +DJ]U(») 


for tiro Bos, sebemo Ga 


Lho)= Gl) VUlo) 


-2 -3 1 iiibi il, sa 
qa Co -ilx tsdudo  Y(£) 
Calculando, 
pra = ofosfoE 9] o dal 


u 


ara E R 3 | 


p(n+43)+2 4 A 


X(n)= [al- AT 'Xo + EA b 4 
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X(n)= [al A! xo +PI-A b 


A+3 4 ” 
- = 1 2 o / 
o (a+2(0+1) (p40)(041) + [ara] dg 
cas a ps Re o) 
(n40)041) (0 420)(041) 


29" 45p+ 4 
Klnl= p(n42) (241) 
-A-3 
(o+0) (241) 
2,8p+) EP é 
Yo) = CL -1) iss = po nar 
p (042) (nt1) 


-2t ar 
ylt) = 05-48 2 +42 , tzo 


À furia au Tranofe-incia £ dada por 
G(n) = eT [nT-4] b +d 


blala [4 -1] 443 Í 
atom) faretoel) 


Z A 
(A+) (a +1) (14)(9+41) 


1 


u 


(1 -<] t 
(142) (241 
A 
(Av) (443) 


41-A 
(142) (241) 
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Aplicando uma entrada em degrau unitário 1/s, 


flo) POE 


a (041) (04%) 


' SE 
JH). gs + sc 28 


obtém-se uma saída diferente da anterior. Tal seria de esperar, dado que G(s) produz a saída 


para condições iniciais nulas e no caso anterior tinhamos condições iniciais iguais a [2 -1]". 
4.4. Solução da equação de estado no domínio temporal 
Solução da equação de estado homogénea 


Consideremos a equação de estado de um sistema linear, variante 


x= Ax+Bu. 


A equação homogénea obtém-se fazendo u=0, de que resulta um sistema dito não 


forçado (também chamado autónomo), 
x(t) = Ax(t) 


x(0)=] l Ax(t)dt + x, 


Mas o integral é recorrente em x(t), 
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x(t) = I, Ax(t)dt+ x, = I, A( K Ax(t)dt + x, )dt + x, = 


= [AC], A(Í, Ax(t)dt + x,)dt+ xo)dt + x, am 
=X + [Axe + | [ AMxçde +[] [max de e 
0 00 000 


=X + Ax, [des Aºx, | [de +a'x, [far se 
0 00 000 


f? f tt 
=x, + Axgt+ A?x, TRE ter Ã ão nte 


2 3 k 
ps Als pls 4 Atl Tx 
2! E k! 


Se compararmos esta última expressão com 


eat= 1+ Mt + (a ++ (E a 


poderemos escrever, por analogia, 


x(t) = e“ x(0) 


e por isso e* chama-se exponencial de matriz. 

Note-se que a exponencial da matriz A não é composta pelas exponenciais dos elementos de A ! 

A exponencial de matriz chama-se matriz de transição de estado, D(t), podendo escrever-se 
x(t) = D(t)x(0) 


Conhecida a condição inicial, para calcular o estado em qualquer instante futuro basta 
multiplicar o estado inicial pela matriz de transição de estado nesse instante; isto é, esta matriz 


transita o estado inicial para o instante t, e daí o seu nome. 
Como calcular e! ? 


Computacionalmente poderemos calcular parcela a parcela, 
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t? t t' 
E ET rAprAÃs FÃ ds dAPS= E 
21 3! k! 


até que a parcela seguinte seja praticamente nula (se a série convergir ...). 


Mas seria interessante procurar outra solução, menos trabalhosa e (ouro sobre azul!) que pelo 


caminho nos permitisse também extrair algumas propriedades do sistema. 


Essa solução existe, e baseia-se na estrutura própria da matriz. Como se sabe a estrutura 


própria da matriz é composta por 
- valores próprios, 
- vetores próprios à esquerda, 
- vetores próprios à direita. 


Considere-se novamente a equação homogénea 
x(t) = Ax(t) 
Uma solução possível para este sistema é aquela em que xe x têm a mesma direção no 


espaço de estados, diferindo apenas na sua amplitude por um fator de proporcionalidade À. 


A solução será então de forma 


E =Ax 
Substituindo acima, obtém-se 
Ax=Ax <=>" [AI-A]x=0 


Se um valor de x * O satisfaz esta equação, é uma solução não trivial, e teremos 


necessariamente 
Q0)=det [A I-A]-0S |[AI-A|=0 


Sendo A de ordem n, Q(A) é o polinómio característico de A e a equação polinomial 


resultante de chama-se a equação característica de A e terá a forma 
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QM)=1"+a A+ .+al+a=0 


As raízes 4; da equação característica chamam-se os valores próprios da matriz 4. Uma raiz 
pode ser simples (distinta), ou repetida com uma multiplicidade p. Uma raiz pode ainda ser 
real ou complexa, aparecendo neste último caso sempre em pares conjugados (dado que os 
coeficientes da equação característica são reais). Quando as raízes são distintas, o polinómio 


Q( À) pode escrever-se 


U4)=(4-11) (41-12)... (4 - An) 


O produto dos valores próprios de uma matriz A, sejam ou não distintos, é igual ao seu 
determinante. A sua soma é igual ao traço da matriz (soma dos elementos da diagonal 


principal). 
Aos vetores x não nulos tais que 
Ax= AX 
chamam-se vetores próprios à direita de 4, porque multiplicam por A à direita. 


Se os valores próprios são distintos, a cada um corresponde um vetor próprio, e os n vetores 


próprios assim obtidos são linearmente independentes. 


Se existe um valor próprio de multiplicidade p, o número de vetores próprios a ele associados 


linearmente independentes varia entre 1 e p, dependendo das propriedades da matriz A. 
Note-se que se x é vetor próprio, ox, o e IR, é-o também. Prova-se pela definição. 


Chamam-se vetores próprios à esquerda (ou recíprocos) da matriz A os vetores wi, 


(associados aos valores próprios 2;, 1= 1, ..., n) tais que (multiplicando à esquerda por 4), 
wlA=Awl o wI[A;I-A]]-0 
Se wj é um vetor próprio recíproco, também o é awi, a e IR. 


Pode provar-se que a solução da equação de estado homogénea, usando a estrutura própria é 


dada por, no caso em que os valores próprios são todos distintos, 
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o At 
= j 
xX(t) = |D &.v,w; |x(0) 
ja 


=etvw x0) ret. vw x(0). +... rev, wi x(0) 


Este resultado é cheio de significado: ele diz-nos que a trajetória temporal do estado, para uma 
dada condição inicial, é uma soma ponderada de exponenciais dos valores próprios da matriz 
A (os termos e). São os vetores próprios que estabelecem os coeficientes de ponderação. 
Quer dizer que a estrutura própria da matriz é como que o código genético da dinâmica do 


sistema (autónomo). 


O que acontece se um dos valores próprios da matriz de estado é positivo ? A sua exponencial 
tende para infinito com t e portanto a parcela respetiva também ! Logo a soma é infinita ! Isto 
é, o estado tende para infinito (pelo menos num dos seus componentes) e por isso o sistema é 


instável. 


Reencontramos aqui a questão da estabilidade às condições iniciais, tal como no Capítulo 3. 
Os valores próprios de A são também chamados os modos do sistema. A esta solução da 


equação de estado chama-se solução modal por explicitar os modos. 
Qual a relação entre estes modos e os que estudámos no Capítulo 3 ? 


Pode-se demonstrar que são exatamente os mesmos! Isto é, o polinómio característico do Cap. 
3 é o mesmo que o polinômio característico da matriz de estado. As raízes do polinómio 


característico são os valores próprios da matriz de estado. 
Tudo o que se disse no Capítulo 3 sobre os modos se pode repetir aqui. 
Solução da equação de estado completa 


Agora temos 


x = Ax+ Bu 
e encontrar a solução não é tão simples. 


Para isso considere-se a propriedade da derivação de matrizes variantes 
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MONO ireoneo+ Meo Meo 


em que a derivada de uma matriz é a matriz constituída pelas derivadas dos elementos da 


matriz inicial, M = [mi 


[e *x60] e Xt)-Aet xt) =e [xo - Axo | 


mas como X(t) = Ax(t) + Bu(t) 


S[e*x00] = e“ Bu(t) 


integrando ambos os membros da equação 


e x(t) = | e““Bu(r)dr + k 


0 


Sendo k uma (matriz) constante de integração. Se substituirmos t = O na expressão 
0 
e“ x(0) = ] e“Bu(r)dr+K & H0=K 
0 


Multiplicando ambos os lados por e4t e agrupando os termos, obteremos, 
t 


etex(t)=e! | e“Bu(r)dr + e x(0) 
0 


IX(b) = e“'x(0) + | e“DBu(r)dr 
0 


XD=X,+X. 
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Pode constatar-se que x(t) é igual à soma de duas contribuições, 


x(t) = xzi (b) + xzs (b) 


uma devida às condições iniciais (zi-zero input) e outra devida à entrada (zs-zero inicial state). 


4.5 Sistemas não-lineares, estados de equilíbrio e linearização 


Muitos sistemas naturais, em particular os biológicos e fisiológicos, são de facto descritos por 
equações diferenciais não lineares. Coloca-se então a questão de aproximar as equações não 
lineares por equações lineares, se se pretender aplicar a teoria dos sistemas lineares. Chama-se 


linearização a essa operação e faz-se em torno dos estados de equilíbrio. 
4.4.1 Estado de equilíbrio 


Em muitas aplicações práticas as condições ideais de operação são caracterizadas por 
entradas, estados e saídas constantes no tempo. Suponhamos que o sistema é descrito pelo 


modelo não linear 


xe f(x, u) 


A representação matemática de uma operação de linearização faz-se passando pela noção de 


ponto de equilíbrio ou estado permanente. 


Definição 4.6 


Estado de equilíbrio. O vetor xs (n x 1) é um ponto de equilíbrio, estado de equilíbrio, ou 


estado permanente, correspondendo à entrada constante u(t) = us, se e só se, 
f(x, b) U,) = 0 


A correspondente saída de equilíbrio ou saída de regime permanente é o vetor yg (r x 1) tal 


que 


y, = 9(x,,u,) 
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Exemplo 1 


Seja o sistema definido pelas equações de estado e de saída 


X=M+*% -2u 


X2 =X -X 
Pá 
y =x +1/2x 


Para calcular o estado de equilíbrio, aquele em que o estado não varia, faz-se 


ou, equivalentemente, iguala-se a zero a parte direita 
xy tx2-2u=0 

x1-x29=0 
y=x9 + 1/2x2, 


e resolve-se em ordem a x, e x,. 
x1 tx9-2u=0 2x1 = 2Uç X1=uUs Xg = Us 
é o ares 


Ys= uç(l + 1/2uç) 


Note-se que a cada valor de entrada u = us, corresponde apenas um valor de estado de 


equilíbrio; por exemplo para us=2, x15s=2, X252, Ys6 
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Exemplo 2 


Considere-se o caso seguinte e o cálculo dos seus estados de equilíbrio, 


X =X +X% -2U 


x =X (X -X) 


x1=0 ORE A o x =0>x=2u 
> 
e 


x (x, - x,) =0 —4=0>x=,>X=4U 


x2=0 
verifica-se que tem dois estados de equilíbrio para cada valor us 


l 2 


Dado um us, qual o estado de equilíbrio que se obtém, x1s ou x25? 


Depende das condições iniciais x,(0) e x,(0). Podemos simular no Matlab o que acontece. O 
programa PPLANES (https://www .mathworks.com/matlabcentral/fileexchange/61636-pplane) 
faz o seguinte: dado o valor de u, para cada valor de x1(0) e x2(0) integra as duas equações de 
estado, calculando assim as trajetórias xi(t) e x2(t). Fazendo isso para um grande número de 
valores das condições iniciais, obtém-se a figura, no caso de u=1, chamada curvas de fase 


(Phase-Plane, PPlane). O programa aceita duas variáveis de estado, x e y. 
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Figura 4.3. Exemplo 2, Curvas de fase. As setas indicam o sentido de progressão da trajetória. 


A notação é x = x1, y=X2 ; note-se que aqui y não é necessariamente a saída do sistema. 


Existe também uma versão Java para o pplane em https:/Awww.cs .unm.edu/-joel/dfield/ 


Exemplo 3 

Seja 
X =X +UX x +ux,=0 x(u-D=0 => x=0 ou u-1=0 
X =X +% x +4=0 x=—* 


Agora temos duas situações, conforme o valor de u, 


0 
1) uz1,x=x/-0, = 
0 


l 
» hE|R 


u=1, x2qualquere xj=-x> Xg=A | 
-] 


Assim para u=1 existe uma infinidade de pontos de equilíbrio. As condições iniciais 


determinam qual o obtido. 
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u= u=1 


Figura 4.4. Exemplo , curvas de fase. Para u=2 a origem é o único ponto singular (instável). 
Para u=1 todos os pontos situados na reta a sombreado são singularidades instáveis. 


Exemplo 4 
Para o caso seguinte 

mM =X -x X =» 

x 14 (=x) +; I+0u +x) =0> x = hj (imaginário) 
Se o sistema for real a solução obtida corresponde à inexistência de estados de equilíbrio. 
Exemplo 5 
No caso geral linear, 


x = Ax+ Bu 

x=0= Ax+Bu=0= Ax, =-Bu, 
e dois casos se podem dar: 

a) A é não singular, solução única x = - Al Bug 


b) A é singular: não existe solução ou existe um número infinito de soluções. 


Por exemplo, no caso seguinte, 
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o A -1 0 
X= X+ u 
| 1 + 


de estado de equilíbrio, para um dado us será 


di OR PR : RES 10 
X Uç— ou Xe =— u 
ERA a ese quo QUA 
À a E 
se Us = A = 
0 1 — -1 


não existe solução possível. 
1/2 1-17|X% 1/2 
Se Us = E = 
1 1-1 1/2 


existe um número infinito de estados de equilíbrio definidas por 


x1s - X2s= 1/2. 


Note-se que quando existe uma solução xç, a singularidade de A quer dizer que tem um (pelo 
menos) valor próprio nulo e portanto um vetor não nulo satisfazendo Aw = 0. Então A (xg + 
Aw) = Axg + AAw = Axç = - Bus. Quer dizer que e A e |R,xg + Aw é também um regime 


permanente e portanto o número de estados de equilíbrio é infinito. 
4. 4.2. Métodos de linearização: aproximação pela Série de Taylor 


O processo de linearização (de um modelo) de equações não lineares na vizinhança de um 
estado de equilíbrio definido pelo tripleto xs, us, Yg» consiste na aproximação dos vetores de 


funções fe g por funções lineares, na vizinhança desse ponto de equilíbrio, 
Linearização pela Série de Taylor 


A linearização analítica baseia-se na aproximação das funções fj e gi por uma série de Taylor 
(na vizinhança dos pontos de equilíbrio), sendo Axi e Aui os desvios em relação ao estado de 


equilíbrio. 
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f(x + AX Xy, + AX... + x + AXU, Pe + Au e 


f (ça Maas + Xnso Uiço Usgo ee ds ADA qu A, 


7 OX, = O 


+ termos de ordem superior 


Para as funções de saída, de igual modo 


TX + AXys Ko, + AX. + Xas + AX, Uy, + AXjsU,, + AU,,...,5U + AU, ) = 


9 9 
Ei(Xiso Nosso + + Xnso Uso Uoçoee das) +D) ROSA +5 A 


1 OX, xa OU, 
+ termos de ordem superior 
Desprezando os termos de ordem superior e considerando, 
of; of | of, oh, 
Ox, Ox, Ou, OU, 
Di Ra = Fy BS | isa =P 
seda Si cálculada em 
E df Xs, Us 
Ox, Ga Rage ôu, ou, | 
ER E 
Ox, Ox, ou, ou, 
C=| = Gy D=| .. =Guí 
og og calculada em og og calculada em 
p ... —P Xs, Us —P ... DP Xs, Us 
Ox, Ox, ou, ou, 


A notação FX simplifica a compreensão da matriz: os índices de x evoluem por linha, e uma 


linha é a transposta de um vetor (neste caso x). O mesmo se aplica para as outras. 


As matrizes das primeiras derivadas das funções em ordem a cada um dos seus argumentos 


chamam-se Jacobianos. 


Substituindo na série de Taylor truncada na 1º ordem obtém-se um modelo linear 


Xs = f(x,,u,) y,=9(x,u,) 
(xs+ Ax) = f(x, + Ax,u, + Au) Ys + Ay = g(x, + Ax,u, + Au) 
x+Ax= f(x,u,)+ AMX + BAU y.+Ay = g(X,,u,) + CAx + DAu 
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Como xs=f(x,u) e Y=g(x,,u,) da expressão acima resulta a equação de estado e de 


saída na forma padrão 


Ns SAE BA 
Ay = CAx + DAu 


cuja validade só pode ser afirmada quando são desprezáveis os termos de ordem superior a 1 


da série de Taylor, 1. e., quando são pequenos os desvios do regime permanente. Para 


simplificar a notação, poderemos representar o sistema linearizado eliminando os A, e tendo 


sempre presente que as variáveis (de estado, de entrada, de saída) representam desvios do 


regime permanente. As matrizes A, B, €, e D são os Jacobianos de fe g (em ordem ax e au). 


Exemplo 5.Pode-se verificar por substituição que o sistema dado por 


x =-(x,-1)+x, +x, -2u 


E 2 2 FAR 
ex (el) RX uu, 


X3 =X, +X, -X, -U, 


ySx(ltx)+u, 


Yo =X, TX -U, 


l 
para u, = É tem um ponto de equilíbrio em 


Linearizando em torno desse ponto obtém-se 


Será esse estado de equilíbrio estável ? 


xg=|1 


Calculando os valores próprios de A, no Matlab 
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>> A=[011;310;11-1] >> v=eig(A) 

A= v= 
0 1/1 2.6274 
3 10 -1.3137 + 0.4211i 
o -1.3137 - 0.4211i 


conclui-se que é instável, dado que um dos valores próprios tem parte real positiva. O que 


acontece quando uma trajetória se aproxima dele ? Veremos no Cap.5. 


Conclusão 


Neste capítulo estudámos alguns elementos sobre o espaço de estados, uma representação 
temporal de um sistema. Aplica-se de igual modo aos sistemas lineares e não lineares, 


variantes ou invariantes. 


No caso linear obtém-se uma representação matricial. As propriedades dinâmicas do sistema 
são dependentes dos valores próprios da matriz de estado, tal como são dependentes dos polos 


da função de transferência na representação no domínio complexo. 


Os sistemas não-lineares podem ter zero, um ou vários estados de equilíbrio para a mesma 


entrada. Alcançam um ou outro conforme as condições iniciais. 


Aproximando as funções de estado e de saída pela série de Taylor, desprezando os termos de 


ordem superior à primeira, obtém-se um sistema linear. 
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5.1. Pontos singulares nos sistemas não-lineares 


Os pontos singulares num sistema são aqueles em que se anulam as derivadas das variáveis de 
estado. Chamámo-los pontos de equilíbrio no Cap. 4 e por vezes também se chamam regime 
permanente; quanto o ponto singular é estável, estas denominações fazem sentido. Caso 
contrário é mais sugestivo chamar-lhe ponto singular, dado que as trajetórias não acabam 


neles, mas são por eles repelidas. 
5.1.1 Análise no caso linear 


Um sistema linear, salvo para a situação particular em que det(A)=0, tem uma só 
singularidade. Se a entrada for nula, ela será a origem. Vamos considerar este caso, sem perda 
de generalidade (como sabemos no caso linear a entrada desloca a singularidade no espaço, 
mas não altera a sua natureza). A estabilidade do sistema é dada pelos valores próprios da 


matriz de estado. 
Há vários tipos de singularidades, umas instáveis, outras estáveis. 
Vejamos um exemplo de cada tipo. 


1) Nó estável (poço-sink, attractor): quando os valores próprios de A são negativos e 


distintos. A trajetória tende para o nó sem oscilações. 


“12 0 1 São º 
= X + u I ] | I l 

(0) —3 1 A SAS Se Cd a GA PS À > + 
R | | | 
I » I I ) | 

05h 4 — = — += 4 ——— > — 4 [> — — — — +- + 
) I I I I 
I I I I | 
! | WU, | | 

» 0-1 — 1555 == + [re cem TE dps] 
| ) I I I 
I I I q I I 
; E I ET I Ne J I 

valores próprios: -2, -3 Da da Ra SD 
z e ; 

bi Vo PEA, ol cd, LIAN NS NNM 4] 
ada I I h A] 

ii) Foco estável (poço-sink, atractor): quando as raízes são complexas conjugadas 


com partes reais negativas. A trajetória tende para o foco oscilando. 
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valores próprios SALA TADÃS q E à VA Nie +] 


-0.5000 + 1.6583i : É, 


-0.5000 - 1.65831 


Hi) Nó instável (fonte-source, repelling): ambos os valores próprios são reais 


positivos. A trajetória é repelida, sem oscilações 


. 2 0 1 Visáy ai di 
Re x+) ju ENTAO. 
0 3 | RRUR ANA ATi Aa ”| 


valores próprios: 2, 3 belo cl a e q 


iv) Foco instável (fonte-source, repelling): quando as raízes são complexas 


conjugadas com parte real positiva. A trajetória é repelida com oscilações. 


e 01 0 A 
X= X+ u Ed RO LAPA sLN pod 
a dp dd a 


valores próprios: 


0.5000 + 1.65831 1a 


0.5000 - 1.65831 
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v) Centro: quando os valores próprios são ambos imaginários. A trajetória move-se 


em torno do centro descrevendo uma curva fechada. O raio da curva depende das 


condições iniciais. 


s ope)! 


valores próprios A 


+1,731 


-1,731 


Um oscilador linear tem a equação diferencial 


2 
EO ao 


dt 


em que a é uma constante positiva. Escolhendo as variáveis de estado x1=y e x2=dy/dt, vem, 


X1 =X 


X2 =—4.X 


cujo plano de fase é 


semelhante ao anterior, que é um oscilador com a=3 (atente-se na matriz de estado). 
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vi) Ponto sela: quando um dos valores próprios é positivo e o outro é negativo. 


valores próprios: 1, -2. 


Chama-se ponto sela porque se cortarmos uma sela de cavalo transversalmente obtemos uma 
curva sem ponto de estacionaridade. Se a cortarmos transversalmente obtemos uma curva com 
um ponto de estacionaridade. Um ponto sela atrai a trajetória até certo ponto (efeito do valor 
próprio negativo, mas depois repele-a, efeito do valor próprio positivo). No ponto sela não 


passa nenhuma trajetória, embora apossamos aproximar tanto quanto quisermos, através das 


condições iniciais. 
5.1.2 Sistemas não lineares: várias singularidades possíveis e estabilidade local 


Exemplo: Modelo de Lotka-Volterra 


x, 2x (presas) x, à y (predadores) 

1 2 
dx,(t E 

- e ax (D)—bx (Dx,(t) x =x, —bx,x, 
O cx (e) + pit (O x =-04,+ pa, 


Calculando os pontos singulares, obtém-se 


0=ax (1) bx (x, (0) 
0=-ex, (1) + px, (x (1) 


| La-bx (Dx (t) = 04 x(t) 0 | [a-bx,(1)]=0 us 
> = EN 
[-c+ px(Dlo(D)=0nx(0)=0  [-c+px(D]=0 x(t)= 


jo cla 


p 
outro ponto singular: x(t)=0 e x(t)=0 


Linearizando as equações de estado obtém-se o Jacobiano 
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a-bl pe O =hE 
da of |a-bx  —bx E b pl. p 
Oo PX, — —C+pxX, ea a E a 
p'b = —c+p— o df) 
pb P, Pp p aa 


calculando os valores próprios de A, o seu polinómio característico é 


A bl 
AI-Al= Pi=gsphol=2?+ac 


a bp 
pe A 
Es 


cujas raízes são: 
AI-A|=052"+ac=0>1=0+ jvac 


Para uma dada população inicial de cada espécie, qual será a evolução delas? A singularidade 


é um centro (porquê?). 
Usando o PPLANES obtém-se as curvas de fase seguintes. 


presa'= (a-b predador) presa a=04 b=0.01 
predador ' = (p presa - c) predador c=0.3 p=0.005 


predador 


20 + — 


3 
T 
I 
i 
Ê 


T 
o Eg>—— 


São bem visíveis os dois pontos singulares: um ponto sela (a origem, valores próprios 0,4 e 
-0,3) e um centro (valores próprios +0,34641 e -0,34641). A mancha representa trajetórias 


muito próximas, já na bacia de atracão da origem. 
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Exemplo: O pêndulo rígido (Khoo, 236)) 


suporte, 
sem atrito, 
B=0 


L, comprimento 


m, 
massa total 


m.g 


torsão, nula 
K=0 


posição angular 


a 
m.g.L.senô 


Um pêndulo rígido é composto por um disco metálico pesado suspenso por uma corda 
(considerada sem peso) a um ponto de fixação vertical. Aplicando a segunda lei de Newton ao 
movimento do pêndulo na direção tangencial ao disco (sistema mecânico de rotação, J = m), 


obtém-se 


J.aceleração = binário 


2 
J o = -mgLsen0 , J=ml? é o momento de inércia 
2 
adia =- dia x senôQ = —k x send 
dt mL 


2 


dO uxseng=0 k=8 
dt L 


Pode-se implementar em Simulink e observar o seu comportamento. 


XY Graph 
0000 
oo 
Scope 
Signal Velocidade Posição teta 
Generator 
MATLAB 
rad2deg1 Function 
ipi MATLAB 
PPI2Pi| Function 
H—p 
Scope1 
MATLAB 
deg2rad Function 
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Quando o pêndulo inicia o movimento numa posição inferior a 180º com velocidade nula, 
verifica-se o movimento pendular normal (oscilatório). Se a velocidade inicial é não nula e se 
a posição inicial for próxima de 7 (por exemplo 0,997) o pêndulo gira sempre no mesmo 


sentido, passando pela posição superior e continuando. 


Escolhendo as variáveis de estado x1=0 e x2= qw teremos 


Xi =X 


x2 =-—k.senx, 


Com o PPLANES pode ver-se que existem muitas singularidades devido à posição angular 


cíclica com 27 (na figura k=1). 


As singularidades são de dois tipos: centros e pontos sela. Os centros correspondem ao 
movimento pendular e têm as coordenadas da posição vertical inferior; as selas correspondem 


ao movimento circular e têm as coordenadas da posição vertical superior. 


Exemplo com 4 singularidades 


x =(x1)-u y=» 


x =(x,-2)/-u 
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1) u=0 


0= (Gs dy aç=1 


0= (xs = 2) 5 &s =2 


Tem uma só singularidade, o ponto [1 2]!. O Jacobiano é aí a matriz nula. Há por isso dois 
valores próprios na origem e o sistema é instável (isto é, o sistema linearizado em torno da 
origem tem dois valores próprios na origem). Traçando as curvas de fase (figura seguinte), 
obtém-se um comportamento estranho: as trajetórias do 3º quadrante são atraídas para o ponto 
singular e depois repelidas para o 1º quadrante. As do 2º e 4º quadrantes comportam-se como 


se a origem fosse um ponto sela. 


x'=(x-1)2+u u=0 
y'=(y-2P+u 
7 T T T T T T T 
A I I I I I I I I 
LAVA AE HM SM AMU UU LIL LO A 
1 | 1 I 1 I I I I 
I Á I | I I I I ) 
I dá I I I I I I I 
Ei pod = cilada Pur SEO auacE Sc Pride Lap pd psi cu ii a CR il e ca E 
I I ) I I I I í I 
I I 4 I I I I I I 
BE ca dee a o o is e po MAS AZ DA 27 eo e 0] 
I I | I I Í I I I 
k I I I I 4 | I I 
! ! I I 1144 i I I I 
o a ia RR ERRORS ROGER PERES Fe | 
I I I t I ! I I I 
I | I ! I I I I I 
1EDTTS Je A HASTA IDA DOS A 2-+4=E À 
I I I I I I I I I 
I I I I I I | I I 
I I I I I I | ) I 
PSA POR TI TT AP RITI AA O TVA] 
I I I I I I I I I 
I ) I I I | I I I 
E ala Gac oa da o faca Ao lag Cega E eae] ua neo odio Ros ao ces dna 
I I I ] I | I I I 
t I | I t t | ! I 
4 0.5 0 05 1 15 2 25 3 
x 


(,12-u=0 SM 1) =uoS(x, -D=+u 
(4-2) -u=0 S(y-D=uS(y- D=+Vu 


o (x, -D=tVu OS x, 


II 
pa 
[ae 
[em 

[9] 


o (x-D=tJuo x =2t/u us=16> x =[2,3] 


x5 =[0,1]" 
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alo Safe Sao) ao] 


nó instável ponto sela ' ponto sela nó estável 


x'=(x-12-u u=1 
y'=(y-2P-u 


As curvas de fase estão traçadas 


na figura ao lado. 

Qual o estado alcançado, a partir 

de uma dada condição inicial ? 

O que acontece nas regiões de fronteira ? 


Note-se também que neste caso, não linear, 


a introdução de uma entrada exógena u altera 
qualitativamente o número e o tipo de singularidades. 

Qual a estabilidade do sistema não-linear ? 

Aqui temos que introduzir a noção de estabilidade local, diferente da estabilidade global. 
Qual a estabilidade local do sistema (na vizinhança de cada singularidade)? 


Depende da singularidade. Neste caso temos um nó estável, um nó instável e dois pontos sela. 


Graficamente é fácil de se concluir sobre a natureza de qualquer um deles. 


Para se estudar analiticamente da estabilidade de cada singularidade procede-se do seguinte 


modo: 


-Jineariza-se o sistema não linear em torno da singularidade (calculam-se os Jacobianos 


nas singularidades), 
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- calculam-se os valores próprios da matriz de estado obtida em cada caso. A 
estabilidade do sistema linear obtida é a estabilidade do sistema não-linear na vizinhança da 


singularidade (note-se que a linearização só é válida numa vizinhança da singularidade). 


Definição. Sistema localmente estável: um sistema diz-se localmente estável (na vizinhança 
de uma singularidade) se, afastado da singularidade por uma perturbação, a ela regressa por si 


próprio. 


Acontece quando a singularidade é um nó estável ou um foco estável. Não acontece no caso 
de um nó instável ou de um foco instável. Para um centro temos a condição limite de 
estabilidade, em que o sistema, excitado, oscila sem fim numa dada trajetória periódica 


(fechada). 


Sempre que temos um sistema não linear com várias singularidades, temos que aplicar a noção 


de estabilidade local a cada uma delas. Num sistema linear a estabilidade é global. 
5.1.3 A incerteza nas condições inicias e o caos 


Vejamos agora o que acontece na fronteira das regiões de convergência para uma dada 
singularidade. Por exemplo, no gráfico PPLANES, o segmento de reta y = 3, entre x = -2 ex = 


0, é uma fronteira. As trajetórias acima divergem para infinito, as abaixo convergem para a 


origem: 
x (estou u=1 wt=rodf£oy aid 
sea y'=(y-2P-u 
abl o. Lo Lo 1 III] T T T T T T T T T T T 
É i | ] ] Í I I | I I I Ec À rd Rae cu Rd da a 6 lu o Ea a E ps a pd ecc Er bi 1 
I I I I I I I I | I I I I I I I I I I I | i 
ad RR RS RR RN TR O aos 4-D p= Rook oOASSAs +++ HS 
I I I I I I I I I | I I I I I I I I I I | | 
qr nr enero qo qo Tera We 4j--Ch>rtoOrho=Os=dt=atTarTrar Ar 
I I I I I I I I I I I | I I I I I I I | | I 
eq TH STO TR TO TE DS TS TEITS aM-s-cr>o-> tornada, ar ps A-A+ 
I I I I I I I I I I I I I I I I I I I I I | 
Mer rr To TETE TS ISIS TST a02h4->-+>>DE>>D[>o 4 DD +04 
I I I I I I I I | b | I I I I I I I I I I | 
» VE SE ESSES cI ts et >» 3p4>D4>DDE>=DE>5D"D4 0 D5>DD ED Ho--- +- 
I I I I I I I I I | I I I I I I I I I I ! I 
DOR faça susoa fc aj ds Sinai asa faia ash fores o E us SS TS a egets-cdtecpecfoedcsqprlopodpoecpo AM 
I I I I I I I I I I I I I ! I I I I I I | I 
206 fd Po Agi dp lã plo Aptos patos prato dp tn age dotado doq SAL Ap Lp LApL o LA 
I I I I I I I I I I I I I I I I I I I I | i 
2Mh4->D4>D>D+>>DE> DDD DD 040 4>>+>5 +45 ess ILHA po DaLapo so Mgd 
I I I I I I I I I | I I I I I I I I I I | 
ag js Hs qo quo ds pasa as ao af abs aposto ER: 2 RARA RES A TR PR pp TO e UT ET Pp 
I I I I I I I I I | | I I I I I I I I I | ! 
pipe phosboçthesbosdss dass aspas elpa 28h 4-4 bobos d+ trate dj 
l l t 1 l t t t l l l t t t 1 L 1 L L l L | 
-2 -1.8 -1.6 -14 1.2 1 -0.8 0.6 -04 0.2 0 2 1.8 1.6 1.4 1.2 41 0.8 0.6 -0.4 0.2 o 
x x 


Reduzindo a escala do eixo vertical, reencontramos sempre uma fronteira bem definida. Na 


figura seguinte o intervalo vertical é [2.9999 3.0001]. 
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x'=(x-12-u u=1 
y'=(y-2P-u 
I I I I I I I I I I I 
3.0001 + — +—-+— H=>4 HDD DproDLp4 4a > S + 
| | | | | | | | | | 
I | | | | | | | | | 
3.0001 - — - — + — ipa de o Rd Cad Lada oi E ne a RU co AR IA CS a 
| | | | | | | | | | | 
| I I | I I I I I I I 
S000 pps pq psp SPrre Tr TNT are Prarpro 
| | | | | | | | | 
| ! I ! | ! I | | | | 
AGORA EAR TND EG a EG ER O O GGI AR TG PT 1 ARG ORE CR TA 
| I | | | | | | | | | 
Fe ER ON E ERA ORE LR SA DIA O A a DO ETA OLL oa LES E A UA rp 
I I I I I I I I | I 
| | | | | | | | | | | 
Ea É PERES SS E ee Ep 
I ! ! ! I I I | I | 
I | | | | | | I | | 
La E UNR UE PR REU POR E EE EE o ME OR APAE 
| | | | | | | | I | 
| | | | | | | | | | 
Busto tis fd id pl do piso po bo podia Pose poa A too] 
| | | | | | | | I | 
I | | | | | | | | I 
2.9999  — + fls obiaso ppbactlica ql ossasilis guobes coube qua dust bo grs actos pos galos aoilgps sl eoslpo o one) 
| | | | | | | | | | 
| | | | | | | | | | 
29999 —- = 7-4 Rr EITA E E RR SR ORE TOR Sar 
| | | | | | | | | ! 
| | I I I I I | I | | 
219008 us pop neoon oos oretas ao re eai cata recuds ond poa 
| | | | | | I ! I t | 
2 1.8 1.6 1.4 1.2 1 -0.8 0.6 0.4 0.2 0 
x 


Coloca-se agora a seguinte questão: se as condições iniciais são medidas por um instrumento, 
e se o erro do instrumento é superior a 0,0001, será possível prever a evolução deste sistema 


não-linear? 


Não é. Este facto ilustra uma das características do comportamento caótico de alguns sistemas 
não lineares: a impossibilidade de prever o seu futuro devido à incerteza nas condições 
iniciais. Uma diferença infinitesimal pode fazer a diferença. Daí o famoso efeito de borboleta 
primeiro utilizado por E.N. Lorenz em 1992 (Hilborn): “Predictability: Does the Flap of a 
Butterfly's Wings in Brazil set off a Tornado in Texas” que depois se repetiu em múltiplas 
versões geográficas “um bater de asas de borboleta em S. Francisco (i.e., uma pequeníssima 
variação das condições iniciais) pode provocar uma tempestade em Pequim (isto é, o sistema 
tende para valores elevados)”. Isto não quer dizer que o sistema seja estocástico. O sistema é 
de facto determinístico, mas a incerteza (da medição) nas condições iniciais impede-nos de 


prever a sua evolução. 
O futuro de um sistema caótico é indeterminado ainda que o sistema seja determinístico. 


As trajetórias que definem uma fronteira chamam-se separatrizes e separam o plano de fase 


em regiões de modos dinâmicos diferentes. 


O conjunto de pontos onde se originam trajetórias que levam a um certo ponto singular 
chama-se bacia de atração desse ponto (basin of attraction), por analogia com o conceito de 


bacia hidrográfica. 
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5.2. Ciclos limite e indução de frequências de oscilação 
Oscilador de van der Pol (Khoo) 


Em finais da década de 20 do séc. XX, van der Pol (e van der Mark) propuseram o pêndulo 
rígido como um modelo mecânico da atividade oscilatória do coração humano. Este oscilador 


é dado pela equação diferencial não linear de segunda ordem (Khoo, 240): 


Note-se que se c=0, obtém-se o oscilador linear. O termo em dx/dt exprime um amortecimento 


(atrito). Temos assim que o oscilador de van der Pol tem amortecimento não linear. 


A fim de se analisar as suas trajetórias de fase, aplique-se a transformação de Lienard, 


Es = a, SE 
Pede 3 de 3 
por outro lado 
3 ) 2 
ox x dy ldix ade de ldix 6 od 
cdt 3 dt cdt dt dt cdt dt 
I| d?x », dx 1 
= c(l—-x =—"—X 
e | E 


que poderemos estudar no plano de fase ou no Simulink. O autor segue aqui uma técnica 
interessante de definir as condições iniciais usando uma régua ajustável (slider gain, de Math 
Operations)) para as variar facilmente. As condições iniciais dos integradores são definidas 
como externas na janela de diálogo do bloco. Introduzindo uma excitação externa na equação 


de x, obtém-se o diagrama de blocos da figura. 
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+» >» x o D 1 MATLAB 
*os Function 
Produto Integrador X'3/3 
3 0.612 e] > 
Coefficient, C X0) Constant2 vs 
RR (amortecimento) 
MATLAB 
Function 
y x 
-1/C 
Produto1 
Integradort 
E] 2000 
oo 
XY Constante1 A) FÊ 
Curva de lase Excitação extema 
x 
> y 
To Workspace 


To Workspace1 


Oscilador de van der Pol com excitação externa (Khoo, 244). 


Com o PPLANES obtém-se o plano de fase seguinte, para c=3. 


x'=c(y-x/3+x) c=3 
y'=-x/e 
T T T T T T T T T T T 
I I I I I I I I I I I 
8 pe espe pe e th pe A o a Den ie AE 
I I I I I I I I i I I 
I I I I I I I I I I I 
4 Ee Ss so SAS Suas e agents, ga VE DOR sf 
I I I I I I I I I I I 
I I I I I I | I I I I 
l I I I I I I I | I I 
HC a q ST SS SG = SSI CNS EIS 
I I I I I I I I I I I 
I I I ] T I T I I I | 
APRE E SOR fo Re mf aj a ip a e DRE mi Re pg 
I I I I I I I VS I I I 
I I I | T I I nÍ | I I 
| E a e mina PD A E A ra 
I I I I I I I ”Y I I I 
| | | | id | | í | | I 
I I I [A I T I I I I 
Ee O-k=5="""[--D 0000 q/500poncag- o 5005/7000 qu = 7-4 
I I I | I I I LN I I I 
; I I ! o I I La I I I I 
Ah-4-D>0D0 0201-000 ME===p = qn eq pon 4 00 >> + —A 
I I I I I I I I I I I 
I I I | I I I t I I I 
E DO ERR Re A ir AA DR CR RD e e o a De or ra e pe DR 
I I I I I I I I I I ! 
I I I I I I I I I I I 
I I 1 I I I I I I I I 
TO | es gs sau ses ras aj a E A] 
I I I I | | ! I T I i 
I I I I I I I I I I I 
Apto += pt Gi o É jet | Sd pe a | | A sl 
I I I I I I I I I I I 
I I I I I I t I I I I 
I I t | I I I I I I I 
A RR O E RR RR 
[ l l Ii L l l [ l l L 
5 4 3 2 1 0 1 2 3 4 5 
x 


Todas as trajetórias vão dar à trajetória a grosso na figura. Tem-se assim, sempre, uma 
oscilação, e a trajetória respetiva é um ciclo limite estável (limit-cycle), pois é o limite de todas 


as trajetórias e elas aí permanecem. Mesmo se a condição inicial está dentro do ciclo limite, a 
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trajetória parte da condição inicial em direção ao ciclo limite e aí permanece. Se observarmos 
a evolução temporal de x, obtém-se (do Simulink), com a condição inicial x(0)=0,612 e 


y(0)-0,367, c=3, a curva, 


2.5 T T T T 


5 10 15 20 25 30 35 40 45 50 


onde se vê bem a oscilação característica de um ciclo limite. Estas oscilações são parecidas 
com a batida cardíaca. Têm período T=9 (calcula-se na resposta) e por isso a sua frequência 
natural (quando não há excitação externa, apenas condições iniciais não-nulas) fr=1/T=0,11 


Hertz. 


Se acrescentarmos uma excitação externa (entrada) sinusoidal na equação de y, 


dy = E = RS RD) ft) 
c 


dt cc 
e aplicarmos vários valores de f, obteremos os seguintes efeitos: 


excitação f=0,113Hz comportamento periódico 
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excitação f-0,01 Hz comportamento não periódico 


nona 0 Na 
T 


| 
0 


1 15 20 so E ao as 


1o 


As 20 30 35 ao as 


1o 15 20 30 Bs ao as 


so 
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nO ns 0 NA 
T 
] 


Verificamos uma propriedade interessante dos sistemas oscilatórios não lineares. 


Eles têm uma certa frequência natural, aquela que aparece quando a excitação externa é nula. 
Se excitarmos o sistema com uma frequência muito afastada (para mais ou para menos) dessa 
sua frequência natural, o comportamento resultante contém uma mistura de componentes que 
resultam da interação entre a frequência de excitação e a frequência natural; dessa mistura 


resulta em geral um sinal aperiódico. 


Mas à medida que a frequência de excitação se aproxima da frequência natural do oscilador 
não linear, a partir de certa altura o oscilador “adota” a frequência de excitação ou, dito de 


outro modo, a excitação induz uma frequência de oscilação do sistema. 


Esta propriedade é muito importante na prática: ela está na base dos pacemakers cardíacos. 
Também se verifica esta propriedade na adaptação dos ritmos biológicos ao ciclo dia-noite, o 
acoplamento entre a respiração e a pressão sanguínea, bem como a sincronização dos 


geradores fisiológicos de padrões centrais durante a caminhada e a corrida (Khoo, 244). 


Outros fenómenos fisiológicos podem ser modelizados com recurso a osciladores não lineares. 
Por exemplo o oscilador de Poincaré é apropriado para simular as arritmias cardíacas (ver 


Khoo, 246). 
5.3 Bifurcações e caos 


Considere-se o modelo de crescimento populacional de uma espécie (também chamada 


equação logística, Khoo, 274). Ele tem a forma que deduzimos no Cap. 2, em modo discreto, 
Xu = AxI-x,)= A(x, =) x e [0,1] 
em que A exprime a taxa de fecundidades da população. 


Vimos que tem dois estados de equilíbrio, 
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e a convergência para um ou para outro depende do estado inicial. Quando A<1 a única 


singularidade possível é a origem, xk=0 (a espécie extingue-se). 


Para o estudarmos com algum pormenor, dado que daí poderemos extrair muita informação, 


implemente-se em Simulink, usando o bloco de memória que faz o mesmo que um atraso 


puro: 
Populacao 
Introduzindo diversos valores de A e simulando, obtém-se o seguinte: 
A=2,8 A=3 
1 1 t T r t 
0.9+ 0.9+ ã 
0.8+ ã 0.8+ ã 
J o O.7+- 1 
É a Ea E tasas sis goi 
ê 0.6+ | É 06) 1 
0.5+ + 0.5+ ã 
0.4 + A 0.4 + ã 
035 2 4 6 8 10 12 14 16 18 20 ts 1000 2000 3000 2000 5000 6000 7000 8000 9000 1000 
Geraçao Geraçao 
A=2,8; atinge um regime final A=3; tende para um regime final 
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; | ne ess 
Desi 1 o.9- 
0.8» 1 0.8 
Fe) 0.7 + 1] fe) 0.7 
É 06! | ê 0.64 
0.5 + J 0.5 
0.4+ A 0.4 
0.3 1 IL IL 1 1 1 [ IL d 0.3 1 1 1 1 1 1 1 L L 
2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20 
Geracao Geraçao 
A=3,3,: oscila, período 1 geração A=3,5,: oscila, período 2 
(ciclo limite estável de período 1) (ciclo limite estável de período 2) 
A=3,55 A=3,6 
1 T T T 1 
0.9 0.9 | 
dn | | 1 08 | pol 
2 0.7 + 9 0.7L | 
ê 0.6 + ê 0.6 + | 
0.5 0.5 | 
0.4 + 4 04» J 
0.3 E L E [ E L L ; 5 0.3 1 l I 1 L a 1 L 1 [ 
Cr SSD vs O 2 4 6 8 10 12 14 16 18 20 
A=3,55,: oscila, período 8 Geração 
(ciclo limite estável de período 8) A=3,6: caos (não há período) 


O sistema apresenta um comportamento estranho: o seu período depende de A, a taxa de 
fecundidade. À medida que A aumenta, a população é periódica de período dependente de A: 


Assim A em 
13; 3,34495] oscila com período 2 
]3,4495; 3,54408] oscila com período 4 
]3,54408; 3,56440] oscila com período 8 
Faça-se agora a seguinte construção gráfica: 


Para cada valor de A calculam-se os valores da população x que compõem o respetivo ciclo. 


Estes valores calculam-se simulando no Simulink e lendo depois no Workspace os valores do 
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vetor popula que compõem um período. Por exemplo para A=3,3 teremos x=0,8236 ou 
x=0,4794; para A=3,5 será x e (0,3828; 0,5009; 0,8269; 0,8750 +. Para A=0,55, x e 40.3548, 
0.81265, 0.54049, 0.88168, 0.37034, 0.82781, 0.50601, 0.88737). E assim sucessivamente. 
Para valores de A inferiores a 3 só teremos um valor, o do regime permanente; para valores 
inferiores a 1, o regime permanente á a origem. Graficando obtém-se a figura seguinte, 


chamada o diagrama de bifurcação da função logística. 


Diagrama de bifurcação do modelo de crescimento populacional. 


-— 


Eua 
Do ooo oa a 
od mom 4 mam 


Ampliando a parte final, entre 3,4 e 4, obtém-se uma figura em que se repetem, a uma outra 
escala, os padrões da figura original. Se nesta segunda ampliássemos novamente, iríamos 
obter de novo os padrões da figura inicial. Isto é, aumentando o detalhe de observação 


encontram-se sempre os mesmos padrões. Esta é uma característica dos fractais. 


ME 
[um O O O 1 
[o E Rm = 
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(de http://www .vanderbilt.edu/AnS/psychology/cogsci/chaos/workshop/BD .html) 


En, 


TA 


Temos assim, com a variação de A, uma duplicação do período (de oscilação), até que deixa 
de se observar um comportamento periódico. Este comportamento, em que não há repetição de 
padrões, é caótico. Um ponto (no eixo de A) em que se iniciam dois caminhos chama-se 


bifurcação. Para se obter o diagrama de bifurcação procede-se do seguinte modo: 


Para cada valor de A calculam-se os valores da população x que compõem o respetivo ciclo. 
Por exemplo para A=3,3 teremos x=0,8236 ou x=0,4794; para A-3,5 será x e (0,3828; 
0,5009; 0,8269; 0,8750 ). Para A=0,55, x e (0.3548, 0.81265, 0.54049, 0.88168, 0.37034, 


0.82781, 0.50601, 0.88737). E assim sucessivamente. Graficando obtém-se as figuras acima. 
Os números de Feigenbaun 


Seja An o valor de A em que o período 2” dá origem ao perído 2"*!. O número 


= AA 
É A o A 
chama-se o delta de Feigenbaum. 
n 2" An Fo 
l 2 3 
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p. 4 3,449490 
3 8 3,544090 — 4,75147991543339 
4 16 3,564407 | 4,65619924201411 
5 32 3.568750 — 4,66842830882361 
6 64 3,569692 — 4,66452304394492 
7 128 3,569891 4,68844221104788 
8 256 3,569934 
9 512 3,569943 
10 1024 3,5699451 
1 2048 — 3,569945557 
mio É da 3,569945672 presa 
acumulação 


Nesta tabela (por exemplo em http://to-campos.planetaclix.pt/fractal/caos.html) o ponto de 


acumulação é o ponto em que deixa de haver período e se inicia o comportamento caótico. 


Feigenbaun estudou a evolução do seu delta para vários sistemas, no Laboratório Nacional de 


los Alamos (USA), e para todos encontrou em 1975 o seguinte resultado (Hilbomn 48): 


limó, = 4,66920161... 


n—500 


Segundo Hilborn (p. 48) este número está destinado a figurar no Panteão dos números da 


Física (tal como o mr, o número de ouro (1+/5 tas etc.) e tem sido considerado uma 


constante universal. 


5.4 Exemplos de sistemas fisiológicos com comportamento caótico 
Regulação da densidade de neutrófilos no sangue (de Khoo) 


O número de glóbulos brancos no sangue de um doente com leucemia mielóide crónica 
(chronic myeloid leukemia, CML) tem grande flutuações em torno de valores elevados. Essas 


flutuações são aproximadamente periódicas com um ciclo de 30 a 70 dias. Essas flutuações 
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foram atribuídas a alterações da dinâmica do sistema de controlo fisiológico que regula o 
equilíbrio entre a produção e a destruição dos neutrófilos (um dos tipos de glóbulos brancos) 
que circulam pelo corpo humano. Para modelizar esse sistema de regulação foi proposta a 


seguinte equação diferencial não linear com atraso puro (Khoo, 278) 


dx ; ; 
EN = taxa de lançamento de neutrófilos no sangue - taxa de morte dos neutrófilos 
t 


BO xt) 
O" +x(). 


Mas há um tempo T, de maturação, que demora a lançá-los no sangue, depois de produzidos. 


A taxa de produção na medula óssea segue uma lei do tipo 


Por isso teremos um atraso puro, ou de transporte, na equação diferencial. 
A taxa de morte dos neutrófilos é proporcional à sua concentração no sangue. 
Assim, 


dx BO'x(t-T,) DE <S side BO"x(t-T,) 
dt O" +x(t-T,) dt 9" +x(t-T,) 


x(t) £ densidade de neutrófilos (um dos tipos de glóbulos brancos) no sangue no instante t 
T, É atraso puro é o tempo de maturação (desde a produção até ao 
lançamento na corrente sauguínea) 
y É taxa de destruição das células devido a diversos fatores 
0 en determinam a relação entre a taxa de produção de neutrófilos e a 
densidade passada de neutrófilos 


é um fator de escala 


Sejam 
E = n=10, D=02. 4=01, 
Calculando os pontos singulares 
dx  0,2x(t-T,) 
dt 1+x(t-T)º 
como se trata de regime permanente, em que o atraso puro deixou de ter influência, 
0,2x(t) 

Ie x(t)? 


—0,Ix(t)=0 


-0,1xX(0) =0> 0,2X(0-0,IX(D-0, 10" =0> 0, 1X(Od-X(0")=0 


| xX(0)=0 Rn =0 
= 
I=x(0" =0 | x()=1 


conclui-se que tem dois estados de equilíbrio: a origem e o ponto 1 (a que corresponde uma 


taxa de produção de neutrófilos de 0,1. Implementando no Simulink, 


(DADC/FCTUC/2022 Licenciatura em Engenharia Biomédica/MCPF 241 


Dinâmica dos Sistemas Biológicos e Fisiológicos 


Cap. 5 Sistemas Não-Lineares e Caóticos 


To Workspace 


taxa de produção 


taxa de morte 


Beta 


MATLAB u? 
Produto Function Atraso puro 
1M(1+u) Potência 
de teta 


Atraso, Td 


Vejamos o que acontece para diversos valores de atraso puro (ou atraso de transporte) 


Densidade de neutrofilos (unidades normalizadas) 


Conccentraçao de neutrofilos 


0.8 


0.6 


0.4 


0.2 


Td=2 dias Td=4 dias 


11 T T T T T 


4 
Concentraçao de neutrofilos (unidades normalizadas) 
o 
o 
7 
1 


' ' ' ' ' ' : ' 1 ' 
100 200 300 400 500 60 0 100 200 300 400 500 60 
Tempo (dias) Tempo (dias) 


Td=6 dias Td=12 dias 
T t T t T 1.4 T T T T T 


0.6+ 4 


Concentração de neutrofilos 


0.4- 1 


' 1 ' : 4 1 ' : 4 1 
100 200 300 400 500 60C 0 100 200 300 400 500 600 
Tempo (dias) Tempo (dias) 
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Td=18 dias Td=25 dias 


0.8+ ! 0.8+ J 


0.6+ J 0.6+ E] 


Concentração de neutrofilos 


04 J 044 1 


Concentraçao de neutrofilos 


0.2+ d 0.2+ d 


(o) 1 00 200 Eno o 200 500 60C 0 100 200 sa Fra 400 500 60C 
Para pequenos valores de atraso de transporte, Ta , o sistema é estável em torno do regime 
permanente x=1. Nestas circunstâncias se por qualquer razão há uma variação brusca na 
densidade (por perda de sangue, por exemplo), o sistema de regulação faz com que seja 
rapidamente reposto o regime permanente. O atraso (que exprime o tempo de maturação dos 
neutrófilos) tem no entanto, à medida que assume valores significativos, um efeito 
desestabilizador, provocando oscilação aos 6 dias, oscilação que aumenta de amplitude com o 
atraso, até que surge um comportamento caótico para valores muito elevados de Ta. A 
oscilação com Tq=6 tem um período de 20 dias (aproximadamente), dentro do intervalo 


[17,28] que se observa nos humanos vítimas da doença chamada neutropenia cíclica. 


Para valores de Ta superiores a 20 dias, temos a situação dos doentes de leucemia crónica, 
expressa pelo comportamento caótico do sistema (flutuações aparentemente aleatórias, mas 


num sistema determinístico). 


Se experimentarmos com duas condições iniciais x(0) = 1,24 e x(0)=1,25 obtemos o 
comportamento típico de um sistema caótico: até 900 dias, as duas trajetórias pouco divergem, 
mas depois seguem percursos bem diferentes, mantendo-se ambas dentro de uma região 


limitada. 
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1.6 


O caos 


144 


1.2 + 


Concentraçao de neutrofilos 


l Í 
500 1000 1501 
Tempo (dias) 


Variabilidade cardiovascular (de Khoo) 
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No ciclo cardíaco verificam-se flutuações na duração do intervalo entre batidas, na pressão 


arterial, e na saída cardíaca. E possível interpretar essas flutuações como as de um sistema 


caótico, nomeadamente através do modelo de Cavalcanti e Belardinelli (1996) que estabelece 


que essas variabilidades espontâneas na batida cardíaca e na pressão sanguínea resultam do 


comportamento caótico que ocorre no sistema de controlo do barorreflexo. 


Poderemos representar este sistema de controlo pelo diagrama de blocos, admitindo que o 


volume de ejeção cardíaca é constante (versão simplificada do modelo 


Bellardinelli): 


Volume de ejeção 
constante Saída 
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Vs 


cardíaca 
Coração Mecânica 


Q=V4T circulatória P 
Pressão 
arterial 


Atraso 
Ta 


T 
Período 


cardíaco 
Atraso dos Barorrecetores Barorreflexo 


e do nódulo sinoauricular 
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Dados Vs e To coração produz a saída volumétrica Q que se sujeita à mecânica do sistema 
circulatório, saindo do coração a uma pressão P. As células sensíveis à pressão (barorreflexo) 
enviam ao cérebro o valor da pressão, e este envia ao nódulo sinoauricular (um ponto nos 
músculos do coração) o período T necessário para manter uma situação saudável. O 
barorreflexo a o nódulo sinoauricular levam algum tempo a reagir, que é representado pelo 


atraso puro Ta. Agindo sobre os músculos da aurícula, produz-se o período T. 


A mecânica circulatória, para uma dada saída (volume de sangue) cardíaca, determina o nível 
da pressão sanguínea resultante. Esta mecânica pode modelizar-se por um modelo Windkessel 
de 3 elementos (relembrem-se as analogias do Cap.2) que relaciona a pressão à saída do 


coração com o caudal sanguíneo. 


Arterial 
Left Compartment ; 
Ventricle Cardiac Peripheral Periphery 
—— Output Dana Flow pm 
. o e , Pá 
Ventricular E Volume Ed Venous 
Presse. dA, Prássure , i Pressure 
Dá R E : ad i a 
e Aorhic "oi E mm 
Valve 1 i Es 
Corgpliance Cónductance 
Ra í Re 


em que 


P pressão à saída do ventrículo esquerdo 

Pc pressão auxiliar para cálculo 

Ra resistência (fluídica) da válvula da aorta 

Rc resistência equivalente do sistema circulatório 

Cc complacência equivalente do sistema circulatório 
Q(t) fluxo sanguíneo que sai do ventrículo esquerdo 


g(t) fluxo sanguíneo que entra no sistema circulatório 


Ts é a duração da sístole e o ciclo cardíaco tem uma duração T. 
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O volume de ejeção (stroke volume) é a quantidade de sangue que sai durante um ciclo, isto é, 
T 
V= Í Q(Ddt 
0 


Durante a diástole a válvula da aorta está fechada, e por isso 
Q()=0, Ts<t< T 


Como Cce Rc estão em paralelo, a impedância equivalente é 


E 
7 - SG age 86 = R. 
Ri+ 1 sRC+1l sRC +1 
sC sC 
Agora 
P=(R,+2Z,)Q 
RRC 
EO psp Ro R+R+SRRC 
Q(s) sRC.+1 I+sRC. 


Ou, de outro modo, 
A(O = ql) + a(t) 
P-P=QR, 


1 t 
R=7](Q-ndt=Ra 


“BP P-RQ 

Rr, R, 
Tai 1 P-RQ 

P=— —qQdt+OR = “dt + QR 
e q)dt + QR, 1? E Jdt+ QR, 


Derivando ambos os lados da equação, 
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dt C R, “dt C R, dt 
RC, Ee QR,-P+RQ+RRC, ao 
dt dt 


dP 1 [o E her do 1 [CE-TARO), dQ 


e finalmente, 


RL. a P-RRLC. o +(R, + RJQ 
dt dt 


Aplicando a transformada de Laplace a esta equação diferencial obtém-se o mesmo resultado 


que pelo método da impedância equivalente: 


sR.CP(s)+P(s)=sRR.C.Q(s)+(R, + RJQ(s) 
(sR.C. +DP(s)=(sR R.C. +(R,+RJQ(s) 
P(s) sRRC.+H(R+R) 

Q(s) — SRC,+41 


Usam-se os valores C=1,333 ml mm Hg'!, Rc=0,900 mm Hg s ml"! R4=0,039 mm Hg s ml! 
Para completar o diagrama de blocos inicial é necessário: 

1) coração Q=Vs/T 

11) barorreflexo que dá a relação entre a pressão P e o período T 
Esta relação é dada por 


É a de 
T(P) =Toin fo uPiB, 
l+ye 
e exprime o efeito da perceção pelo sistema nervoso central do valor da pressão e do 


consequente ajuste do ritmo cardíaco. Toma-se Tmin=0,66 s Tmax= 1,2s Pe=89 mm Hg, q=31 


y=6,7.1053 


11) o atraso puro exprime o tempo de reação do sistema nervoso a alterações da 


pressão. 


Juntando todos os blocos, implementa-se no Simulink. O ficheiro pode ser descarregado de 


ftp://ftp iece.org/uploads/press/khoo, Cap 10, cvvarlimdl. Alguns blocos são subsistemas. 


Clicando neles podem-se ver em detalhe. A técnica dos subsistemas permite compor 
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diagramas complexos com maior legibilidade. Para criar um subsistema, selecionam-se os 
respetivos elementos e no menu “Edit” aciona-se “Create Subsystem”. O Simulink cria 


automaticamente as suas ligações de entrada (“In”) e de saída (“Out”) (de Khoo, 282) 


To Workspace Art Blood Pressure 


Transfer Fen 
(with initial states) 


Heart Rate 
HR Gain 
Variable 
T(t-Tau) Transport Delay 


Cardiac Output 
Q 


mts Tau 


Constant 
777 


Stroke Volume, Vs 
(Fixed) 


(nota: o T é o tempo de atraso Ta) 


Bloco subsistema T(P) Bloco subsistema Heart 


in 1 


Division Denominator 


Denominator (sub-subsistema) 
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— 6.7340M13 


gamma 


Product! 


Math — x 
Function 


Product e 
x 


Product2 171 


Division (sub-subsistema). 


(Dex tu (2) 


out 1 Product 1/T in 2 


“Heart Rate (HR)”, em batidas por minuto, é dado por 60/T, sendo T o período em segundos. 
Podem observar-se as evoluções temporais de P, de HR e as trajetórias no plano (Q,P), ou seja, 
as curvas de fase. O parâmetro mais influente é o atraso puro T. Por simulação pode ver-se que 
à medida que T aumenta , o sistema torna-se oscilatório e com período que cresce com Ta (ou 


seja Tau no diagrama Simulink). 


Ta-0,5 Plano de fase 
94 89 
92 + 5 88.98| 
90 + e 
; 88.96 
VA É , x 
ssr + E 
88.94 | E 
86, + 
88.92, RO 
sf + Fá 
el q 88.9, 
8o | É es.88- 
8h | 88.86 | 
| j j r ' j j j 
76 ! ! L ! ! L ! ! [ 
0 = O CABO 90 25 SO a dO 4 dO M4 MU5 UG UT MAB MU 95 95.1 95: 
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Ta-1 
100 1 ' 1 : 1 E 1 ' 1 100 
95 + 1 95 | q 
; 
90 | 7 90 | 5 
EE 1 es | 
8 1 so 1 
75 , L j : j : j : j E l 
(o) 5 10 15 2 25 3 3 40 45 5 bo 70 30 0 100 PET o 
Ta =2 
105 
105 RR 
100 + À 
100 É 1 
95 + 1 
95 + | 
90 + À E | 
85. J es | 
o 1 8o |. a 
75 ) 75 | + 
TO + 4 70 | 
65 L L 1 1 L L 1 1 L 65 1 1 1 L L 
0 5 10 15 20 25 30 35 40 45 50 60 70 80 90 100 110 120 


O modelo completo de Cavalcanti e Belardinelli (publicado em IEEE Trans. on Biomedical 
Engineering, vol. 43, nº 1, pp 982-989, 1996) inclui duas malhas de retroação que interagem 
uma com a outra: uma representa o efeito do barorreflexo na batida cardíaca (heart rate) tal 
como no modelo simplificado, e a outra representa o efeito do barorreflexo na contratilidade 
do coração, que por sua vez afeta o volume de ejeção (stroke volume) que aqui é variável e 
regulado pelo sistema nervoso. Da interação entre estas duas malhas de realimentação surge o 


caos, com veremos. Teremos o diagrama de blocos seguinte: 
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Agora o nódulo sinoarterial cumpre duas ordens emanadas do cérebro: uma para T e outra 


Saída 
cardíaca 
Q 
Coração Mecânica 
Q=V4T circulatória P 
Pressão 
arterial 


PÉ 
: Período | 
: cardíaco | 


i Atraso dos Barorrecetores Barorreflexo 


; e do nódulo sinoauricular 


Volume de: 


Atraso : ejeceção 


Ta 


i Atraso dos Barorrecetores | : Barorreflexo 
i e do nódulo sinoauricular i 


para a contratilidade (de que resulta o volume de ejeção Vs). 
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A variação de T com P, na primeira malha de retroação, mantém-se como no modelo anterior. 


É uma função de tipo sigmoidal (em forma de S). 


Quanto ao volume de ejeção Vs, ele é também uma função de tipo sigmoidal, mas com 


diferentes formas e parâmetros: 


l 
V (P) E ani x ES 
el 
H, 
S max = 86 
P,=25 


8 =72 (unidades apropriadas) 


v 


Implementando no Simulink, ou descarregando de 


(https://bcs.wiley.com/he-bes/Books?action=resource&bcsld=1 1 138&itemId=1119055334&resourceld=43992, 


27/07/2022) 
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Art Blood Pressure 


Cardiac Output 


Transfer Fen 
(with initial states) 


Math 
Function 


T(t-Tau) 


Variable 
Transport Delay 


Constant Delay, Tau 


Cardiac Output 
Q 


Vs Stroke Volume 


(Nota: o tempo de atraso Ta é o Tau do diagrama) 


Subsistema Vs(P) Sub-subsistema Bottom 
One| 1 | 
Vsmax na 
Pio o Ei 
Division in 1 out. Sum2 - Division 
Bottom 
[º pre 
Sub-subsistema Division 
in 1 
OE mes 
out 1 Product JT in 2 
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O tempo de atraso, em sistemas com retroação (feedback) é muito problemático: normalmente 


tende a desestabilizar o sistema. Vejamos a sua influência neste caso, por simulação repetida 


no Simulink para T=0,5;1:1,8:2,5. 


Ta-0,5: obtém-se regime estacionário (foco estável) 


100 


90 


85 


80 


75 
0 


92 


90 


ss 


86 


84 


82 


8o 


78 


76 
(o) 


Período cardíaco T (s) 


Plano de fase 


Saída Q (ml s'!) 
92 


Tempo (s) 


4 90 5 
a ss A 
1 ao: Fá | 
pa sat + 
* 82 A 
d so Ss J 


>. Pressão P (mm Hg) 
d 78 à J 


L 1 1 76 


r r r 1 
35 40 45 50 80 85 90 95 100 105 


obtém-se um ciclo-limite (centro) 


Q 
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100 
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100 


75 E E o A 


T os 
90 | 
85 + 
8o |- 
75H ; s 4 
70 ! r ! ! ! ! r ! ; 
o) 10 20 30 40 50 60 70 Bo 90 


Constata-se que para pequenos valores 


hn 
100 105 


de Ta o sistema alcança um regime estacionário, 


correspondendo a um foco estável. Aumentando Ta a dinâmica do sistema muda 


qualitativamente, aparecendo oscilações (Ta =1) e bifurcações (Ta =1,8) levando à duplicação 


do período de oscilação. Aumentando Ta ainda mais o sistema atinge o caos (Ta =2,5); pode-se 


observar no plano de fase que, neste caso, as trajetórias são todas distintas, não aparecendo 


repetições de padrões (no intervalo de tempo simulado, que generalizamos). A figura seguinte 


representa a trajetória do comportamento caótico no espaço tridimensional (HR- Heart Rate, 


Vs, P). Note-se o seu caráter errático. 


Plano de fase tridimensional (HR, Vs, P) 


100 


o 
q 


Cond. 
inicial 


(o) 
a) 


85 


80 . 


Pressao arterial P (mm Hg) 


75. 
90 


80 


70 
Volume de ejecçao Vs (ml) 60 


; 100 
90 
80 


: 70 
60 Ritmo cardicao HR (batidas/min) 


50 50 


(Esta figura obtém-se escrevendo na janela de comando do Matlab, após a simulação no Simulink, 


“>>plot3(HR,Vs,P”) 
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Estes dois exemplos de sistemas fisiológicos mostram que o caos pode aparecer em sistemas 
de ordem relativamente baixa quando existem tempos de atraso e malhas de retroação inter 
atuantes. Do ponto de vista fisiológico, dado que estes modelos são simplificações da 
complexa realidade, é natural que os resultados obtidos nem sempre estejam em concordância 


com a realidade. Mas eles permitem explicar muito do que realmente acontece. 


Estes exemplos mostram também a utilidade do que estudámos em todos os capítulos 
anteriores para o estudo dos sistemas fisiológicos: sistemas análogos, obtenção das equações 
diferenciais, função de transferência, espaço de estados, curvas de fase, etc., são temas 


essenciais para o estudo dos sistemas fisiológicos. 
5.5 Conclusão. 


Neste capítulo usou-se a representação no espaço de estados para estudar a dinâmica de 
sistemas fisiológicos, fortemente não lineares. O plano de fase (para sistemas de ordem 2) 
permite definir a estabilidade de cada ponto singular. O comportamento caótico, muito 


frequente em sistemas fisiológicos que contêm atrasos puros e retroação, está associado a 
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sucessivas bifurcações, que definem um caminho para o caos, produzidas por variações de 
parâmetros do sistema. O número de Feigenbaum, considerado uma constante universal, 


permite encontrar as bifurcações. 
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